Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,132,077 G→A A111V (GCG→GTG)  yghQ ← putative transport protein YghQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,132,0770GA93.3% 178.2 / NA 60A111V (GCG→GTG) yghQputative transport protein YghQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/32);  minor base C (0/4);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCGCCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCA  >  NC_000913/3131837‑3132317
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
cccGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGcgcc                                                                                                                                                                                                                                         <  1:21365/251‑1 (MQ=255)
           aaattCAAAATCCCAGGCTAAGACTCTCCCCGCTGCTGGCGGAAAAGGTTTCGTCGCCTGTAGTACAGCAAGTAAAGCTAAGCGATCTACCGCACGCAGAATGCCGGGCGGCGCGGAGGAAGCCATGGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGCCCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCCCCGCGCCGCTGACGatat                                                                                                                                                                                                                              <  2:215743/247‑1 (MQ=255)
            aaGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATc                                                                                                                                                                                                                             <  1:300305/251‑1 (MQ=255)
                 aaaaTACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCg                                                                                                                                                                                                                        <  1:165888/251‑1 (MQ=255)
                      aCCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCgaga                                                                                                                                                                                                                     >  1:3699/1‑249 (MQ=255)
                         aGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGa                                                                                                                                                                                                                 >  1:558246/1‑250 (MQ=255)
                         aGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGa                                                                                                                                                                                                                 >  1:254678/1‑250 (MQ=255)
                             tacaacgCTCCCCGCTTCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAAGCGAAGCGATCTCCCGCACGCAGAAGGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCTCTGCCAGCCAAAAACGTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGTCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAAt                                                                                                                                                                                                             <  2:57201/246‑1 (MQ=255)
                                gacgCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGACGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGCAAGGAATGGCAGTAAGGCAATGCCACCGCAAAACACCACCGCGCCGCTGACGAACTACAGCGAGAAGGcaaagg                                                                                                                                                                                                            >  2:6562/1‑246 (MQ=255)
                                  cgCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAgacg                                                                                                                                                                                                       <  1:291270/251‑1 (MQ=255)
                                     cccccTCTGCGCGCAGAAAAGTTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTTCCGCACGCAGAATGCCGGTCGTCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGTCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACa                                                                                                                                                                                                    <  2:257483/250‑1 (MQ=255)
                                       cccGCTGAGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAAGCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCGGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACAtt                                                                                                                                                                                                  <  2:476034/251‑1 (MQ=255)
                                             gcgcgcAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCgg                                                                                                                                                                                              <  1:3209/249‑1 (MQ=255)
                                              cgcgcAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAAt                                                                                                                                                                                           >  1:472991/1‑251 (MQ=255)
                                               gcgcAGAAAAGGTTTCGCCGCCGTCTGTACAGCAATTAAAGCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAAtt                                                                                                                                                                                          <  2:406891/251‑1 (MQ=255)
                                                   aGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATtgctg                                                                                                                                                                                      >  2:210870/1‑251 (MQ=255)
                                                    gAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATtgctgc                                                                                                                                                                                     <  1:47046/251‑1 (MQ=255)
                                                             ccgtcgCCTGCTGTACAGCAATGAAATCGAAGCGATCGACCGCACGCAGAATGCCGGTCGTCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTGTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGTCAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGattatt                                                                                                                                                                             <  1:352098/249‑1 (MQ=255)
                                                                      tctgTACAGCTAATAAATCGACGCGATCTACCGCACGCAGAATGCCGGTCCGCGTGAGTGAAGCAATTGAAGGAACGAGCGTGTAAGAGAGCGCTGCCAGCCAACATCTTTAGTAATCCCGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCCCCACCGCGCCGATGACTATATCCATCTATTAGGAAAATGAGACGACATTGCTGATTTGCTTCGGATTATTTTTGGTTa                                                                                                                                                                     <  2:424407/248‑1 (MQ=255)
                                                                      gctgTACAGCAATTAAATCGAAGCGATCTACCTCCCGCAGAATGCCGGGCGGCGGTGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAAGGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTa                                                                                                                                                                     <  2:424394/249‑1 (MQ=255)
                                                                      gcggTACAGCAATTAAAGCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGGCCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCAACGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAAt                                                                                                                                                                   <  2:3699/248‑1 (MQ=255)
                                                                       ctgTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATg                                                                                                                                                                  >  1:300482/1‑251 (MQ=255)
                                                                         gTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATg                                                                                                                                                                  <  1:536140/249‑1 (MQ=255)
                                                                              gCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGgtgt                                                                                                                                                           >  1:413376/1‑251 (MQ=255)
                                                                               cAATTAAATCGAAGCGGTCTACCGCACTCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTTCAATATATCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCCCCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGtt                                                                                                                                                          <  2:379451/251‑1 (MQ=255)
                                                                                      aTCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACt                                                                                                                                                   <  1:480739/251‑1 (MQ=255)
                                                                                          aaGCGATCTACCGCACGCAGAATGCCGGTCGGCGGGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGTCCTAATGAAGGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCATCGAGAAGGAAAATGAGACGACATTGCGGAAGTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAac                                                                                                                                               <  1:39355/251‑1 (MQ=255)
                                                                                               aTCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATtgctgc                                                                                                                                                                                     <  1:322580/208‑1 (MQ=255)
                                                                                               aTCTACCGCACGCAGAACGACGGTCGGCGAGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAACAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATtgctgc                                                                                                                                                                                     >  2:322580/1‑208 (MQ=255)
                                                                                                             gAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGt                                                                                                                                                                       >  1:186544/1‑208 (MQ=255)
                                                                                                             gAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGt                                                                                                                                                                       <  2:186544/208‑1 (MQ=255)
                                                                                                               aTGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGaaa                                                                                                                          <  1:103834/251‑1 (MQ=255)
                                                                                                                       cGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTtaataa                                                                                                                  <  1:164041/251‑1 (MQ=255)
                                                                                                                          cGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAAt                                                                                                               <  1:209689/251‑1 (MQ=255)
                                                                                                                                ggaAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAAATTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTg                                                                                                          >  1:487786/1‑250 (MQ=255)
                                                                                                                                 ggacgccatagaaggaaggggctggCAAAAGAGCGCTGCCAGCCAAAAACTTGGGTCATCCAGACCTAATGAAAGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAAGCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCGGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGa                                                                                                         <  2:274777/227‑1 (MQ=255)
                                                                                                                                       cATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACt                                                                                                   >  1:434841/1‑251 (MQ=255)
                                                                                                                                             aGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGt                                                                                             <  1:412030/251‑1 (MQ=255)
                                                                                                                                                   gggctttcAATAGAGCGCGACCAGCCAAAAACTTTGGTCATACAGAACAAAGGAATGGGAAAGGTATAGCATTAAGGCAATGACACCGACAATCCCCACCGCTCCGCGGACGATATCCAGCGAGAAGGAAAATGAGACTTCATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCTCCAGGTGCCATGGCTGAAACTTAATAAAAGCGCTGATCGACTTGGCGTACGAt                                                                                        <  2:344270/243‑1 (MQ=255)
                                                                                                                                                          cAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAAt                                                                                <  1:442495/251‑1 (MQ=255)
                                                                                                                                                          cAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGATGGTGTTCCGTCCTGAAACACCAGTTGCCATGTATGAAACTTAATAACATCGCTGATCGACTTGGCGTACGATTGCACaa                                                                                 >  2:125952/1‑250 (MQ=255)
                                                                                                                                                                   gcTGCGAGGCAAAAAATTTGGTCAGCCAGACCTAATGAAGGGGAAAGGAATGGCAGTAAGGCAAGGCCACCGACAATCCCCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAAGTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTAc                                                                       <  2:254678/251‑1 (MQ=255)
                                                                                                                                                                      gccagtcaAAAACTTTGGTCACCCGGACCTAATGAATGTGAAAGGACGGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGGTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCTCCAGTACGcc                                                                    <  2:267251/251‑1 (MQ=255)
                                                                                                                                                                               aaaCTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCACTGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATAGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGCACCACCAGTTGCCATGTCTGAAACTTAATAAACTCGCTGATCGCCTTCGCGAACGACTGCACCATCACCAGTACGCCAAACATgcg                                                           >  2:172017/1‑249 (MQ=255)
                                                                                                                                                                                    ttGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATggcgggc                                                       <  1:517655/250‑1 (MQ=255)
                                                                                                                                                                                         cATCCAGAACAAATGAATGGGAAAGGAATGGCAGTAATGCAATGCCACCGACAAGCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAAGAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCAt                                                 <  2:518294/250‑1 (MQ=255)
                                                                                                                                                                                          aTCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAACCTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATggcgggcttcctcc                                                >  1:124226/1‑244 (MQ=255)
                                                                                                                                                                                                      tgaatgGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGGCATCCCTTTACCGGcac                                   >  1:59676/1‑251 (MQ=255)
                                                                                                                                                                                                            ggAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGCCTTGCAGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGTcacccg                                >  1:287503/1‑245 (MQ=255)
                                                                                                                                                                                                             gAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACgcg                              >  1:530354/1‑247 (MQ=255)
                                                                                                                                                                                                              aaaGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAgc                           <  1:6562/251‑1 (MQ=255)
                                                                                                                                                                                                                aGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAgcgcgc                         >  1:381166/1‑251 (MQ=255)
                                                                                                                                                                                                                aGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCAACAGTACGCCAAACATGGCGGGCGTCCACCCCTTTAcggccccccaggcagc                          >  2:242980/1‑236 (MQ=255)
                                                                                                                                                                                                                       gCAGTACGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCCCCAGTTg                                                                                                                                      >  2:375283/1‑136 (MQ=255)
                                                                                                                                                                                                                       gCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACAGTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTg                                                                                                                                      <  1:375283/136‑1 (MQ=255)
                                                                                                                                                                                                                             aggaaaTGCCCCCGACAAGCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATGGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAg                                                                          <  2:472027/186‑1 (MQ=255)
                                                                                                                                                                                                                             aGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAg                                                                          >  1:472027/1‑190 (MQ=255)
                                                                                                                                                                                                                              ggaaaTGCCACCGACAATCACCACCGAGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGGCGACATTGCGGAATTGCTTCGGATTATTGGTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCGTTACCGGCACACGAGAGCGCCAACAGACCCAGCa           <  2:343490/248‑1 (MQ=255)
                                                                                                                                                                                                                                 aaTGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACTCCATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGTCGTACGCTTGCACAATCACCAGTACGCCAAACATGGCGGGCGTTATCCCCTTACCGGCACCCgagggccccacc                     >  1:298895/1‑229 (MQ=255)
                                                                                                                                                                                                                                       aCCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCCAACGAGAGCGACAACAGACCCAGCAAGGCGCTcc  >  1:160603/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
CCCGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACTTTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCGCCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCA  >  NC_000913/3131837‑3132317

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