Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,477,153 T→A 8.3% K95M (AAG→ATG) ‡ rlmH ← 23S rRNA (pseudouridine(1915)‑N(3))‑methyltransferase RlmH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,477,1530TA 8.3% 927.3 / 50.7 306K95M (AAG→ATG) ‡rlmH23S rRNA (pseudouridine(1915)‑N(3))‑methyltransferase RlmH
Reads supporting (aligned to +/- strand):  ref base T (120/159);  new base A (10/16);  total (130/176)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

AAGTCAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTCGGCG  >  NC_002947/5477004‑5477302
                                                                                                                                                     |                                                                                                                                                     
   tccgcggcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGc                                                                                                                                                    <  1:464180/147‑1 (MQ=255)
   tCAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGc                                                                                                                                                    >  1:5225479/1‑150 (MQ=255)
   tCAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAAGTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGc                                                                                                                                                    <  2:5049396/150‑1 (MQ=255)
     aGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCa                                                                                                                                                  <  1:4190713/150‑1 (MQ=255)
      gcggcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCAt                                                                                                                                                 <  2:2634811/150‑1 (MQ=255)
        ggcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGa                                                                                                                                              <  2:3351745/150‑1 (MQ=255)
         gcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCTCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGa                                                                                                                                              <  2:5202748/150‑1 (MQ=255)
         gcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGa                                                                                                                                              <  1:5656979/150‑1 (MQ=255)
         gcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGa                                                                                                                                              <  2:6074365/150‑1 (MQ=255)
         gcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGa                                                                                                                                              <  1:227441/150‑1 (MQ=255)
              gccGAGCAACGTTGGTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCg                                                                                                                                         >  1:4319794/3‑150 (MQ=255)
              aGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCg                                                                                                                                         >  2:5060043/1‑150 (MQ=255)
                cGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAg                                                                                                                                       <  2:674580/150‑1 (MQ=255)
                cGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAg                                                                                                                                       <  1:1596105/150‑1 (MQ=255)
                 gACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCCAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                      <  2:3635968/150‑1 (MQ=255)
                 gACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                      >  2:5461751/1‑150 (MQ=255)
                 gACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                      <  2:4938875/150‑1 (MQ=255)
                    cAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTg                                                                                                                                   <  2:5845440/150‑1 (MQ=255)
                     aaCGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGa                                                                                                                                  >  1:5347120/1‑150 (MQ=255)
                     aaCGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGa                                                                                                                                  <  1:2882262/150‑1 (MQ=255)
                        gTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACaa                                                                                                                               >  2:4145206/1‑150 (MQ=255)
                            cTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCg                                                                                                                           <  2:6683140/150‑1 (MQ=255)
                             tCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGt                                                                                                                          >  2:3730778/1‑150 (MQ=255)
                             tCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGt                                                                                                                          <  1:4207920/150‑1 (MQ=255)
                                gcgcgcgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTcc                                                                                                                       >  1:6641827/1‑150 (MQ=255)
                                gcgcgcgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTcc                                                                                                                       <  1:4324197/150‑1 (MQ=255)
                                 cgcgcgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTccc                                                                                                                      <  1:1432297/150‑1 (MQ=255)
                                     cgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAgg                                                                                                                  <  1:819549/150‑1 (MQ=255)
                                        gcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTg                                                                                                               >  2:5619635/1‑150 (MQ=255)
                                        gcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTg                                                                                                               >  1:4951504/1‑150 (MQ=255)
                                           aaCCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCAc                                                                                                            <  2:6386968/150‑1 (MQ=255)
                                             ccTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACtt                                                                                                          <  2:2785750/150‑1 (MQ=255)
                                              cTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACttt                                                                                                         >  1:5382968/1‑150 (MQ=255)
                                                   gcgcCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCa                                                                                                    <  2:2349288/150‑1 (MQ=255)
                                                   gcgcCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCa                                                                                                    <  1:1168602/150‑1 (MQ=255)
                                                     gcCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGc                                                                                                  >  2:3502317/1‑150 (MQ=255)
                                                        aGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATg                                                                                               >  2:1524260/1‑150 (MQ=255)
                                                             cTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTc                                                                                          >  2:2029986/1‑150 (MQ=255)
                                                             cTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTc                                                                                          >  1:3776732/1‑150 (MQ=255)
                                                             cTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTc                                                                                          >  2:5649517/1‑150 (MQ=255)
                                                              tCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCg                                                                                         >  2:239833/1‑150 (MQ=255)
                                                              tCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCg                                                                                         >  2:2283499/1‑150 (MQ=255)
                                                                cGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGcc                                                                                       >  2:4443514/1‑150 (MQ=255)
                                                                       cccACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGa                                                                                >  1:3163762/1‑150 (MQ=255)
                                                                       cccACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGa                                                                                <  1:2405995/150‑1 (MQ=255)
                                                                        ccaccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGAc                                                                               >  2:4802059/1‑150 (MQ=255)
                                                                        ccaccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGAc                                                                               >  1:2251777/1‑150 (MQ=255)
                                                                        ccaccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGAc                                                                               >  1:1970679/1‑150 (MQ=255)
                                                                           ccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGa                                                                            <  2:793387/150‑1 (MQ=255)
                                                                            caTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGAt                                                                           >  2:6603337/1‑150 (MQ=255)
                                                                             ctCAAATTCCCCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCGCCTGatccgtc                                                                          <  2:5737581/149‑7 (MQ=255)
                                                                              tCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCATCTCGGTCGCCAGCTTCTCGGTACTCCATTGCTTGCCATGGACCTCGAGTGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCa                                                                         <  1:4834284/150‑1 (MQ=255)
                                                                              tCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCa                                                                         >  1:6807508/1‑150 (MQ=255)
                                                                              tCAAAGTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGGACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCa                                                                         <  2:6138335/150‑1 (MQ=255)
                                                                                aaaTTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAgg                                                                       >  1:6322385/1‑150 (MQ=255)
                                                                                 aaTTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggc                                                                      >  1:1617535/1‑150 (MQ=255)
                                                                                      aCCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACt                                                                                                           <  2:6308489/108‑1 (MQ=255)
                                                                                      aCCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACt                                                                                                           >  1:6308489/1‑108 (MQ=255)
                                                                                      aCCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggcgggcg                                                                 >  1:6798502/1‑150 (MQ=255)
                                                                                        cGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGAc                                                               <  2:2151564/150‑1 (MQ=255)
                                                                                        cGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGAc                                                               <  1:31290/150‑1 (MQ=255)
                                                                                        cGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCGCCTGATCCGTCAGGAGGGCGaa                                                               <  1:200652/150‑2 (MQ=255)
                                                                                        cGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCGCCTGATCCGTCAGGAGGGCGaa                                                               <  2:5137697/150‑2 (MQ=255)
                                                                                        cGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCGCCTGATCCGTCAGGAGGGCGaa                                                               <  2:2609053/150‑2 (MQ=255)
                                                                                                 aTCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCAtt                                                      >  2:4290227/1‑150 (MQ=255)
                                                                                                  tCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTc                                                     <  1:5628061/150‑1 (MQ=255)
                                                                                                   ccAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCt                                                    <  1:2499956/150‑1 (MQ=255)
                                                                                                    cAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCtt                                                   <  1:1842855/150‑1 (MQ=255)
                                                                                                    cAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCtt                                                   <  1:3686993/150‑1 (MQ=255)
                                                                                                    cAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCGGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCtt                                                   <  1:3392844/150‑1 (MQ=255)
                                                                                                      ggCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACACTGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGc                                                 <  1:5780407/150‑1 (MQ=255)
                                                                                                         gcCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCAc                                              <  1:2608040/150‑1 (MQ=255)
                                                                                                           cAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACgg                                            <  2:1263714/150‑1 (MQ=255)
                                                                                                           cAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACgg                                            <  2:4709009/150‑1 (MQ=255)
                                                                                                            aGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACggg                                           <  1:4891325/150‑1 (MQ=255)
                                                                                                            aGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTAAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACggg                                           >  1:1058854/1‑150 (MQ=255)
                                                                                                                  ccAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCa                                     <  1:1299848/150‑1 (MQ=255)
                                                                                                                   cAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAg                                    >  2:1828277/1‑150 (MQ=255)
                                                                                                                     gCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCg                                  <  2:2978939/150‑1 (MQ=255)
                                                                                                                      cTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCgg                                 <  2:3959530/150‑1 (MQ=255)
                                                                                                                      cTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCgg                                 <  1:6144880/150‑1 (MQ=255)
                                                                                                                        cGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGa                               <  1:6206645/150‑1 (MQ=255)
                                                                                                                          gTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGAtt                             <  1:5447964/150‑1 (MQ=255)
                                                                                                                           tCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGAttt                            >  1:2100773/1‑150 (MQ=255)
                                                                                                                             gCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTcc                          >  1:3269527/1‑150 (MQ=255)
                                                                                                                                      tCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTcc                 >  2:3109910/1‑150 (MQ=255)
                                                                                                                                       tcgaggtCCATGGCAAGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCa                >  1:2750949/7‑150 (MQ=255)
                                                                                                                                        ggTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAg               <  1:3918113/150‑1 (MQ=255)
                                                                                                                                        ggTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAg               <  1:605627/150‑1 (MQ=255)
                                                                                                                                         gaggtCCATGGCAAGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGc              <  1:6551105/146‑1 (MQ=255)
                                                                                                                                         gTACTCCATGGCTGGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGc              <  2:4596304/150‑1 (MQ=255)
                                                                                                                                            cTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGAc           >  1:345787/1‑150 (MQ=255)
                                                                                                                                             tCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACa          <  1:6656787/150‑1 (MQ=255)
                                                                                                                                             tCCATGGCAAGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACa          >  1:5899239/1‑150 (MQ=255)
                                                                                                                                             tCCATGGCAAGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACa          >  1:5119681/1‑150 (MQ=255)
                                                                                                                                              ccATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAg         <  2:1204279/150‑1 (MQ=255)
                                                                                                                                              ccATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAg         <  1:1150548/150‑1 (MQ=255)
                                                                                                                                               cATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGc        >  2:775557/1‑150 (MQ=255)
                                                                                                                                                   gCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTcgg    >  1:1703116/1‑150 (MQ=255)
                                                                                                                                                     ttGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTcggcg  <  1:5060043/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGTCAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTCGGCG  >  NC_002947/5477004‑5477302

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: