Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,216,463 G→A 5.4% R8C (CGC→TGC)  PP_3693 ← transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,216,4630GA 5.4% 949.2 / 26.2 293R8C (CGC→TGC) PP_3693transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base G (128/149);  new base A (8/8);  total (136/157)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 308 total aligned reads displayed.

AATATTGATTACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGAGCAGGGATTTGACCGACCTCC  >  NC_002947/4216314‑4216612
                                                                                                                                                     |                                                                                                                                                     
 atatTGATTACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCgg                                                                                                                                                      <  2:1371362/150‑1 (MQ=255)
  tatTGATTACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGa                                                                                                                                                     <  2:3331028/150‑1 (MQ=255)
   atTGATTACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGaa                                                                                                                                                    <  2:6007319/150‑1 (MQ=255)
   atTGATTACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGaa                                                                                                                                                    <  1:3882329/150‑1 (MQ=255)
     tGATTACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAAtt                                                                                                                                                  <  1:960738/150‑1 (MQ=255)
       atgaCATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCg                                                                                                                                                <  2:6923088/147‑1 (MQ=255)
       aTTACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCg                                                                                                                                                <  2:6578375/150‑1 (MQ=255)
       aTTACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCg                                                                                                                                                <  1:6080480/150‑1 (MQ=255)
        ttACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCgg                                                                                                                                               <  2:6813580/150‑1 (MQ=255)
        ttACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCgg                                                                                                                                               >  1:3045752/1‑150 (MQ=255)
        ttACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCAGAATTCgg                                                                                                                                               >  2:7170107/1‑150 (MQ=255)
           cATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCa                                                                                                                                            <  1:5721362/150‑1 (MQ=255)
               aaGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACg                                                                                                                                        >  1:4874366/1‑150 (MQ=255)
               aaGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACg                                                                                                                                        >  2:3988083/1‑150 (MQ=255)
                 gCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGGAGGTTGGTTTCGTATTCGATTTCGTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAATGATTGTTCGAGCTTGCGGTATTCGGCGAgaggag                                                                                                                                      >  2:5645837/1‑150 (MQ=255)
                  cAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACgaga                                                                                                                                     >  1:4651617/1‑150 (MQ=255)
                     gCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACAt                                                                                                                                  <  2:1903795/150‑1 (MQ=255)
                       tGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTa                                                                                                                                <  1:2525239/150‑1 (MQ=255)
                         cACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGt                                                                                                                              <  2:6215206/150‑1 (MQ=255)
                         cACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTGTGCAGCTCGGAGCTGCCCTTCATGGCTTCATGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGt                                                                                                                              <  1:2185236/150‑1 (MQ=255)
                          aCCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTg                                                                                                                             <  1:188458/150‑1 (MQ=255)
                           cccATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGa                                                                                                                            <  2:7061180/150‑1 (MQ=255)
                           cccATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGa                                                                                                                            >  1:199258/1‑150 (MQ=255)
                            ccATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGAt                                                                                                                           >  2:4157940/1‑150 (MQ=255)
                            ccATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGAt                                                                                                                           >  1:6617150/1‑150 (MQ=255)
                             cATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATc                                                                                                                          <  1:6500194/150‑1 (MQ=255)
                             cATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATc                                                                                                                          <  1:5625087/150‑1 (MQ=255)
                             cATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATc                                                                                                                          <  1:2289354/150‑1 (MQ=255)
                             cATACTCAGCAAGGAGGTTCCGCAGCTTGGTTACGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATc                                                                                                                          <  2:5518956/150‑1 (MQ=255)
                                aCTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTc                                                                                                                       <  1:2933443/150‑1 (MQ=255)
                                aCTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTc                                                                                                                       >  2:3135216/1‑150 (MQ=255)
                                  tCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCTGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCaa                                                                                                                     >  1:3168981/1‑150 (MQ=255)
                                  tCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCaa                                                                                                                     <  2:4753382/150‑1 (MQ=255)
                                  tCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTCCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCaa                                                                                                                     >  1:2477177/1‑150 (MQ=255)
                                     gCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGa                                                                                                                  <  2:4537895/150‑1 (MQ=255)
                                       aaGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCATGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCAGAATTCGGGCAGACGAGACATTAGTGATCCTCAATGAtt                                                                                                                >  1:7125361/1‑150 (MQ=255)
                                       aaGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGAtt                                                                                                                >  2:2624545/1‑150 (MQ=255)
                                        aggaggTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGAttt                                                                                                               >  2:3768394/1‑150 (MQ=255)
                                         ggaggTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTa                                                                                                              >  1:2213875/1‑150 (MQ=255)
                                         ggaggTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTa                                                                                                              >  2:4263199/1‑150 (MQ=255)
                                         ggaggTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTa                                                                                                              <  1:1074134/150‑1 (MQ=255)
                                           aggTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGa                                                                                                            <  2:5842405/150‑1 (MQ=255)
                                                ccGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGc                                                                                                       <  1:5489936/150‑1 (MQ=255)
                                                ccGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGc                                                                                                       <  2:1292606/150‑1 (MQ=255)
                                                         ggTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAAtt                                                                                              <  2:2668823/150‑1 (MQ=255)
                                                          gTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTg                                                                                             >  2:5604696/1‑150 (MQ=255)
                                                           tttCGAATTCGATTTCCTTTTGCTGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGa                                                                                            >  1:2110345/1‑150 (MQ=255)
                                                           tttCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGa                                                                                            <  2:1416903/150‑1 (MQ=255)
                                                           tttCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGa                                                                                            <  1:4509048/150‑1 (MQ=255)
                                                           tttCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGa                                                                                            <  2:2360315/150‑1 (MQ=255)
                                                           tttCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGa                                                                                            <  1:3768394/150‑1 (MQ=255)
                                                           tttCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAg                                                                                           >  2:3680414/1‑150 (MQ=255)
                                                             tCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGt                                                                                          <  2:5745824/150‑1 (MQ=255)
                                                             tCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGt                                                                                          >  1:6124858/1‑150 (MQ=255)
                                                                aaTTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATc                                                                                       <  1:6871942/150‑1 (MQ=255)
                                                                  ttCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCt                                                                                     >  2:6439935/1‑150 (MQ=255)
                                                                  ttCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCt                                                                                     >  2:4823560/1‑150 (MQ=255)
                                                                   tCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTa                                                                                    >  1:4341702/1‑150 (MQ=255)
                                                                       tttCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAAtt                                                                                              <  2:5178215/136‑1 (MQ=255)
                                                                       tttCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAAtt                                                                                              >  1:5178215/1‑136 (MQ=255)
                                                                       tttCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAAcc                                                                                >  2:3794926/1‑150 (MQ=255)
                                                                       tttCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAAcc                                                                                >  1:6351972/1‑150 (MQ=255)
                                                                          ccTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAgg                                                                             >  2:1799867/1‑150 (MQ=255)
                                                                          ccTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAgg                                                                             >  1:1051137/1‑150 (MQ=255)
                                                                            ttttGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGaa                                                                           >  2:123549/1‑150 (MQ=255)
                                                                            ttttGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGaa                                                                           <  2:3018511/150‑1 (MQ=255)
                                                                            ttttGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGaa                                                                           >  1:1661020/1‑150 (MQ=255)
                                                                            ttttGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGaa                                                                           <  1:273116/150‑1 (MQ=255)
                                                                                  aGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGAc                                                                                                                                         >  2:1644447/1‑82 (MQ=255)
                                                                                       ggAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTg                                                                <  2:782722/150‑1 (MQ=255)
                                                                                             gCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGt                                                          <  1:4236707/150‑1 (MQ=255)
                                                                                               ccTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTg                                                        <  2:2884828/150‑1 (MQ=255)
                                                                                               ccTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTg                                                        >  1:5415463/1‑150 (MQ=255)
                                                                                               ccTTCAGGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTg                                                        <  2:1518346/150‑1 (MQ=255)
                                                                                                cTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTgg                                                       >  2:5125851/1‑150 (MQ=255)
                                                                                                 ttCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGt                                                      >  2:4676836/1‑150 (MQ=255)
                                                                                                    aTGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTtc                                                   >  1:4153105/1‑150 (MQ=255)
                                                                                                    aTGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAAACAGGAAGCAAGGTTCTGCTACGTTGGtttg                                                   >  1:3948606/1‑149 (MQ=255)
                                                                                                     tGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTtct                                                  >  1:3415259/1‑150 (MQ=255)
                                                                                                     tGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCAGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTtct                                                  <  1:2796630/150‑1 (MQ=255)
                                                                                                       gCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTctct                                                <  1:872140/150‑1 (MQ=255)
                                                                                                         ttCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCttt                                              >  2:298670/1‑150 (MQ=255)
                                                                                                         ttCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCttt                                              <  2:2808508/150‑1 (MQ=255)
                                                                                                         ttCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCttt                                              >  1:2586520/1‑150 (MQ=255)
                                                                                                             aGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCttttttt                                          <  1:1349552/150‑1 (MQ=255)
                                                                                                                   gCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTg                                    >  2:2095365/1‑150 (MQ=255)
                                                                                                                   gCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTg                                    >  2:4605628/1‑150 (MQ=255)
                                                                                                                     gAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGaa                                  <  1:7041681/150‑1 (MQ=255)
                                                                                                                           cgcAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCg                            >  1:2065470/1‑150 (MQ=255)
                                                                                                                           cgcAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCg                            >  1:1168494/1‑150 (MQ=255)
                                                                                                                                ccAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAGGCTTCGCTCGa                       >  1:5800659/1‑150 (MQ=255)
                                                                                                                                ccAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGa                       >  1:5030457/1‑150 (MQ=255)
                                                                                                                                 cAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGAg                      >  1:1458333/1‑150 (MQ=255)
                                                                                                                                  aaaGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTAGAAGCTTCGCTCGAGc                     >  2:7208093/1‑150 (MQ=255)
                                                                                                                                    aGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGAGCAg                   <  1:7312384/150‑1 (MQ=255)
                                                                                                                                            ccAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGAGCAGGGATTTga           <  2:1735145/150‑1 (MQ=255)
                                                                                                                                            ccAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGAGCAGGGATTTga           <  2:6354442/150‑1 (MQ=255)
                                                                                                                                                 ttGCAGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGAGCAGGGATTTgaccgac      >  1:6262253/1‑150 (MQ=255)
                                                                                                                                                   gCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGAGCAGGGATTTGACCGAcct    <  2:4862838/150‑1 (MQ=255)
                                                                                                                                                     ggAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGAGCAGGGATTTGACCGAcctcc  <  2:6334102/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATATTGATTACATCAAGCAGGCTGCACCCATACTCAGCAAGGAGGTTCCGCAGCTTGGTTTCGAATTCGATTTCCTTTTGCAGCTCGGAGCTGCCCTTCATGGCTTCAAGCTCAGCGAGCTGCGCAGCCAAAGATTGTTCCAGCTTGCGGAATTCGGCCAGACGAGACATTAGTGATCCTCAATGATTTAGATGGGCATGTAAATTGAGTATCCTAAACCAGGAAGCAAGGTTCTGCTACGTTGGTTTCTCTTTTTTTGACTTGAAGCTTCGCTCGAGCAGGGATTTGACCGACCTCC  >  NC_002947/4216314‑4216612

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