Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,226 T→G 12.3% intergenic (‑369/+126) ybjE ← / ← aqpZ putative transporter/aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2260TG12.3% 260.4 / 6.0 133intergenic (‑369/+126)ybjE/aqpZputative transporter/aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (69/34);  new base G (4/15);  total (82/51)
Fisher's exact test for biased strand distribution p-value = 2.56e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

CATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGG  >  NC_000913/915077‑915373
                                                                                                                                                     |                                                                                                                                                      
caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGACGGCGATTGGCCGCCCCCt                                                                                                                                                        >  1:4057599/1‑150 (MQ=17)
 aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGACGGCGATTGGCCGCCCCCT‑g                                                                                                                                                      <  2:4057599/150‑1 (MQ=17)
    gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑gc‑gc                                                                                                                                                  <  2:2861275/150‑1 (MQ=32)
        tgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCT‑GC‑GCGCTc                                                                                                                                              >  1:2520259/1‑150 (MQ=17)
          ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCA‑GC‑GCGCGCCt                                                                                                                                            >  1:2469663/1‑150 (MQ=12)
          ggCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCt                                                                                                                                            >  2:3629131/1‑150 (MQ=32)
           gCTGAGGGCGCGTTTTACAGGTATGTTACCTCGCGCCCGACCGGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCGGGGGGCGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCtt                                                                                                                                           <  2:76603/150‑1 (MQ=255)
             tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTgc                                                                                                                                         >  2:1162011/1‑150 (MQ=32)
                    gcgTTTTACAGGGACGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTg                                                                                                                                  >  2:757118/1‑150 (MQ=17)
                     cgTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCCGCGGCGTTTTTGTCCGCCAGAGCCGGGGGGTGCCGGGGGCGGCGGCGTTTGGCCGCCCCCG‑GC‑GCGCTCCTTGCGCCAGTgg                                                                                                                                 <  2:2428483/150‑1 (MQ=11)
                      gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGc                                                                                                                                >  1:1711562/1‑150 (MQ=33)
                       ttttACTGGGTTGTTACCTCGCGCCCTACGGGTTTGCCGTAGGGTATATGTGGTGGGCCCCGCGCCGCGGGTTTTTTTTTCGCCAGAGCGGGGGGGGGCTGGGGGGGGCGGCTTTGGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCa                                                                                                                               <  2:3337805/150‑1 (MQ=2)
                          tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAata                                                                                                                            >  1:4010060/1‑150 (MQ=37)
                           aCAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCG‑GC‑GCGCTCCGTGCGCCCGGGGCAaaat                                                                                                                           >  1:2861589/1‑150 (MQ=11)
                             aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGGTTGGCCGCCCCCA‑GC‑GCGCGCCATGCGCCCGGGGCAAAATG‑t                                                                                                                        >  1:356473/1‑150 (MQ=11)
                             aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑t                                                                                                                        >  2:3398667/1‑150 (MQ=39)
                             aGGGAAGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATAAG‑t                                                                                                                        >  2:3397925/1‑150 (MQ=255)
                               ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTg                                                                                                                      >  2:1707651/1‑150 (MQ=40)
                                gATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGc                                                                                                                     >  2:1364582/1‑150 (MQ=255)
                                 ttGTTCCCTCGCGCCCGCCCGGTTTGCGGTATGGTTTATGTCGAGGGCACGGCGCTGGGGCTTTTTTGTTCGCCTGAGCCGGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCt                                                                                                                    <  2:684297/149‑1 (MQ=14)
                                   gTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTa                                                                                                                  >  1:3119019/1‑150 (MQ=255)
                                      aCCTGGCGCCCGACCCGGTTGCCGTAAGGGATACGTGGTGGGCACCGCGCCGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCTGGGGGCGGCGGCGTTTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCt                                                                                                               <  2:2767948/150‑1 (MQ=14)
                                         tcccccccGACCCGTTAGCCTTAAGGTATAAGCCGAGGGCACCGCCCAGCGGCTATTTTGTTCGCCAGAGCCCGGGGGTCCAGGGGGCGGCGGCCATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCAt                                                                                                            <  2:2760146/145‑1 (MQ=14)
                                         tCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGCCGCCCCCA‑GC‑CCGCGCCTTGCGCCAGAGGCAAAATG‑TTGCTTATCTCAt                                                                                                            >  1:21263/1‑150 (MQ=14)
                                            cgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATT‑TTTCGTAGCTCATGaa                                                                                                         >  1:3164258/1‑150 (MQ=255)
                                              cccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCA‑GC‑GCGCTCCTTGCGCCAGTGGCAATATT‑TTGCTTAGCTCATGTAAg                                                                                                       >  1:2450508/1‑150 (MQ=255)
                                               ccGACCCGGTAGCCGTAAGGGATAATTCTAGGGCACCGCGCAGCTGAGATTTTGTTAGCCAGAGCCCGGGGGTGCATGGGGCGGCGGGGATTTGGCGCCCCCT‑TC‑GCGCTCATTTCGCCAGTGTCAATATG‑TTTATTAGCTAATTaaatg                                                                                                      >  1:1704562/1‑148 (MQ=14)
                                                      ggTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGcaa                                                                                               >  2:2334941/1‑150 (MQ=255)
                                                          gCCTTATGGTATATGTCGTGGGCACCGCGCTGCGGCGTTTTTGTTCGCCTGAGCCGGGGGGTGCTGGGGGCGGGGGGGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAg                                                                                           <  2:4100307/150‑1 (MQ=18)
                                                           ccGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGGCCGCCCCCA‑GC‑GCGCTCCATGCGCCAGTGGCAATATG‑TTGCGTATCTCATGAAAGGAGCGCAACAAga                                                                                          >  1:2176570/1‑150 (MQ=17)
                                                            cGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCACGGGGCGGCGGCGATTGGCCGCCCCCG‑GC‑GCGCACCATGCGCCACTGGCAATATT‑TTTCTTATCTCATGAAAGGAGCGCAACAAaat                                                                                         >  1:3741944/1‑150 (MQ=17)
                                                              tacgggATCCGTCGCGGGCCCCGCGCTGCGGCGTTTTTGTTGGCCTGCGCCGGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                                       <  2:2388362/147‑1 (MQ=18)
                                                              tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                                       >  2:617873/1‑150 (MQ=255)
                                                              tAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCA‑GC‑GCCCTCCGTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACCAgatga                                                                                       >  2:1416999/1‑150 (MQ=255)
                                                                  ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGaa                                                                                   >  2:2327313/1‑150 (MQ=255)
                                                                   gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCCTGAAAGGAGCGCAACAAGATGATGaaa                                                                                  >  2:333390/1‑149 (MQ=255)
                                                                   gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                                  >  2:4020883/1‑150 (MQ=255)
                                                                   gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAt                                                                                  >  2:1502339/1‑150 (MQ=255)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                            >  2:1713713/1‑150 (MQ=255)
                                                                         tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                            >  2:4083084/1‑150 (MQ=255)
                                                                           gAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGGCCGCCCCCA‑GC‑GCGCGCCGTGCGCCCCGGGCAATAAT‑TTTCTTATATCATGTAATGAGCGCAACAACAATATGAATCaaaaaaa                                                                          >  1:1319427/1‑145 (MQ=14)
                                                                            aGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGATTGGGCGCCCCCA‑GC‑GCCCGCCGTGCGCCCCTGGCAATATGTTTGGTTTG‑TCATTAAAGGAGCGCCACAAGAAGTTGAATCAAAAtata                                                                         >  1:387917/1‑150 (MQ=14)
                                                                             gggCACCGCGCAGCGGCGTTTTTGTTCGCCAGCGCCGGGGGGTGCCGGGGGGGGCGGCGTTTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                                        <  2:318495/150‑1 (MQ=21)
                                                                                cACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaa                                                                     >  2:2933319/1‑150 (MQ=255)
                                                                                   cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAacaacaa                                                                  >  2:2797807/1‑149 (MQ=255)
                                                                                    gcgcAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATc                                                                 >  2:91635/1‑150 (MQ=255)
                                                                                         gcggctATTGTGTCCGCCAGAGCCCGGGGGGGCTGGGGGCGGCGGCGATGGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaa                                                            <  2:3782303/150‑1 (MQ=25)
                                                                                           ggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaa                                                          >  2:2307172/1‑150 (MQ=255)
                                                                                              gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAAAAACATCTTAAAAAAAg                                                       >  2:1154966/1‑150 (MQ=255)
                                                                                                 tttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGTcc                                                    >  2:3782820/1‑150 (MQ=255)
                                                                                                    gTTCGCCTTCGCCGGGGGGTGCTGGGGGCGGGGGCGTTTGGCGGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGa                                                 <  2:1298794/150‑1 (MQ=37)
                                                                                                    gTTCGCCAGAGCCGGGGGGTGCTGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGa                                                 <  2:2725772/150‑1 (MQ=255)
                                                                                                     ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGAc                                                >  2:2773973/1‑150 (MQ=255)
                                                                                                     ttCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGAc                                                >  2:20917/1‑150 (MQ=255)
                                                                                                           agagCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACg                                          >  2:2599152/1‑150 (MQ=255)
                                                                                                               ccggggggTGCTGGGGGGGGCGGCGTTTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                                      <  2:3044484/147‑1 (MQ=255)
                                                                                                               ccggggggTGCGGGGGGGGGCGGTGTTGGGCCCCCCCCG‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAg                                      <  1:2377072/147‑1 (MQ=38)
                                                                                                                cggggggTGGTGGGGGCGGGGGCGTTTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                                     <  2:1711562/148‑1 (MQ=255)
                                                                                                                   ggggTGCTGGGGGGGGGGGCGTTGGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCt                                  <  2:3013409/150‑1 (MQ=38)
                                                                                                                    gggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                                 >  2:2095279/1‑150 (MQ=255)
                                                                                                                    gggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt                                 >  2:1562714/1‑150 (MQ=255)
                                                                                                                       tGCTGGGGGGGGGGGCTTTTGGCCGCCCCCG‑CC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                              <  2:564680/150‑1 (MQ=255)
                                                                                                                       tGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTTCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTg                              >  1:3696089/1‑150 (MQ=255)
                                                                                                                         ctgggggcgggggcGATGGGCCGCCCCCC‑CC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                            <  1:3339276/148‑1 (MQ=255)
                                                                                                                         cAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                            >  2:3460602/1‑150 (MQ=255)
                                                                                                                         cAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                            >  2:801939/1‑150 (MQ=255)
                                                                                                                         cAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAGTACGCCAGGCCTTCTGCg                            >  2:3638249/1‑150 (MQ=255)
                                                                                                                          aGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                           >  1:1947606/1‑150 (MQ=255)
                                                                                                                            ggggCGGGGGCTATGGCCCCCCCCCG‑CC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCAGGAAAGGAGCGCAACAAGATGATGAATCAACATATAACACCATCTTAAAAAAAGGCCGGACATTACGCCAGGCCTTCTGCGtta                         <  2:2264274/150‑1 (MQ=255)
                                                                                                                            ggggCGGCGGCGCGGGCCCGCCCCCG‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                         <  1:1658882/150‑1 (MQ=39)
                                                                                                                            ggggCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                         >  2:889081/1‑150 (MQ=255)
                                                                                                                             ggggggCGGGGTTTGGCCCCCCCCC‑CCGGCGCTCCTTGCGCCTGGGGCATTTTG‑TTGGTTCGCTCATGAACGGCGGGCACAAAGATGTGGAACCACAATATACAAACAACTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtta                         <  2:356473/150‑1 (MQ=255)
                                                                                                                             gggCGGGGGCTTTGGGCCCCCCCCG‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                        <  1:1707651/150‑1 (MQ=255)
                                                                                                                             gggCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa                        >  1:16064/1‑150 (MQ=255)
                                                                                                                              ggggggggcgATGGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat                       <  2:977770/144‑1 (MQ=255)
                                                                                                                              gggggcggcgTTTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat                       <  2:977771/147‑1 (MQ=255)
                                                                                                                                gggcggctTTTGCCCGCCCCCG‑GC‑GCGCTCCTTTCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta                     <  2:2635928/149‑1 (MQ=255)
                                                                                                                                cggcggcgATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatta                     >  1:1510473/1‑150 (MQ=255)
                                                                                                                                 ggcggcgATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaattaa                    >  1:1813095/1‑150 (MQ=255)
                                                                                                                                 ggcgcctATGGGCCGCCCCAT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaattaa                    <  2:650813/150‑1 (MQ=255)
                                                                                                                                  gcggcgATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                   >  2:1116340/1‑150 (MQ=255)
                                                                                                                                  gcggcgATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                   >  2:1116223/1‑150 (MQ=255)
                                                                                                                                  gcggcgATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTtaattaat                   >  1:1837099/1‑150 (MQ=255)
                                                                                                                                     gctttgggCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAc                <  2:4044980/144‑1 (MQ=255)
                                                                                                                                      ctttgggCCGCCCCCG‑CC‑GCGCTCCTTGCGCCAGTGGAAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg               <  1:1636425/145‑1 (MQ=255)
                                                                                                                                      ctaTTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACg               <  2:1332244/148‑1 (MQ=255)
                                                                                                                                      cgATTGGCCGCCCCCG‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACACCAGGCCTTCTGCGTTAATTAATCACg               <  1:4052408/150‑1 (MQ=255)
                                                                                                                                       gATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAATCAGCGCAAAAAGAAGATGAATCAACAAATAACAACATCTAAAAAAAAGGCCTGACAATACGCCAGCCCTTCACCGATAAGAAATCAccc              >  2:2245234/1‑148 (MQ=255)
                                                                                                                                       gATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc              >  1:2681396/1‑150 (MQ=255)
                                                                                                                                       gATTGGCCGCCCCCT‑GC‑GCGCTCCTTGAGCCAGTTGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAACTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc              <  1:3414632/150‑1 (MQ=255)
                                                                                                                                           ggCCGCCCGGG‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtttt          <  2:2861589/150‑1 (MQ=255)
                                                                                                                                           ggCCGCCCCCG‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCtttt          <  2:2891701/150‑1 (MQ=255)
                                                                                                                                             ccgccCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGCGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAATAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAAGGAATCACGCTTTTcc        >  2:4101339/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc        >  1:1310871/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc        >  1:4086902/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc        >  2:1123847/1‑150 (MQ=255)
                                                                                                                                             ccgccCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTcc        >  2:2696791/1‑150 (MQ=255)
                                                                                                                                               gccCCCG‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcag      <  1:2944163/150‑1 (MQ=255)
                                                                                                                                               ccccccTGGC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCca       <  2:483666/149‑1 (MQ=255)
                                                                                                                                                  cccT‑GC‑GCGCTCCGTGCGCCAGTGGCAATATG‑TTGCTTAGCTCCTGAACGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag   <  2:1030717/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCT‑GC‑GCGCTCCTTGCGCCAGTGGCAATATG‑TTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGG  >  NC_000913/915077‑915373

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