Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,389,337 A→G C12R (TGC→CGC)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,389,3370AG97.8% 326.7 / ‑6.4 91C12R (TGC→CGC) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (43/46);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAGGTTCAGCGTCTG  >  NC_000913/4389190‑4389483
                                                                                                                                                   |                                                                                                                                                  
cTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGcc                                                                                                                                                  >  1:1273900/1‑150 (MQ=255)
  tttCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCa                                                                                                                                                 <  1:79779/149‑1 (MQ=255)
  tttCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGCTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCa                                                                                                                                                 <  2:63590/149‑1 (MQ=255)
    tCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGc                                                                                                                                              <  2:1122024/150‑1 (MQ=255)
      gTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCAt                                                                                                                                           <  2:311080/151‑1 (MQ=255)
        tCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATc                                                                                                                                          >  2:726308/1‑150 (MQ=255)
          cTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGCCTTCCGGCTGAGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGCGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACTCCCCCCAACTGAGGCGCCCGGCCATc                                                                                                                                          >  1:1261354/1‑148 (MQ=255)
          cTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATc                                                                                                                                          >  2:1032209/1‑148 (MQ=255)
            cAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGaa                                                                                                                                      >  1:712811/1‑150 (MQ=255)
            cAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGa                                                                                                                                       <  2:1234668/149‑1 (MQ=255)
            cAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGa                                                                                                                                       <  2:624425/149‑1 (MQ=255)
             aaGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGaaaa                                                                                                                                    <  1:180413/151‑1 (MQ=255)
             aaGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGaaa                                                                                                                                     >  2:218068/1‑150 (MQ=255)
                  cATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGAt                                                                                                                               >  1:769559/1‑151 (MQ=255)
                   aTTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGAt                                                                                                                               >  1:825151/1‑150 (MQ=255)
                    tttAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATa                                                                                                                              >  1:1244124/1‑150 (MQ=255)
                    tttAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATa                                                                                                                              >  2:104077/1‑150 (MQ=255)
                       aaGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATc                                                                                                                           <  2:300458/150‑1 (MQ=255)
                       aaGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCa                                                                                                                          <  2:900603/151‑1 (MQ=255)
                        aGGCAGCCTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAAt                                                                                                                            >  1:1177158/1‑148 (MQ=255)
                        aGGCAGACTGGCGAGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGCGGCACAAGGCCATCAGCAAAGCGATaa                                                                                                                             >  1:933179/1‑147 (MQ=255)
                           cAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGgcg                                                                                                                      >  1:943232/1‑151 (MQ=255)
                              aCTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCAc                                                                                                                   >  1:1115628/1‑151 (MQ=255)
                                tGGAGCCCCTCTACGACTTCCGGCTGGGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCa                                                                                                                    <  1:615601/148‑1 (MQ=255)
                                 ggAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACg                                                                                                                  <  2:503044/149‑1 (MQ=255)
                                 ggAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACg                                                                                                                  <  2:908565/149‑1 (MQ=255)
                                 ggAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCa                                                                                                               <  1:959730/151‑1 (MQ=255)
                                       ccTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGtt                                                                                                          <  1:290071/151‑1 (MQ=255)
                                          ctACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCCATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGttt                                                                                                         <  1:1248019/149‑1 (MQ=255)
                                           tACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCt                                                                                                      >  2:358541/1‑151 (MQ=255)
                                            aCGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTc                                                                                                        <  1:844755/148‑1 (MQ=255)
                                            aCGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCtt                                                                                                     >  2:157574/1‑151 (MQ=255)
                                             cGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTc                                                                                                        <  2:1231509/147‑1 (MQ=255)
                                               aCTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTg                                                                                                   >  1:312140/1‑150 (MQ=255)
                                                  tCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATg                                                                                                <  2:537890/150‑1 (MQ=255)
                                                  tCCGGCTGGGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATgg                                                                                               >  2:77938/1‑151 (MQ=255)
                                                   ccGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATgg                                                                                               >  2:268831/1‑150 (MQ=255)
                                                     ggCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGAtt                                                                                            >  1:1017712/1‑151 (MQ=255)
                                                         gtgCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAg                                                                                          <  2:657446/149‑1 (MQ=255)
                                                           gCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGAcc                                                                                       <  1:121358/150‑1 (MQ=255)
                                                           gCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGAcc                                                                                       <  2:951462/150‑1 (MQ=255)
                                                               gTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTc                                                                                  <  1:1189848/151‑1 (MQ=255)
                                                                 ttCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCtt                                                                                >  2:470115/1‑151 (MQ=255)
                                                                 ttCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCtt                                                                                >  1:1001383/1‑151 (MQ=255)
                                                                   cgccgcTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCtttt                                                                              >  2:1069026/1‑151 (MQ=255)
                                                                   cgccgcTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCtttt                                                                              >  2:647125/1‑151 (MQ=255)
                                                                     ccgcTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCtttttt                                                                            >  2:1199176/1‑151 (MQ=255)
                                                                      cgcTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTg                                                                           >  2:244316/1‑151 (MQ=255)
                                                                      cgcTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTg                                                                           <  1:162524/151‑1 (MQ=255)
                                                                          tttGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTgtcgt                                                                       >  2:609915/1‑151 (MQ=255)
                                                                             gCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCaa                                                                    <  2:467735/151‑1 (MQ=255)
                                                                             gCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCa                                                                     <  1:682907/150‑1 (MQ=255)
                                                                                 gCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCa                                                                 >  2:242176/1‑150 (MQ=255)
                                                                                 gCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCa                                                                 <  1:92193/150‑1 (MQ=255)
                                                                                  cTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAAt                                                               >  2:258066/1‑151 (MQ=255)
                                                                                   tCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATg                                                              >  1:306596/1‑151 (MQ=255)
                                                                                    ccAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAATTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAAt                                                               <  1:741745/149‑1 (MQ=255)
                                                                                    ccAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATgg                                                             >  2:114176/1‑151 (MQ=255)
                                                                                        ttCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTg                                                         >  2:204656/1‑151 (MQ=255)
                                                                                         tCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTgg                                                        >  1:490113/1‑151 (MQ=255)
                                                                                          ccTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTgg                                                        <  2:288217/150‑1 (MQ=255)
                                                                                          ccTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGc                                                       >  2:256258/1‑151 (MQ=255)
                                                                                             gAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGcgtc                                                    <  2:960540/151‑1 (MQ=255)
                                                                                                 gATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCgtgt                                                <  2:968363/151‑1 (MQ=255)
                                                                                                 gATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCgtgt                                                >  1:999903/1‑151 (MQ=255)
                                                                                                   tttgtttgCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTc                                              >  2:617581/1‑151 (MQ=255)
                                                                                                    ttgtttgCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTc                                              <  1:1117622/150‑1 (MQ=255)
                                                                                                      gtttgCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGc                                           <  2:1279818/151‑1 (MQ=255)
                                                                                                         tgCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCtt                                         <  2:799392/150‑1 (MQ=255)
                                                                                                         tgCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCtt                                         <  1:1125831/150‑1 (MQ=255)
                                                                                                              tCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGAtct                                   >  2:1279973/1‑151 (MQ=255)
                                                                                                              tCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGAtct                                   >  2:1180066/1‑151 (MQ=255)
                                                                                                               cGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGAtctc                                  >  1:705743/1‑151 (MQ=255)
                                                                                                                    gCCGTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATctct                                 >  1:1124927/1‑147 (MQ=255)
                                                                                                                       gTCGCGGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTcagcag                          <  2:873816/151‑1 (MQ=255)
                                                                                                                         cgcgGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAgg                         >  1:697137/1‑150 (MQ=255)
                                                                                                                         cgcgGCGTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAgg                         >  1:1165074/1‑150 (MQ=255)
                                                                                                                          gcggcgTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGaa                       <  2:404088/151‑1 (MQ=255)
                                                                                                                          gcggcgTACGCCCCCCAACTGAGGCGCCAGGCCAGCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGaa                       <  1:1032170/151‑1 (MQ=255)
                                                                                                                           cggcgTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAACTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGa                        >  2:902641/1‑149 (MQ=255)
                                                                                                                               gTACGCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCggg                  <  1:385469/151‑1 (MQ=255)
                                                                                                                                   gCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGc                >  2:817891/1‑149 (MQ=255)
                                                                                                                                   gCCCCCCAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAg              <  2:158705/151‑1 (MQ=255)
                                                                                                                                    ccccccAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAgg             <  2:770910/151‑1 (MQ=255)
                                                                                                                                      ccccAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAGGtt           <  2:1283440/151‑1 (MQ=255)
                                                                                                                                       cccAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAGGTTc          <  1:1213762/151‑1 (MQ=255)
                                                                                                                                        ccAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAGGTTCa         <  2:946436/151‑1 (MQ=255)
                                                                                                                                        ccAACTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAGGTTCa         <  2:520010/151‑1 (MQ=255)
                                                                                                                                          aaCTGAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAGGTTCAGc       >  1:1227828/1‑151 (MQ=255)
                                                                                                                                              gAGGCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAGGTTCAGCGTc    <  1:937227/150‑1 (MQ=255)
                                                                                                                                                ggCGCCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAGGTTCAGCGTCTg  <  2:1012221/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTTTTCGTTCCTCAAGCGCATTTAAGGCAGACTGGAGCGCCTCTACGACTTCCGGCTGTGCGGGTTTCGCCGCTTTTGCCTGCTCCAGTTCCTGAGTGATTTGTTTGCTATCGGGGGCCGTCGCGGCGTACGCCCCCCAACTGAGGCACCAGGCCATCAGAAAAGTGATAATCAGGCGCACGTCAGCGTTTCCTTTGATGGATTAGACCTGGTCTTTTTTGTCGTCAACCAATGGGCTGGCGTCGTGTTCTGCTTCGATCTCTTCAGCAGGAAGCGGGGCAGGTTCAGCGTCTG  >  NC_000913/4389190‑4389483

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