Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,205,660 T→A 5.4% Y355* (TAT→TAA yehL → putative AAA(+) MoxR family ATPase YehL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,205,6600TA 5.4% 486.7 / 11.2 147Y355* (TAT→TAAyehLputative AAA(+) MoxR family ATPase YehL
Reads supporting (aligned to +/- strand):  ref base T (72/67);  new base A (2/6);  total (74/73)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 156 total aligned reads displayed.

TTAGCGAATCGCGCGGGGGAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGAAGGCCCGGCTG  >  NC_000913/2205511‑2205809
                                                                                                                                                     |                                                                                                                                                     
 tAGCGAATCGCGCGGGGGAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATc                                                                                                                                                      <  1:1375894/150‑1 (MQ=255)
  aGCGAATCGCGCGGGGGAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGGCTGCTCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCa                                                                                                                                                     <  1:2379923/150‑1 (MQ=255)
        tCGCGCGGGGGAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAGGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCg                                                                                                                                               <  2:391068/150‑1 (MQ=255)
         cgcgcgGGGGAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATc                                                                                                                                                      >  1:1364341/1‑142 (MQ=255)
          gcgcgGGGGAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCGGGCTTGATAATAAGCCTGcc                                                                                                                                             <  1:500407/150‑1 (MQ=255)
          gcgcgGGGGAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGcc                                                                                                                                             >  2:2373740/1‑150 (MQ=255)
              gggggAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCg                                                                                                                                         <  1:990700/150‑1 (MQ=255)
                gggAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCGGGCTTGATAATAAGCCTGCCagtgag                                                                                                                                       <  2:1004801/150‑6 (MQ=255)
                gggAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCGGGCTTGATAATAAGCCTGCCagtgag                                                                                                                                       <  2:1404176/150‑6 (MQ=255)
                gggAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGcc                                                                                                                                       <  2:1536582/150‑1 (MQ=255)
                gggAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCGGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGcc                                                                                                                                       <  2:516857/150‑1 (MQ=255)
                  gAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTg                                                                                                                                     <  2:934509/150‑1 (MQ=255)
                  gAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTg                                                                                                                                     <  2:2417850/150‑1 (MQ=255)
                    gCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGcc                                                                                                                                   >  1:291155/1‑150 (MQ=255)
                    gCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGcc                                                                                                                                   >  1:1817652/1‑150 (MQ=255)
                     ccTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCg                                                                                                                                  <  1:1104947/150‑1 (MQ=255)
                     ccTGCCGACCTGGTGGAATGTATTGCCGCAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCg                                                                                                                                  >  2:214189/1‑150 (MQ=255)
                      cTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGt                                                                                                                                 >  2:1853409/1‑150 (MQ=255)
                       tGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTg                                                                                                                                <  1:1601090/150‑1 (MQ=255)
                       tGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTg                                                                                                                                <  1:1225715/150‑1 (MQ=255)
                        gCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGa                                                                                                                               >  2:1669548/1‑150 (MQ=255)
                               tggtggAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGCTCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAgc                                                                                                                        >  2:1422962/1‑149 (MQ=255)
                                ggtggAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGAt                                                                                                                       >  2:1605565/1‑150 (MQ=255)
                                 gtggAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATg                                                                                                                      <  2:1114100/150‑1 (MQ=255)
                                 gtggAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATg                                                                                                                      >  1:516231/1‑150 (MQ=255)
                                 gtggAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATg                                                                                                                      <  1:1483888/150‑1 (MQ=255)
                                  tggAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGc                                                                                                                     >  1:985997/1‑150 (MQ=255)
                                   ggAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCGAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCa                                                                                                                    >  2:576644/1‑150 (MQ=255)
                                   ggAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCa                                                                                                                    <  1:913808/150‑1 (MQ=255)
                                     aaTGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATg                                                                                                                  <  1:783585/150‑1 (MQ=255)
                                        gTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATgagc                                                                                                               <  1:543508/150‑1 (MQ=255)
                                          aTTGCCGGAACCATCGTCAAAGATAACTAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGa                                                                                                                        <  2:537544/139‑1 (MQ=255)
                                             gCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGcc                                                                                                          >  2:331648/1‑150 (MQ=255)
                                             gCCGGAACCATCGACAAAGATAACGAACAAGATCGCGCGAGTCTGAGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTCAGGAAAGATGCATGAGAGAGcc                                                                                                          >  2:82853/1‑150 (MQ=255)
                                             gCCGGAACCATCATCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGcc                                                                                                          >  1:2367706/1‑150 (MQ=255)
                                              ccGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCg                                                                                                         <  1:1384622/150‑1 (MQ=255)
                                                 gAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTa                                                                                                      <  1:93412/150‑1 (MQ=255)
                                                  aaCCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTaa                                                                                                     <  1:184867/150‑1 (MQ=255)
                                                   aCCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAAt                                                                                                    <  1:1645544/150‑1 (MQ=255)
                                                   aCCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAAt                                                                                                    <  1:32040/150‑1 (MQ=255)
                                                   aCCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAAt                                                                                                    <  1:2422202/150‑1 (MQ=255)
                                                     cATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGcc                                                                                                          <  1:1436032/142‑1 (MQ=255)
                                                     cATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTg                                                                                                  >  1:1765401/1‑150 (MQ=255)
                                                        cGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcg                                                                                               <  1:1026715/150‑1 (MQ=255)
                                                         gTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcgg                                                                                              <  1:135751/150‑1 (MQ=255)
                                                           cAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcggca                                                                                            <  2:1116996/150‑1 (MQ=255)
                                                             aaGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcggcatc                                                                                          >  1:562065/1‑150 (MQ=255)
                                                                aTAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcggcatccgg                                                                                       >  1:189355/1‑150 (MQ=255)
                                                                aTAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcggcatccgg                                                                                       <  2:1673638/150‑1 (MQ=255)
                                                                  aaCGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcggcatccggca                                                                                     <  2:255934/150‑1 (MQ=255)
                                                                   aCGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcggcatccggcat                                                                                    >  1:1958951/1‑150 (MQ=255)
                                                                   aCGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcggcatccggcat                                                                                    >  2:2449978/1‑150 (MQ=255)
                                                                    cGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTcggcatccggcatc                                                                                   >  1:701809/1‑150 (MQ=255)
                                                                      aagaagATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGcatcat                                                                                 <  2:2349372/150‑1 (MQ=255)
                                                                          agaTCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTc                                                                             <  2:124705/150‑1 (MQ=255)
                                                                                cgcgcgTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCt                                                                       >  1:115136/1‑150 (MQ=255)
                                                                                        tGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCgg                                                               >  2:276668/1‑150 (MQ=255)
                                                                                          cgcCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCCCCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCt                                                             >  2:2297081/1‑150 (MQ=255)
                                                                                          cgcCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCt                                                             <  2:217016/150‑1 (MQ=255)
                                                                                             cGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGt                                                          <  1:2101212/150‑1 (MQ=255)
                                                                                              gTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCATGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTg                                                         >  1:1825762/1‑150 (MQ=255)
                                                                                                 aCTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGaaa                                                      >  1:2282680/1‑150 (MQ=255)
                                                                                                  cTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGTCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAAt                                                     <  2:562065/150‑1 (MQ=255)
                                                                                                  cTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAAt                                                     >  1:330976/1‑150 (MQ=255)
                                                                                                  cTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAAt                                                     >  1:412181/1‑150 (MQ=255)
                                                                                                    tCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCg                                                   >  1:1270444/1‑150 (MQ=255)
                                                                                                      gAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGTATGATAGAGCATTTAATTGTCGGCATACGTAATCATATTCAGTCCTTCTCCCTGCTGGTGTAATCGtt                                                 <  2:865732/150‑1 (MQ=255)
                                                                                                       aaCAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTa                                                <  1:165941/150‑1 (MQ=255)
                                                                                                           gCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCg                                            >  2:2274708/1‑150 (MQ=255)
                                                                                                               gtCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAg                                        <  1:1938817/150‑1 (MQ=255)
                                                                                                               gtCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAg                                        >  1:22647/1‑150 (MQ=255)
                                                                                                                tCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGc                                       <  2:1876798/150‑1 (MQ=255)
                                                                                                                 cGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGcc                                      <  1:1012049/150‑1 (MQ=255)
                                                                                                                 cGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGcc                                      >  1:1182444/1‑150 (MQ=255)
                                                                                                                 cGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGcc                                      >  1:2487469/1‑150 (MQ=255)
                                                                                                                  gCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCa                                     >  1:2372297/1‑150 (MQ=255)
                                                                                                                    cacacaTAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGc                                   >  2:769816/1‑150 (MQ=255)
                                                                                                                     acacaTAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCg                                  <  1:1468477/150‑1 (MQ=255)
                                                                                                                       acaTAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGc                                >  2:426061/1‑150 (MQ=255)
                                                                                                                       acaTAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGc                                >  2:1153244/1‑150 (MQ=255)
                                                                                                                           aaaGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACg                            <  1:214189/150‑1 (MQ=255)
                                                                                                                           aaaGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACg                            >  2:2005102/1‑150 (MQ=255)
                                                                                                                           aaaGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGAGAGCAAAGATCCATGACCGCGCCGATAATTGACGCCATCCGGCATCATACTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGCAAGCCAGCGCCCAcc                            >  2:1654955/1‑149 (MQ=255)
                                                                                                                             agaagCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGAt                          >  2:1442633/1‑150 (MQ=255)
                                                                                                                               aagCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATAc                        >  1:1874802/1‑150 (MQ=255)
                                                                                                                                    cACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGtt                   <  1:1393951/150‑1 (MQ=255)
                                                                                                                                      gtGGCGGGCTTGATAATAAGCCTGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGttga                 >  2:1294140/2‑150 (MQ=255)
                                                                                                                                      gtGGCGGGCTTGATAATAAGCCTGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGACTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGttga                 >  1:128716/2‑150 (MQ=255)
                                                                                                                                      cTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGttga                 <  1:733988/150‑1 (MQ=255)
                                                                                                                                       tGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGttgat                <  1:1304961/150‑1 (MQ=255)
                                                                                                                                        ggCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGttgatt               >  2:1915809/1‑150 (MQ=255)
                                                                                                                                           aGGCTTATTATCAAGCCCGCCACCGCCTGCCATGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGaa            <  2:1825762/150‑1 (MQ=255)
                                                                                                                                            ggCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGAAg           <  2:701809/150‑1 (MQ=255)
                                                                                                                                            ggCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGAAg           >  2:1001985/1‑150 (MQ=255)
                                                                                                                                              cTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGAAGGc         >  2:2131089/1‑150 (MQ=255)
                                                                                                                                               ttattaTCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGAATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGAAGGcc        >  2:1963066/1‑150 (MQ=255)
                                                                                                                                                  ttatCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGAAGGCCCgg     >  2:1447157/1‑150 (MQ=255)
                                                                                                                                                  ttatCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGAAGGCCCgg     <  1:331648/150‑1 (MQ=255)
                                                                                                                                                    atCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGAAGGCCCGGCt   >  2:2273875/1‑150 (MQ=255)
                                                                                                                                                     tCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCCGCGGCCACGATACGTGTTGATTGAAGGCCCGGCTg  >  1:1886120/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTAGCGAATCGCGCGGGGGAGCCTGCCGACCTGGTGGAATGTATTGCCGGAACCATCGTCAAAGATAACGAAGAAGATCGCGCGCGTCTGCGCCGTTACTTCGAACAGCGTGTCGCCACACATAAAGAAGCTCACTGGCAGGCTTATTATCAAGCCCGCCACCGCCTGCCGTGAGGAAAGATGCATGAGCGAGCCGTTAATTGTCGGCATCCGGCATCATAGTCCGGCCTGCGCCCGGCTGGTGAAATCGTTAATCGAAAGCCAGCGGCCACGATACGTGTTGATTGAAGGCCCGGCTG  >  NC_000913/2205511‑2205809

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