Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,899,866 T→A 5.1% K238* (AAA→TAA) ‡ yieL ← putative hydrolase YieL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,899,8660TA 5.1% 510.9 / 11.5 156K238* (AAA→TAA) ‡yieLputative hydrolase YieL
Reads supporting (aligned to +/- strand):  ref base T (68/80);  new base A (5/3);  total (73/83)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 160 total aligned reads displayed.

CCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGAAC  >  NC_000913/3899719‑3900015
                                                                                                                                                    |                                                                                                                                                     
ccAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATAAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTAc                                                                                                                                                     >  2:1687263/1‑150 (MQ=255)
 cAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACg                                                                                                                                                    <  1:1554585/150‑1 (MQ=255)
 cAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACg                                                                                                                                                    >  2:1659692/1‑150 (MQ=255)
        tGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTcc                                                                                                                                             <  2:2220207/150‑1 (MQ=255)
        tGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTcc                                                                                                                                             <  1:1786564/150‑1 (MQ=255)
         ggTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCt                                                                                                                                            <  1:2333397/150‑1 (MQ=255)
         ggTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCt                                                                                                                                            >  2:1166647/1‑150 (MQ=255)
         ggTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCt                                                                                                                                            >  1:1186764/1‑150 (MQ=255)
          gTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTg                                                                                                                                           >  1:1959884/1‑150 (MQ=255)
          gTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACGTAAAGTCttttatc                                                                                                                                           <  1:917940/150‑5 (MQ=255)
             cccGCCTTGTGATAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAgg                                                                                                                                        <  2:1884669/150‑1 (MQ=255)
              ccGCCTTGTGAATGACCTGCCAGCGCGCGTCCATCGGCATCTTTACGGACATTAAAACTCTTGCTAATCAGCTTGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCTGATAAAAGACTTTACTTCTTTCCTGAGGc                                                                                                                                       <  2:1627649/150‑1 (MQ=255)
                 ccTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACg                                                                                                                                    >  1:1278218/1‑150 (MQ=255)
                 ccTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACg                                                                                                                                    <  1:605859/150‑1 (MQ=255)
                    tgtgAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGaaa                                                                                                                                 <  2:1769536/150‑1 (MQ=255)
                     gtgAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAAt                                                                                                                                <  2:2070040/150‑1 (MQ=255)
                       gAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCt                                                                                                                              <  1:1245153/150‑1 (MQ=255)
                        aaaGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAACACTTTACGTCTTTCCTGAGGCACGAAATCtt                                                                                                                             <  1:1160023/150‑1 (MQ=255)
                          aGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCg                                                                                                                           <  1:719387/150‑1 (MQ=255)
                           gCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTCCGTCTTTCCTGAGGCACGAAATCTTCgg                                                                                                                          >  2:1565320/1‑150 (MQ=255)
                           gCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCgg                                                                                                                          <  1:1759221/150‑1 (MQ=255)
                            cccTGCCAGCGCGCGGCCATCGTCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCggg                                                                                                                         <  1:216777/150‑1 (MQ=255)
                            cccTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCTTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCggg                                                                                                                         <  1:1097574/150‑1 (MQ=255)
                            cccTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCggg                                                                                                                         <  1:1018150/150‑1 (MQ=255)
                            cccTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGATTGCTAATCAGCGGGATAATATGGTTCATCAGTTCGCGATCTGACGCTTTAGCATTCAGCGGATAAAAGACTTTAGGTCTTTCCGGAGGCACGAAATCTTCggg                                                                                                                         >  2:440999/1‑150 (MQ=255)
                             ccTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGa                                                                                                                        >  1:991878/1‑150 (MQ=255)
                              cTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGaa                                                                                                                       >  1:1479741/1‑150 (MQ=255)
                                      gcgcgGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGcac                                                                                                               >  1:1298409/1‑148 (MQ=255)
                                           gCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGc                                                                                                          <  2:529275/150‑1 (MQ=255)
                                           gCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGc                                                                                                          >  2:989446/1‑150 (MQ=255)
                                           gCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGc                                                                                                          >  1:1607669/1‑150 (MQ=255)
                                            ccATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCa                                                                                                         <  1:1560814/150‑1 (MQ=255)
                                             cATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAGTAATGCCCTTCGCAt                                                                                                        >  1:1393837/1‑150 (MQ=255)
                                              aTCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATc                                                                                                       <  1:661898/150‑1 (MQ=255)
                                                     tCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCt                                                                                                <  1:1326317/150‑1 (MQ=255)
                                                      cTTTACGGACATTAAAACGCTTGCTAATTAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTg                                                                                               >  2:884619/1‑150 (MQ=255)
                                                         tACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTAt                                                                                            <  2:139688/150‑1 (MQ=255)
                                                              aCATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCgg                                                                                        >  2:2162778/1‑150 (MQ=255)
                                                              aCATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCggg                                                                                       <  2:1496233/150‑1 (MQ=255)
                                                                 ttAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATc                                                                                    >  2:2125653/1‑150 (MQ=255)
                                                                 ttAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATc                                                                                    >  2:1521022/1‑150 (MQ=255)
                                                                  tAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATca                                                                                   <  1:73219/150‑1 (MQ=255)
                                                                  tAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATca                                                                                   >  2:628961/1‑150 (MQ=255)
                                                                  tAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATca                                                                                   <  1:1655279/150‑1 (MQ=255)
                                                                   aaaaCGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATcac                                                                                  >  2:114564/1‑150 (MQ=255)
                                                                   aaaaCGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATcac                                                                                  <  1:371008/150‑1 (MQ=255)
                                                                    aaaCGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATcacc                                                                                 >  1:1874706/1‑150 (MQ=255)
                                                                    aaaCGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATcacc                                                                                 <  2:110419/150‑1 (MQ=255)
                                                                      aCGCTTGCTAATCAGCGGGTTAATATCGTTCATCAGTTCTCGATCTGCCGCTTTATCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATcaccac                                                                               <  2:1901160/150‑1 (MQ=255)
                                                                       cGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCaccacc                                                                              <  1:2434590/150‑1 (MQ=255)
                                                                         cTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAg                                                                            <  2:86460/150‑1 (MQ=255)
                                                                           tGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCa                                                                          >  1:1451607/1‑150 (MQ=255)
                                                                                   aGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTa                                                                  <  2:1607669/150‑1 (MQ=255)
                                                                                    gCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTaa                                                                 >  1:2118627/1‑150 (MQ=255)
                                                                                      gggATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAAtt                                                               <  2:1479741/150‑1 (MQ=255)
                                                                                        gATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAAtttt                                                             >  2:511967/1‑150 (MQ=255)
                                                                                          taataTCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTcc                                                           <  1:2457789/150‑1 (MQ=255)
                                                                                           aataTCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTccc                                                          <  2:1031135/150‑1 (MQ=255)
                                                                                           aataTCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTccc                                                          <  2:1268483/150‑1 (MQ=255)
                                                                                            atatCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCt                                                         <  1:1659692/150‑1 (MQ=255)
                                                                                            atatCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCt                                                         <  2:551045/150‑1 (MQ=255)
                                                                                            atatCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGGCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCt                                                         <  1:1054838/150‑1 (MQ=255)
                                                                                              atCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTc                                                       <  2:733632/150‑1 (MQ=255)
                                                                                                 gTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCagc                                                    <  1:2112701/150‑1 (MQ=255)
                                                                                                 gTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCagc                                                    >  1:338904/1‑150 (MQ=255)
                                                                                                      tcaGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCagcaagca                                               <  2:2052518/150‑1 (MQ=255)
                                                                                                       caGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAg                                              >  2:1461001/1‑150 (MQ=255)
                                                                                                        aGTTCGCGTTCTGGCGCTTTAGCATTCTGCGGATAAAAGTCTTTACGTCTTTCCTGAGGCACGAATTCTTCGGGATTAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAgg                                             <  1:961771/150‑1 (MQ=255)
                                                                                                             gcgATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATc                                        <  1:1986724/150‑1 (MQ=255)
                                                                                                             gcgATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATc                                        <  2:2118627/150‑1 (MQ=255)
                                                                                                              cgATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATcc                                       >  2:1321380/1‑150 (MQ=255)
                                                                                                               gATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCa                                      >  2:2287108/1‑150 (MQ=255)
                                                                                                                aTCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCAt                                     <  2:1959884/150‑1 (MQ=255)
                                                                                                                  cTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGa                                   >  2:727012/1‑150 (MQ=255)
                                                                                                                   tGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGAt                                  >  1:310516/1‑150 (MQ=255)
                                                                                                                    gCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGAtt                                 >  1:1251824/1‑150 (MQ=255)
                                                                                                                    gCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGAtt                                 >  2:2168602/1‑150 (MQ=255)
                                                                                                                            aaGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCTGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCAATGATTTGCGGGa                          >  2:2372118/2‑150 (MQ=255)
                                                                                                                              gCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATa                        <  1:511967/150‑1 (MQ=255)
                                                                                                                               cATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATAc                       <  2:1286973/150‑1 (MQ=255)
                                                                                                                                aTTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACg                      >  1:838926/1‑150 (MQ=255)
                                                                                                                                  tCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGcc                    <  1:1379081/150‑1 (MQ=255)
                                                                                                                                  tCAGCGGACAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGcc                    <  2:249800/150‑1 (MQ=255)
                                                                                                                                      ctcaggaaaGACGTAAAGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTg                >  2:1745530/7‑150 (MQ=255)
                                                                                                                                      cGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTg                >  1:2159748/1‑150 (MQ=255)
                                                                                                                                       tcaggaaaGACGTAAAGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTgg               >  1:1325589/6‑149 (MQ=255)
                                                                                                                                         aggaaaGACGTAAAGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATc             >  1:2105008/4‑150 (MQ=255)
                                                                                                                                         aTAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATc             <  1:1865659/150‑1 (MQ=255)
                                                                                                                                          tAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCg            >  2:1572635/1‑150 (MQ=255)
                                                                                                                                          tAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCg            >  1:602381/1‑150 (MQ=255)
                                                                                                                                           aaaaGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGa           >  1:439417/1‑150 (MQ=255)
                                                                                                                                            aaaGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGAt          >  1:1312592/1‑150 (MQ=255)
                                                                                                                                              aGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAg        >  2:1607328/1‑150 (MQ=255)
                                                                                                                                               gACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGc       >  2:271383/1‑150 (MQ=255)
                                                                                                                                                 cTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTTTATCCGGGATCAAAACCATCATAGGTATAATTTTCCCTTCAGCAAGCAAGTTATCCATGATTAGCGGGATAAGCCCCTGATCGATAGCgg     >  1:1591060/1‑150 (MQ=255)
                                                                                                                                                 cTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCgg     <  2:1369539/150‑1 (MQ=255)
                                                                                                                                                   ttACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGaa   >  2:2174967/1‑150 (MQ=255)
                                                                                                                                                   ttACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGaa   <  2:1186764/150‑1 (MQ=255)
                                                                                                                                                    tACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGAAc  >  2:563676/1‑150 (MQ=255)
                                                                                                                                                    tACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGAAc  >  1:1420638/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAGCGCCTGGTACCCGCCTTGTGAAAGCCCTGCCAGCGCGCGGCCATCGGCATCTTTACGGACATTAAAACGCTTGCTAATCAGCGGGATAATATCGTTCATCAGTTCGCGATCTGCCGCTTTAGCATTCAGCGGATAAAAGACTTTACGTCTTTCCTGAGGCACGAAATCTTCGGGAATAATGCCCTTCGCATCGGTTTCTGTATCCGGGATCACCACCAGCATCGGTTTAATTTTCCCTTCAGCAAGCAGGTTATCCATGATTTGCGGGATACGCCCCTGATCGATAGCGGAAC  >  NC_000913/3899719‑3900015

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: