Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,944,240 T→G 5.8% intergenic (+41/+106) znuB → / ← ruvB zinc ABC transporter permease/ATP‑dependent DNA helicase, component of RuvABC resolvasome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,944,2400TG 5.8% 468.5 / 13.8 133intergenic (+41/+106)znuB/ruvBzinc ABC transporter permease/ATP‑dependent DNA helicase, component of RuvABC resolvasome
Reads supporting (aligned to +/- strand):  ref base T (72/53);  new base G (3/5);  total (75/58)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 136 total aligned reads displayed.

CGGTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGCC  >  NC_000913/1944094‑1944388
                                                                                                                                                  |                                                                                                                                                    
cGGTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTa                                                                                                                                                    <  1:6207/149‑1 (MQ=255)
  gTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                                  <  3:774483/149‑1 (MQ=255)
   tttAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                                 <  1:992955/149‑1 (MQ=255)
   gatctccTTTTCCGTGGTTTACTATACGCCGGCGGGTCCTTCGGTGGTGCTATGTGCGGCAATGTTATTTATTCTCAGTTTGATGAAAAAGCAGGCCAGCTAATCTACCCCCGAACACATTTGTCGGATGCGGCGCGAGCGCCGAATcc                                                                                                                                                 <  2:472088/144‑1 (MQ=255)
     tAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGa                                                                                                                                               <  3:739491/149‑1 (MQ=255)
      aaCCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGAc                                                                                                                                              >  1:576136/1‑149 (MQ=255)
         cTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGagatc                                                                                                                                           >  1:472088/1‑145 (MQ=255)
         cTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTa                                                                                                                                           >  2:364057/1‑149 (MQ=255)
         cTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTa                                                                                                                                           >  3:434359/1‑149 (MQ=255)
         cTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTa                                                                                                                                           >  2:706555/1‑149 (MQ=255)
         cTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTa                                                                                                                                           >  4:183008/1‑149 (MQ=255)
         cTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATAAGGCGCTCGCGCCGCATCCGAcaaa                                                                                                                                           <  2:200433/149‑4 (MQ=255)
          ttttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTAc                                                                                                                                          >  2:1099603/1‑149 (MQ=255)
          ttttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTAc                                                                                                                                          >  3:288740/1‑149 (MQ=255)
          ttttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTAc                                                                                                                                          >  2:715861/1‑149 (MQ=255)
          ttttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATAAGGCGCTCGCGCCGCATCCGAcaaat                                                                                                                                          <  3:1430046/149‑5 (MQ=255)
          ttttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATAAGGCGCTCGCGCCGCATCCGAcaaat                                                                                                                                          <  4:1053524/149‑5 (MQ=255)
              ccGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGtt                                                                                                                                      <  4:702572/149‑1 (MQ=255)
                    tttACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTAt                                                                                                                                <  1:939253/149‑1 (MQ=255)
                    tttACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTAt                                                                                                                                >  3:987629/1‑149 (MQ=255)
                       aCGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTAtctc                                                                                                                             <  4:1387621/149‑1 (MQ=255)
                         gATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTg                                                                                                                           >  2:1303304/1‑149 (MQ=255)
                          aTACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGAATCTCTgg                                                                                                                          <  4:127897/149‑1 (MQ=255)
                            aCGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTa                                                                                                                        <  3:1323528/149‑1 (MQ=255)
                                 ggcggGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCt                                                                                                                   <  3:1039327/149‑1 (MQ=255)
                                    gggTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGt                                                                                                                <  3:188876/149‑1 (MQ=255)
                                         cGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGAc                                                                                                           <  3:68120/149‑1 (MQ=255)
                                           tCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGAcgc                                                                                                         <  4:1451282/149‑1 (MQ=255)
                                                   gCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCg                                                                                                 >  1:1125021/1‑149 (MQ=255)
                                                     tATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTGAGACGCGAACAGAGTc                                                                                               <  4:432592/149‑1 (MQ=255)
                                                     tATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTc                                                                                               <  1:1020638/149‑1 (MQ=255)
                                                     tATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTc                                                                                               <  2:1312014/149‑1 (MQ=255)
                                                       tgtgCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGc                                                                                             <  4:514453/149‑1 (MQ=255)
                                                          gCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATc                                                                                          <  3:520507/149‑1 (MQ=255)
                                                            ggCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAg                                                                                        <  4:444569/149‑1 (MQ=255)
                                                            ggCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAg                                                                                        <  2:112260/149‑1 (MQ=255)
                                                                 tGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGGata                                                                                   >  2:38141/1‑149 (MQ=255)
                                                                 tGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCata                                                                                   >  4:1239739/1‑149 (MQ=255)
                                                                  gTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCatat                                                                                  <  1:706555/149‑1 (MQ=255)
                                                                  gTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGCCGCATCCGGCactg                                                                                  >  4:403028/1‑146 (MQ=255)
                                                                   ttatttatTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTa                                                                                                                        <  3:985454/110‑1 (MQ=255)
                                                                   ttatttatTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATAtt                                                                                 >  4:237948/1‑149 (MQ=255)
                                                                   ttatttatTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATAtt                                                                                 >  3:878113/1‑149 (MQ=255)
                                                                   ttatttatTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATAtt                                                                                 >  1:885260/1‑149 (MQ=255)
                                                                   ttatttatTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATAtt                                                                                 >  3:1047682/1‑149 (MQ=255)
                                                                     atttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTc                                                                                               >  2:326762/1‑133 (MQ=255)
                                                                     atttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTc                                                                                               <  1:326762/133‑1 (MQ=255)
                                                                      tttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGcc                                                                              >  4:1162117/1‑149 (MQ=255)
                                                                      tttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGcc                                                                              >  2:1434021/1‑149 (MQ=255)
                                                                      tttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGcc                                                                              >  1:1199228/1‑149 (MQ=255)
                                                                          ttCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTg                                                                          >  3:979973/1‑149 (MQ=255)
                                                                          ttCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTg                                                                          >  1:1475423/1‑149 (MQ=255)
                                                                            ctcAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGcc                                                                        <  2:683384/149‑1 (MQ=255)
                                                                            ctcAGTATGATGAAAAAGCAGGCCAGCTAATCTATCCCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACTCGAACAGCGTCGCATCAGGCATATTGCCAGTGcc                                                                        <  3:1292196/149‑1 (MQ=255)
                                                                             tcAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCg                                                                       >  3:588740/1‑149 (MQ=255)
                                                                               aGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGa                                                                     <  1:534910/149‑1 (MQ=255)
                                                                                     atgaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATgcggc                                                               >  3:328818/1‑149 (MQ=255)
                                                                                      tgaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATgcggcg                                                              >  4:1138097/1‑149 (MQ=255)
                                                                                       gaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTAt                                                                                                                                >  4:295838/1‑82 (MQ=255)
                                                                                       gaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGgcgc                                                             >  4:375449/1‑149 (MQ=255)
                                                                                       gaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGAAGGGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTAGTAGTCCTGGTCAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGgcgc                                                             >  4:897197/1‑149 (MQ=255)
                                                                                        aaaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGcgcg                                                            >  4:1330804/1‑149 (MQ=255)
                                                                                           aaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACTGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCgcgagc                                                         >  4:1317378/1‑149 (MQ=255)
                                                                                            aGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCgcgagcg                                                        <  1:222111/149‑1 (MQ=255)
                                                                                                 gCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGACATATTGCCAGTGCCGGATGCGGCGCGAGCGACCaa                                                   <  2:576136/149‑1 (MQ=255)
                                                                                                   cAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATc                                                 <  4:1053763/149‑1 (MQ=255)
                                                                                                   cAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATc                                                 <  1:43560/149‑1 (MQ=255)
                                                                                                    aGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATcc                                                >  1:1070544/1‑149 (MQ=255)
                                                                                                            tATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACg                                        >  3:250316/1‑149 (MQ=255)
                                                                                                             aTCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTgg                                                                                                                          >  3:335011/1‑66 (MQ=255)
                                                                                                             aTCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTgg                                                                                                                          <  4:335011/66‑1 (MQ=255)
                                                                                                             aTCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACgg                                       >  3:588742/1‑149 (MQ=255)
                                                                                                              tCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGc                                      >  1:956483/1‑149 (MQ=255)
                                                                                                                gCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTAc                                         <  3:1328084/144‑1 (MQ=255)
                                                                                                                  tGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCAttt                                  >  4:941287/1‑149 (MQ=255)
                                                                                                                       acaTTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggc                             <  3:183008/149‑1 (MQ=255)
                                                                                                                         gtaggTCGGATAAGGCGCTCGCGCCGCATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTc                                 <  3:1045096/139‑1 (MQ=38)
                                                                                                                         aTTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggcgg                           >  3:229712/1‑149 (MQ=255)
                                                                                                                         aTTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggcgg                           >  4:788915/1‑149 (MQ=255)
                                                                                                                         aTTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggcgg                           >  1:1310212/1‑149 (MQ=255)
                                                                                                                         aTTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGGGCCGGATGCGGCGCGCGCGACCACTCCGACTTACGGCATTTCTggcgg                           >  2:1026969/1‑149 (MQ=255)
                                                                                                                         aTTTGTCGGATGCGGCGCGAGCGCCTTAACCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggcgg                           >  4:1329300/1‑149 (MQ=255)
                                                                                                                          tttGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggcggc                          >  1:336225/1‑149 (MQ=255)
                                                                                                                               cGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTAt                     <  1:1283893/149‑1 (MQ=255)
                                                                                                                                 gataaggcgCTCGCGCCGCATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGc                   >  4:810018/6‑149 (MQ=39)
                                                                                                                                 gataaggcgCTCGCGCCGCATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGc                   >  4:399096/6‑149 (MQ=39)
                                                                                                                                  aTGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGcc                  >  2:1227226/1‑149 (MQ=255)
                                                                                                                                             gcgcCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTcc       <  2:942935/149‑1 (MQ=255)
                                                                                                                                                 cTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCgcgagcga                                                       >  3:1181753/1‑97 (MQ=255)
                                                                                                                                                 cTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCgcgagcga                                                       <  4:1181753/97‑1 (MQ=255)
                                                                                                                                                 cTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGc   >  1:353439/1‑149 (MQ=255)
                                                                                                                                                 cTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGc   >  2:305395/1‑149 (MQ=255)
                                                                                                                                                 cTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGc   >  2:1373971/1‑149 (MQ=255)
                                                                                                                                                 cTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGc   >  2:1423446/1‑149 (MQ=255)
                                                                                                                                                 cTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGc   >  3:1193930/1‑149 (MQ=255)
                                                                                                                                                 cTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGc   >  3:1214577/1‑149 (MQ=255)
                                                                                                                                                 cTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCCGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGc   >  3:339483/1‑149 (MQ=255)
                                                                                                                                                  ttATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGcc  >  1:173286/1‑149 (MQ=255)
                                                                                                                                                  ttATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGcc  >  2:576042/1‑149 (MQ=255)
                                                                                                                                                  ttATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGcc  >  2:902154/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGGTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATTCCACGCC  >  NC_000913/1944094‑1944388

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: