Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,414,841 G→C 7.3% intergenic (+140/‑281) rlmM → / → PP_2114 23S rRNA (cytidine(2498)‑2'‑O)‑methyltransferase RlmM/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,414,8410GC 7.3% 648.5 / 33.0 205intergenic (+140/‑281)rlmM/PP_211423S rRNA (cytidine(2498)‑2'‑O)‑methyltransferase RlmM/transposase
Reads supporting (aligned to +/- strand):  ref base G (99/91);  new base C (12/3);  total (111/94)
Fisher's exact test for biased strand distribution p-value = 5.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 215 total aligned reads displayed.

GCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACA  >  NC_002947/2414692‑2414990
                                                                                                                                                      |                                                                                                                                                       
gCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTGAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATg                                                                                                                                                         >  1:2378377/1‑150 (MQ=255)
gCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATg                                                                                                                                                         >  1:6055953/1‑150 (MQ=255)
 ccACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATgg                                                                                                                                                        <  2:4247727/150‑1 (MQ=255)
 ccACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATgg                                                                                                                                                        >  1:4406228/1‑150 (MQ=255)
 ccACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATgg                                                                                                                                                        >  1:417998/1‑150 (MQ=255)
  cACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGa                                                                                                                                                       >  1:7262426/1‑150 (MQ=255)
  cACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGa                                                                                                                                                       <  1:6329908/150‑1 (MQ=255)
    cGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATg                                                                                                                                                     <  1:5674715/150‑1 (MQ=255)
       tGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGccc                                                                                                                                                  >  2:4458680/1‑150 (MQ=255)
         aGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGt                                                                                                                                                <  1:7012138/150‑1 (MQ=255)
          ggCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGACATGGATGCCCGTc                                                                                                                                               <  2:2364129/150‑1 (MQ=255)
           gCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCa                                                                                                                                              <  1:6519227/150‑1 (MQ=255)
            ccAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTCac                                                                                                                                             <  1:3014972/150‑2 (MQ=255)
             cAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAg                                                                                                                                            >  1:9101879/1‑150 (MQ=255)
                cACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAgcgc                                                                                                                                         <  1:1847522/150‑1 (MQ=255)
                 aCGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                                                                                                                        <  1:583176/150‑1 (MQ=255)
                   gtcgtTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                     <  2:3846893/150‑1 (MQ=255)
                   gtcgtTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                     >  1:6858659/1‑150 (MQ=255)
                   gtcgtTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                     <  2:4687519/150‑1 (MQ=255)
                       ttctgcttctgcGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGcc                                                                                                                                 >  2:3943789/1‑150 (MQ=255)
                       ttctgcttctgcGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGcc                                                                                                                                 >  1:6599596/1‑150 (MQ=255)
                       ttctgcttctgcGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGcc                                                                                                                                 >  2:4109561/1‑150 (MQ=255)
                        ggcgagtcactttttgtcaaccgcgacaaaaagtaaccaaaaaacgctgcgctcccatcatccggcccctgcgctgcgctccggggtcccctcactccGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGccc                                                                                                                                >  2:7925395‑M1/99‑150 (MQ=255)
                             ttctgcGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatgg                                                                                                                           >  1:5800705/1‑150 (MQ=255)
                               ctgcgtctgcgtCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGTGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggat                                                                                                                         >  1:5131141/1‑150 (MQ=255)
                               ctgcgtctgcgtCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggat                                                                                                                         <  2:3223153/150‑1 (MQ=255)
                                 gcgtctgcgtctGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggatgg                                                                                                                       <  1:3762918/150‑1 (MQ=255)
                                  cgtctgcgtctgCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATggg                                                                                                                      <  1:6210784/150‑1 (MQ=255)
                                  cgtctgcgtctgCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATggg                                                                                                                      <  2:6294146/150‑1 (MQ=255)
                                  cgtctgcgtctgCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGGATggg                                                                                                                      >  1:3058470/1‑150 (MQ=255)
                                   gtctgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGAGgggg                                                                                                                     >  2:6049021/1‑150 (MQ=255)
                                    tctgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGc                                                                                                                    >  1:3264113/1‑150 (MQ=255)
                                       gcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTg                                                                                                                 <  1:1794080/150‑1 (MQ=255)
                                       gcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCATGGGGCCtc                                                                                                                 <  1:4729440/150‑2 (MQ=255)
                                       gcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCATGGGGCCtc                                                                                                                 <  2:1999697/150‑2 (MQ=255)
                                        cgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGc                                                                                                                <  2:65419/150‑1 (MQ=255)
                                         gtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCATGGGGCCtca                                                                                                                <  1:9018372/150‑3 (MQ=255)
                                            acgcgacaaaaagtaaccaaaaaacgctgcgctcctatcatccggcccctgcgctacgctccggggttccctcgctccGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatgg                                                                                                                           <  2:8273847‑M1/57‑1 (MQ=255)
                                              cttctgcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCATGGGGCCTCAgcgct                                                                                                          <  2:1322135/150‑2 (MQ=14)
                                                tctgcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTa                                                                                                       <  1:1460030/150‑1 (MQ=255)
                                                tctgcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTa                                                                                                       <  2:2145013/150‑1 (MQ=255)
                                                    cttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTcc                                                                                                   >  1:7563979/1‑150 (MQ=40)
                                                     ttttGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTccc                                                                                                  >  2:8922519/1‑150 (MQ=40)
                                                       ttGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCgg                                                                                                <  2:5401562/150‑1 (MQ=38)
                                                       aagtaaccaaaaaacgctggctcctatcatccggcccctgcgctacgctccggggttccctcgctccGGGCTTGCTCCCGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCgg                                                                                                >  1:1101028‑M1/68‑150 (MQ=255)
                                                         gCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGa                                                                                              <  2:8095162/150‑1 (MQ=37)
                                                             ccaaaaaacgctggctcctatcatccggcccctgcgctacgctccgggggtccctcgctccGGGCTTGCTCCCGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCaa                                                                                          >  1:6046528‑M1/62‑150 (MQ=255)
                                                              aaaaaacgctgcgctcccatcatccggcccctgcgctgcgctccggggtcccctcactccGGCCTTGCGAAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAg                                                                                         >  2:146547‑M1/61‑150 (MQ=255)
                                                                 ttCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGccc                                                                                      >  1:8196240/1‑150 (MQ=37)
                                                                 ttCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGccc                                                                                      >  2:2562942/1‑150 (MQ=37)
                                                                   accgtcggctcccatcatccggcccctacgctacgctccggggtcccctcactctGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                                                   >  1:2968202‑M1/56‑150 (MQ=255)
                                                                    tAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGa                                                                                   <  1:368202/150‑1 (MQ=38)
                                                                       gcgcgcGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGAGCGTACTCCCGGGAGCAAGCCCGGAGCg                                                                                <  1:378945/150‑1 (MQ=38)
                                                                        tggctcctatcatccggcccctgcgctacgctccggggttccctcgctcgGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAg                                                                              >  2:4657222‑M1/51‑150 (MQ=255)
                                                                         gcgcgATAGTTGAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGTAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAg                                                                              <  2:3498168/150‑1 (MQ=21)
                                                                          cgcgATCGTTCCGGCGCCGCAGATTGCGACTTCAGGAGGCCGAGCGTCGGGCTTGCGAAGGGAGGTGACGGGCATGGAGGCCCGTCAAGCGCTGAGGCCCCAGGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAgg                                                                             <  2:3055385/150‑1 (MQ=11)
                                                                          cgcgATCGTTCCGGCGACGCAGATTGCGACTTCCGGTGGCCGCGCGTTGGGCTTGGGACGGGAGGTGACGGGCTTGGATGCCCGTCAAGCGCTGTGGCCCCATGGATGGGGCCTGCTGCGCTTCCTCCCGGGAGCAAGCCCGGAGCGAgg                                                                             <  2:1548737/150‑1 (MQ=11)
                                                                           gcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAggg                                                                            >  1:84759/1‑150 (MQ=37)
                                                                                   tCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCgg                                                                    >  1:8937490/1‑150 (MQ=37)
                                                                                         gACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTa                                                              <  2:4406228/150‑1 (MQ=37)
                                                                                         gACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTa                                                              <  2:2637519/150‑1 (MQ=37)
                                                                                           cGCAGATTGCGACTTCCGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCCTGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAgc                                                            <  2:5053984/150‑1 (MQ=255)
                                                                                            gCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAgcg                                                           >  2:8292273/1‑150 (MQ=37)
                                                                                              cttttGCGACTCCTGGTGGCCGAGCTTTGGGCTTGCGTTGGGTGGTGAGGGGCATGGTTGCCCGTCATGCGCTGCGGCCCCCTGTTTGGGGCCTGCAGCGCTTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCa                                                         <  2:281389/147‑1 (MQ=11)
                                                                                              aGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGACCCCCGGAGCGTAGCGCa                                                         >  1:1584807/1‑150 (MQ=255)
                                                                                                aTTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGTGGTGACGGGCATGGAGGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAgg                                                       <  1:8344069/150‑1 (MQ=255)
                                                                                                  tGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAgggg                                                     >  2:1831681/1‑150 (MQ=37)
                                                                                                     gACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCg                                                  >  1:7450712/1‑150 (MQ=37)
                                                                                                        ttCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgat                                               >  1:7187504/1‑150 (MQ=37)
                                                                                                           aggaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatgat                                            >  1:6270408/1‑150 (MQ=37)
                                                                                                           aggaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatgat                                            >  2:1063523/1‑150 (MQ=37)
                                                                                                            ggaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATa                                           <  2:7285620/150‑1 (MQ=38)
                                                                                                              aggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAgg                                         <  2:6858659/150‑1 (MQ=37)
                                                                                                               ggCCGAGCGTAGGGCTTGCGACGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGa                                        <  2:7462847/150‑1 (MQ=37)
                                                                                                               ggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGa                                        <  2:3743964/150‑1 (MQ=37)
                                                                                                                 ccGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAga                                      >  1:1030845/1‑149 (MQ=37)
                                                                                                                           GGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttt                            <  1:7777185‑M1/150‑1 (MQ=255)
                                                                                                                           GGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttt                            <  2:7163729‑M1/150‑1 (MQ=255)
                                                                                                                            ggtggCGAGGGGTGGTGCGGGGCTTGGCGGCCCGCCACGCGCGGCGGCCCCTTGGTGGGGGCCGGCCGCGCGTACCCCGGGGACCACGCCGGGGGCGCGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGttttt                           <  1:5553562/146‑1 (MQ=9)
                                                                                                                             CTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGACCCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttttt                          <  2:6213800‑M1/150‑1 (MQ=255)
                                                                                                                             CTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttttt                          >  2:5243701‑M1/1‑150 (MQ=255)
                                                                                                                              TTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGTATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTg                         >  2:5784564‑M1/1‑150 (MQ=255)
                                                                                                                              TTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTg                         >  2:8024262‑M1/1‑150 (MQ=255)
                                                                                                                                   gacgggAGGTGACGGGCATGGATGCCCGTCACGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTa                     <  2:8572971/147‑1 (MQ=255)
                                                                                                                                      gggAGGTGACGGGCAGGGATGCCCGCCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCACGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACtt                  <  2:7828923/150‑1 (MQ=255)
                                                                                                                                          ggTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTg              >  1:1334605/1‑150 (MQ=39)
                                                                                                                                           ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGt             >  1:6751305/2‑150 (MQ=255)
                                                                                                                                           ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGt             >  1:2445877/2‑150 (MQ=255)
                                                                                                                                           ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCAGGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGt             >  2:8462255/2‑150 (MQ=255)
                                                                                                                                             gACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcg           >  2:2492610/1‑150 (MQ=39)
                                                                                                                                             gACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcg           >  2:5744219/1‑150 (MQ=255)
                                                                                                                                              aCGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcgc          >  1:8924422/1‑150 (MQ=255)
                                                                                                                                              aCGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcgc          >  2:2760303/1‑150 (MQ=255)
                                                                                                                                                  gCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGttt      >  2:8289955/1‑150 (MQ=39)
                                                                                                                                                  gCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGgtt      >  2:6449953/1‑150 (MQ=255)
                                                                                                                                                  gCATGGATGCCCGACAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGttt      <  2:4779137/150‑1 (MQ=255)
                                                                                                                                                   cATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTg     >  2:7699401/1‑150 (MQ=39)
                                                                                                                                                   cATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTg     >  1:6400065/1‑150 (MQ=39)
                                                                                                                                                    aTGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGa    <  2:6065297/150‑1 (MQ=39)
                                                                                                                                                      ggATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACa  >  2:6685467/1‑150 (MQ=39)
                                                                                                                                                      |                                                                                                                                                       
GCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACA  >  NC_002947/2414692‑2414990

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