Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,257,602 A→T N163I (AAT→ATT)  lexA → transcriptional repressor of SOS regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,257,6020AT100.0% 247.1 / NA 66N163I (AAT→ATT) lexAtranscriptional repressor of SOS regulon
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (35/31);  total (35/31)

CGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCC  >  NC_000913/4257454‑4257746
                                                                                                                                                    |                                                                                                                                                
cgcgTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTa                                                                                                                                                <  2:1062019/151‑1 (MQ=255)
   gTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTaa                                                                                                                                               >  1:429471/1‑149 (MQ=255)
     cAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTc                                                                                                                                           <  2:577706/151‑1 (MQ=255)
      aGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCg                                                                                                                                          >  2:218082/1‑151 (MQ=255)
         gggATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGaa                                                                                                                                        >  2:1087269/1‑150 (MQ=255)
         gggATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAAc                                                                                                                                       >  2:1070388/1‑151 (MQ=255)
          ggATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACt                                                                                                                                      <  1:1028325/151‑1 (MQ=255)
          ggATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACt                                                                                                                                      >  2:1208068/1‑151 (MQ=255)
           gATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTg                                                                                                                                     >  1:829736/1‑151 (MQ=255)
              gTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTg                                                                                                                                  >  1:623823/1‑151 (MQ=255)
                cGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGcc                                                                                                                                >  2:1034168/1‑151 (MQ=255)
                     aaaGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGaaa                                                                                                                           <  2:426364/151‑1 (MQ=255)
                      aaGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGaaa                                                                                                                           >  1:682703/1‑150 (MQ=255)
                       aGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAAt                                                                                                                         <  1:60467/151‑1 (MQ=255)
                         atatCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAg                                                                                                                       <  2:142149/151‑1 (MQ=255)
                          tatCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGc                                                                                                                      >  1:29515/1‑151 (MQ=255)
                          tatCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGc                                                                                                                      >  1:495397/1‑151 (MQ=255)
                           atCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGc                                                                                                                      <  2:37707/150‑1 (MQ=255)
                             cGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGa                                                                                                                    <  1:257369/150‑1 (MQ=255)
                              ggCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGt                                                                                                                  <  1:1070257/151‑1 (MQ=255)
                                cATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGttt                                                                                                                >  1:11101/1‑151 (MQ=255)
                                cATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGttt                                                                                                                <  1:293454/151‑1 (MQ=255)
                                cATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGttt                                                                                                                <  2:640183/151‑1 (MQ=255)
                                cATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGtt                                                                                                                 <  2:165324/150‑1 (MQ=255)
                                  ttATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTa                                                                                                               >  1:1234176/1‑150 (MQ=255)
                                       gatggTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCa                                                                                                          >  2:750418/1‑150 (MQ=255)
                                              aCTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGt                                                                                                  >  1:21501/1‑151 (MQ=255)
                                              aCTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGt                                                                                                  >  2:1167041/1‑151 (MQ=255)
                                                 tGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTg                                                                                                <  2:651670/150‑1 (MQ=255)
                                                                aaaCTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCagag                                                                                <  1:178838/151‑1 (MQ=255)
                                                                     cAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTc                                                                            >  1:642603/1‑150 (MQ=255)
                                                                      cggATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCAc                                                                          >  2:118498/2‑151 (MQ=255)
                                                                           gTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCAtt                                                                      <  2:1241310/150‑1 (MQ=255)
                                                                           gTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTg                                                                     >  2:984576/1‑151 (MQ=255)
                                                                           gTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTg                                                                     <  1:484136/151‑1 (MQ=255)
                                                                               gTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAgg                                                                 <  1:395107/151‑1 (MQ=255)
                                                                               gTAACGGTCAGGACGTTGTCGCACGAATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAgg                                                                 >  2:612667/1‑151 (MQ=255)
                                                                                   cGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCt                                                              <  1:929825/150‑1 (MQ=255)
                                                                                     gTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGc                                                           <  1:81233/151‑1 (MQ=255)
                                                                                     gTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGc                                                           >  2:728006/1‑151 (MQ=255)
                                                                                      tCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGc                                                           >  2:475363/1‑150 (MQ=255)
                                                                                          gTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTg                                                      <  2:1029851/151‑1 (MQ=255)
                                                                                               tGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGt                                                  <  2:770449/150‑1 (MQ=255)
                                                                                                 tCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGtt                                                 >  1:132163/1‑149 (MQ=255)
                                                                                                     aCGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTAtt                                              >  2:426955/1‑148 (MQ=255)
                                                                                                         aTTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACg                                       <  1:1249589/151‑1 (MQ=255)
                                                                                                           tgatgaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACgg                                      <  2:408699/150‑1 (MQ=255)
                                                                                                            gatgaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCg                                    >  1:660942/1‑151 (MQ=255)
                                                                                                               gaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTTGCAAccggga                                   >  2:900647/1‑144 (MQ=255)
                                                                                                                aCGAAGTTATCGTTAAGCGCCTGAAAAATCAGGGCATTAAAGTCGAACTGTGGCCGGAAAATAGCGAGTTTGAACCAATTGTCGTTGACTTTCGTCAGCAGAGCTTGACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGGGACt                                 <  2:955520/150‑1 (MQ=255)
                                                                                                                 cGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGGTGGGGGTTTTGGCAACGGCggcg                                 >  2:1039564/1‑146 (MQ=255)
                                                                                                                  gAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGGGTTATTTGCAACGGCGACt                                 >  2:559104/1‑148 (MQ=255)
                                                                                                                     gTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGt                           <  2:829706/151‑1 (MQ=255)
                                                                                                                      ttACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTTGCAACGGGGACTGGCTGt                           >  2:852375/1‑150 (MQ=255)
                                                                                                                       tACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGGGGTATTCGCAACGGCGACCgg                               >  2:203143/1‑145 (MQ=255)
                                                                                                                          cGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTaaaa                       <  1:517243/150‑3 (MQ=255)
                                                                                                                           gTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACata                     <  2:943653/151‑1 (MQ=255)
                                                                                                                             tAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACatat                    <  2:1186707/150‑1 (MQ=255)
                                                                                                                             tAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTTGCAACGGCGACTGGCTGTAACata                     >  1:395698/1‑149 (MQ=255)
                                                                                                                               aGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTAACatat                    >  2:508869/1‑148 (MQ=255)
                                                                                                                                 cgcCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTg               <  1:367586/151‑1 (MQ=255)
                                                                                                                                    cTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGGGTTATTTGCAACGGCGACTGGGTGTAAAATATCTCTggg             >  2:287499/1‑148 (MQ=255)
                                                                                                                                       aaaaaaCAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAAAGGCGACTGGCTGTAACATATCTCTGAGGcc          >  1:1142923/1‑150 (MQ=255)
                                                                                                                                        aaaaaCAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACcgc        <  1:628741/151‑1 (MQ=255)
                                                                                                                                             cAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCCCAACGGCGACCGGCTGTAAAATATTTCTCAGACCCCGGTg    >  2:523546/1‑150 (MQ=255)
                                                                                                                                              aGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATgcc  <  1:695213/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCC  >  NC_000913/4257454‑4257746

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: