Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,743,157 G→A 7.5% P29L (CCC→CTC)  ribC ← riboflavin synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,743,1570GA 7.5% 285.2 / 7.8 93P29L (CCC→CTC) ribCriboflavin synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (42/43);  new base A (5/2);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

CCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCT  >  NC_000913/1743008‑1743306
                                                                                                                                                     |                                                                                                                                                     
ccACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCgg                                                                                                                                                       >  1:528165/1‑150 (MQ=255)
 cACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCggg                                                                                                                                                      <  1:502952/150‑1 (MQ=255)
  aCTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGt                                                                                                                                                     <  1:1502460/150‑1 (MQ=255)
    ttttAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTaa                                                                                                                                                   >  1:470022/1‑150 (MQ=255)
      ttAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACt                                                                                                                                                 <  2:1536563/150‑1 (MQ=255)
            cGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTccacca                                                                                                                                           <  1:1103808/150‑1 (MQ=255)
            cGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTccacaa                                                                                                                                           >  2:1181108/1‑148 (MQ=255)
              ccAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATATCCGGCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTGCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACcaca                                                                                                                                         <  1:68705/150‑1 (MQ=255)
                aaGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATg                                                                                                                                       <  2:1277680/150‑1 (MQ=255)
                aaGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATg                                                                                                                                       <  1:164029/150‑1 (MQ=255)
                aaGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATg                                                                                                                                       <  2:500165/150‑1 (MQ=255)
                 aGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACTGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGc                                                                                                                                      <  2:1178853/150‑1 (MQ=255)
                  gATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCg                                                                                                                                     <  1:1526257/150‑1 (MQ=255)
                  gATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCg                                                                                                                                     >  1:511744/1‑150 (MQ=255)
                    ttGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTa                                                                                                                                   >  2:770812/1‑150 (MQ=255)
                    ctGGTACTGCGTAACCCCCCTTTCATCACCTCAAACCCCACAAGCCAGCCGCAAATATCCGCCACGGTCACGCAGCAACCGTTATGCGCCACGGAAGCACCCGCCTCCAGGCCGTCCCCCAGGCGGTCGGGTAACTCCACCACACGCGTa                                                                                                                                   <  2:1307405/149‑1 (MQ=255)
                     tGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTAc                                                                                                                                  <  2:369434/150‑1 (MQ=255)
                     tGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTAc                                                                                                                                  <  2:1587530/150‑1 (MQ=255)
                       gTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGa                                                                                                                                >  1:194569/1‑150 (MQ=255)
                          aTGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGaaaa                                                                                                                             >  1:398699/1‑150 (MQ=255)
                             cGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAAttt                                                                                                                          >  1:1582989/1‑150 (MQ=255)
                              gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTg                                                                                                                         <  1:1181108/150‑1 (MQ=255)
                              gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTg                                                                                                                         <  2:716837/150‑1 (MQ=255)
                                         ttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCa                                                                                                              <  2:705761/150‑1 (MQ=255)
                                         ttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCa                                                                                                              >  2:1545441/1‑150 (MQ=255)
                                         ttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCa                                                                                                              >  1:1656438/1‑150 (MQ=255)
                                             tcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCg                                                                                                          <  2:1221922/150‑1 (MQ=255)
                                               aGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGac                                                                                                        >  1:1697653/1‑150 (MQ=255)
                                                          gACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGc                                                                                             >  2:467427/1‑150 (MQ=255)
                                                          gACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGc                                                                                             >  1:1532972/1‑150 (MQ=255)
                                                           aCATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCg                                                                                            <  1:437052/150‑1 (MQ=255)
                                                              tGGTTGCCGTTAATTTCCGTCATGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTg                                                                                         >  2:595623/1‑150 (MQ=255)
                                                              tGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTg                                                                                         >  1:1200821/1‑150 (MQ=255)
                                                               ggTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGc                                                                                        >  2:1512291/1‑150 (MQ=255)
                                                               ggTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGc                                                                                        >  2:635012/1‑150 (MQ=255)
                                                                 ttGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGccc                                                                                      <  1:1545441/150‑1 (MQ=255)
                                                                 ttGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGccc                                                                                      <  2:398699/150‑1 (MQ=255)
                                                                   gCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTg                                                                                    >  1:344450/1‑150 (MQ=255)
                                                                    ccGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTGTGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGt                                                                                   >  2:21163/1‑150 (MQ=255)
                                                                    ccGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGt                                                                                   >  2:1283378/1‑150 (MQ=255)
                                                                     cGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACAGTTTTGCGGTGCCCTGTAc                                                                                 <  2:7385/150‑1 (MQ=255)
                                                                       ttatattCCGTCAAGGTCAGGGAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCGAGGCCGTCGAGCATGTGGTCGGGTAACTCGACCTCATGCGTTCGAAAATGTGGGTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTAct                                                                                <  1:1112570/145‑2 (MQ=255)
                                                                          aTTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATc                                                                             >  1:1331786/1‑150 (MQ=255)
                                                                           tttACGTCACGGTCAGTAAGCAACCTTTATGAGCCACGGAAGCAACGGTATCAAGGCCGTCAAGCATGTGGTCGGGTAACTCCACCACATGCTTAAGAAAATTTGGTTTCTCGTCAATCGACAACAGTTTTGCGGTGCCCTGTACAATcc                                                                            >  2:995924/1‑150 (MQ=255)
                                                                             tCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTTCCCTGTACAATcccc                                                                          <  2:291594/150‑1 (MQ=255)
                                                                             tCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcccc                                                                          <  1:1522630/150‑1 (MQ=255)
                                                                              ccGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCg                                                                         <  2:536617/150‑1 (MQ=255)
                                                                                gTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcccc                                                                          <  2:319886/147‑1 (MQ=255)
                                                                                gTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcccc                                                                          >  1:319886/1‑147 (MQ=255)
                                                                                   aCGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAAc                                                                    <  1:1360002/150‑1 (MQ=255)
                                                                                   aCGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAAc                                                                    <  2:194569/150‑1 (MQ=255)
                                                                                      gTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATa                                                                 <  2:1582989/150‑1 (MQ=255)
                                                                                          ggCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTa                                                             >  1:571231/1‑150 (MQ=255)
                                                                                          ggCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTa                                                             >  1:316672/1‑150 (MQ=255)
                                                                                            cagcaACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACa                                                           <  2:511744/150‑1 (MQ=255)
                                                                                              gcaACCGTTGTGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAAc                                                         <  1:21163/150‑1 (MQ=255)
                                                                                               caACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACt                                                        <  1:1283378/150‑1 (MQ=255)
                                                                                               caACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACt                                                        >  2:1007187/1‑150 (MQ=255)
                                                                                                 aCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTcc                                                      >  1:1082642/1‑150 (MQ=255)
                                                                                                  ccGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCt                                                     <  1:870594/150‑1 (MQ=255)
                                                                                                     ttATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGaa                                                  <  2:1697653/150‑1 (MQ=255)
                                                                                                      tATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGaaa                                                 <  1:1512291/150‑1 (MQ=255)
                                                                                                        tGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATc                                               >  1:1466835/1‑150 (MQ=255)
                                                                                                         gcgcCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCTCCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCa                                              >  2:840443/1‑150 (MQ=255)
                                                                                                          cgcCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAg                                             >  2:309212/1‑150 (MQ=255)
                                                                                                           gcCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTAGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGt                                            <  2:1127764/150‑1 (MQ=255)
                                                                                                             cACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGTGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTa                                          <  2:1420433/150‑1 (MQ=255)
                                                                                                             cACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTa                                          >  2:1398443/1‑150 (MQ=255)
                                                                                                                  aaGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGACAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACa                                     <  2:858836/150‑1 (MQ=255)
                                                                                                                   aGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACAt                                    >  2:1004372/1‑150 (MQ=255)
                                                                                                                      aCCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCt                                 >  1:624000/1‑150 (MQ=255)
                                                                                                                          gTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTc                             <  2:1331786/150‑1 (MQ=255)
                                                                                                                             tCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGc                          >  2:160472/1‑150 (MQ=255)
                                                                                                                             tCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGc                          >  1:391414/1‑150 (MQ=255)
                                                                                                                              ccAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACAAACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGcc                         >  1:1295371/1‑149 (MQ=255)
                                                                                                                                aGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGcaca                       >  1:136479/1‑150 (MQ=255)
                                                                                                                                  gCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAAt                     <  1:1309148/150‑1 (MQ=255)
                                                                                                                                  gCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAAt                     <  2:1200821/150‑1 (MQ=255)
                                                                                                                                  gCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTAGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAAt                     >  2:237528/1‑150 (MQ=255)
                                                                                                                                    cGTCCAGCATGTGGTCGGGTAACTCCACCGCATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAg                   >  2:1021996/1‑150 (MQ=255)
                                                                                                                                    cGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAg                   >  1:925836/1‑150 (MQ=255)
                                                                                                                                      tCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCa                 >  1:260443/1‑150 (MQ=255)
                                                                                                                                      gtacgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCa                 >  1:1292501/5‑150 (MQ=255)
                                                                                                                                      gtacgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCa                 >  1:939261/5‑150 (MQ=255)
                                                                                                                                      gtacgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCa                 >  2:800470/5‑150 (MQ=255)
                                                                                                                                       tacgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAg                <  1:800470/147‑1 (MQ=255)
                                                                                                                                       ccAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGGTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAg                >  1:728622/1‑150 (MQ=255)
                                                                                                                                         cgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGt              >  2:426981/2‑150 (MQ=255)
                                                                                                                                         cgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGt              >  1:1183419/2‑150 (MQ=255)
                                                                                                                                          gCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTg             <  2:1656438/150‑1 (MQ=255)
                                                                                                                                           cATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTgg            <  2:316672/150‑1 (MQ=255)
                                                                                                                                           cATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTgg            <  1:595623/150‑1 (MQ=255)
                                                                                                                                                ggTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTATGACATTCTGTGCAGCACAATAGCAGGTGGaaacc       >  2:1441483/1‑148 (MQ=255)
                                                                                                                                                  tCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGc     <  2:925836/150‑1 (MQ=255)
                                                                                                                                                   cGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGcc    <  2:470022/150‑1 (MQ=255)
                                                                                                                                                     ggTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCt  >  2:1500621/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCT  >  NC_000913/1743008‑1743306

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