Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,477,153 T→A 5.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 5.3% 742.1 / 22.2 227K95M (AAG→ATG) ‡rlmH23S rRNA (pseudouridine(1915)‑N(3))‑methyltransferase RlmH
Reads supporting (aligned to +/- strand):  ref base T (92/118);  new base A (5/7);  total (100/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 229 total aligned reads displayed.

AAGTCAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTCGGCG  >  NC_002947/5477004‑5477302
                                                                                                                                                     |                                                                                                                                                     
aaGTCAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCt                                                                                                                                                       >  2:3327869/1‑150 (MQ=255)
       cggcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATg                                                                                                                                                <  1:2307925/150‑1 (MQ=255)
       cggcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATg                                                                                                                                                <  2:4253922/150‑1 (MQ=255)
       cggcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATg                                                                                                                                                <  1:1482670/150‑1 (MQ=255)
       cggcgACAGCGACCAACGTTGCTCGGCGCGCGCACAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATg                                                                                                                                                <  1:1273921/150‑1 (MQ=255)
         gcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGa                                                                                                                                              <  2:4933531/150‑1 (MQ=255)
         gcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGa                                                                                                                                              <  1:717009/150‑1 (MQ=255)
          cgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGAc                                                                                                                                             >  2:2940488/1‑150 (MQ=255)
          cgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGCCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGag                                                                                                                                             <  2:540514/150‑2 (MQ=255)
           gACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGAcc                                                                                                                                            >  1:4590007/1‑150 (MQ=255)
            aCAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCt                                                                                                                                           <  2:4288180/150‑1 (MQ=255)
              aGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCg                                                                                                                                         >  1:3315494/1‑150 (MQ=255)
               gCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGa                                                                                                                                        <  1:1535754/150‑1 (MQ=255)
                 gACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                      <  2:599176/150‑1 (MQ=255)
                 gACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTACAGCTCTGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                      <  1:200572/150‑1 (MQ=255)
                  aCCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAggg                                                                                                                                     >  2:4826366/1‑150 (MQ=255)
                   ccAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGt                                                                                                                                    <  1:374025/150‑1 (MQ=255)
                   ccAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGt                                                                                                                                    <  1:1344290/150‑1 (MQ=255)
                         ttGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAAt                                                                                                                              <  2:1845923/150‑1 (MQ=255)
                         ttGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAAt                                                                                                                              <  2:1485278/150‑1 (MQ=255)
                           gCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGc                                                                                                                            <  1:3334853/150‑1 (MQ=255)
                               ggCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTc                                                                                                                        <  1:2688251/150‑1 (MQ=255)
                               ggCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTc                                                                                                                        <  2:1482900/150‑1 (MQ=255)
                                    gcgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAg                                                                                                                   >  1:5009572/1‑150 (MQ=255)
                                    gcgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAg                                                                                                                   >  1:1322904/1‑150 (MQ=255)
                                    gcgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAg                                                                                                                   >  1:2522962/1‑150 (MQ=255)
                                     cgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAgg                                                                                                                  >  1:2253968/1‑150 (MQ=255)
                                     cgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAgg                                                                                                                  <  2:4534021/150‑1 (MQ=255)
                                      gcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGc                                                                                                                 >  2:2290191/1‑150 (MQ=255)
                                        gcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATgg                                                                                                                                               <  1:880389/118‑1 (MQ=255)
                                           aaCCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCAc                                                                                                            <  2:3528645/150‑1 (MQ=255)
                                            aCCTCAGGCGCCAGGCCCTCCGGGCAGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGACAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACt                                                                                                           >  2:2573897/1‑150 (MQ=255)
                                             ccTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACtt                                                                                                          >  1:4614516/1‑150 (MQ=255)
                                                 aGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCt                                                                                                      <  2:3657602/150‑1 (MQ=255)
                                                   gcgcCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCa                                                                                                    <  1:5353544/150‑1 (MQ=255)
                                                     gcCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGc                                                                                                  >  1:5170029/1‑150 (MQ=255)
                                                     gcCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGc                                                                                                  <  2:752654/150‑1 (MQ=255)
                                                     gcCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGc                                                                                                  <  1:2200091/150‑1 (MQ=255)
                                                     gcCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGc                                                                                                  <  1:1167278/150‑1 (MQ=255)
                                                        aGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATg                                                                                               >  2:3418816/1‑150 (MQ=255)
                                                         ggCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATgg                                                                                              <  1:3257336/150‑1 (MQ=255)
                                                           cccTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCt                                                                                            >  2:5550396/1‑150 (MQ=255)
                                                             cTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCACGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTc                                                                                          <  1:2950269/150‑1 (MQ=255)
                                                              tCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCg                                                                                         >  2:4460358/1‑150 (MQ=255)
                                                                cGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGcc                                                                                       >  2:5055142/1‑150 (MQ=255)
                                                                   gccgccCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCcctc                                                                                    <  2:5037072/150‑1 (MQ=255)
                                                                    ccgccCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCcctcc                                                                                   <  2:5477130/150‑1 (MQ=255)
                                                                            caTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGAt                                                                           <  1:1820605/150‑1 (MQ=255)
                                                                            caTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGAt                                                                           <  1:3638486/150‑1 (MQ=255)
                                                                             aTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCGCCTGatccgtc                                                                          <  2:5103456/150‑7 (MQ=255)
                                                                              tCAAATTCACCGTGCGCGCATCCGGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCa                                                                         >  1:1796482/1‑150 (MQ=255)
                                                                              tCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCACTCCTGACGGATCa                                                                         >  2:4668524/1‑150 (MQ=255)
                                                                                  aTTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggcg                                                                     <  2:2371161/150‑1 (MQ=255)
                                                                                   ttCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggcgg                                                                    <  2:5158402/150‑1 (MQ=255)
                                                                                   ttCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggcgg                                                                    <  1:4143544/150‑1 (MQ=255)
                                                                                   ttCACCGTGCGCGCAGCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggcgg                                                                    <  1:3850300/150‑1 (MQ=255)
                                                                                        cGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGATCCGTCAGGAGGGCGaa                                                               <  2:141307/150‑2 (MQ=255)
                                                                                         gTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACg                                                              >  2:166207/1‑150 (MQ=255)
                                                                                          tGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGt                                                             >  1:1914069/1‑150 (MQ=255)
                                                                                           gcgcgcATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTc                                                            <  1:2940488/150‑1 (MQ=255)
                                                                                                cATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCAt                                                       >  2:3909023/1‑150 (MQ=255)
                                                                                                 aTCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCcctc                                                                                    <  2:2268970/120‑1 (MQ=255)
                                                                                                   ccAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCt                                                    <  2:2728099/150‑1 (MQ=255)
                                                                                                    cAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCtt                                                   <  1:4501575/150‑1 (MQ=255)
                                                                                                    cAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCGTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCtt                                                   <  2:3862962/150‑1 (MQ=255)
                                                                                                           cAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACgg                                            <  2:3496975/150‑1 (MQ=255)
                                                                                                           cAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACgg                                            <  1:4963228/150‑1 (MQ=255)
                                                                                                             gCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGt                                          >  1:4017249/1‑150 (MQ=255)
                                                                                                                gtactgCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGGCATTGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCTGGCGGGCTACGTCGGCATTCTTGCCACGGGTAtt                                       <  1:5269325/145‑1 (MQ=255)
                                                                                                                gTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTAtt                                       <  2:4637328/150‑1 (MQ=255)
                                                                                                                gTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTAtt                                       >  2:4970191/1‑150 (MQ=255)
                                                                                                                  ccAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCa                                     <  2:251604/150‑1 (MQ=255)
                                                                                                                  ccAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCa                                     <  1:742235/150‑1 (MQ=255)
                                                                                                                   cAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAg                                    >  1:1034605/1‑150 (MQ=255)
                                                                                                                     gCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCg                                  <  2:3104571/150‑1 (MQ=255)
                                                                                                                     gCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCg                                  <  1:310184/150‑1 (MQ=255)
                                                                                                                     gCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCg                                  <  1:323026/150‑1 (MQ=255)
                                                                                                                        cGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGa                               <  2:2693325/150‑1 (MQ=255)
                                                                                                                        cGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCCGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGa                               <  1:1589424/150‑1 (MQ=255)
                                                                                                                         ggTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGAt                              <  2:2058881/150‑1 (MQ=255)
                                                                                                                            cGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTc                           <  1:5483806/150‑1 (MQ=255)
                                                                                                                             gCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTcc                          <  1:5220442/150‑1 (MQ=255)
                                                                                                                             gCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTcc                          <  1:634124/150‑1 (MQ=255)
                                                                                                                               cAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTccac                        >  1:3822863/1‑150 (MQ=255)
                                                                                                                               cAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTccac                        <  1:4339911/150‑1 (MQ=255)
                                                                                                                               cAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTccac                        >  2:1105942/1‑150 (MQ=255)
                                                                                                                               cAGCTGCTCGGTACTCCAGGGCGTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTccac                        <  2:4062163/150‑1 (MQ=255)
                                                                                                                                 gctgctCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTccacca                      <  2:907789/150‑1 (MQ=255)
                                                                                                                                   tgctCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGc                    >  1:5947/1‑150 (MQ=255)
                                                                                                                                     ctCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTc                  >  2:2544109/1‑150 (MQ=255)
                                                                                                                                        ggTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAg               <  1:4219438/150‑1 (MQ=255)
                                                                                                                                          tACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCcgcg             >  2:4192254/1‑150 (MQ=255)
                                                                                                                                          tACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCg             >  1:3492045/1‑150 (MQ=255)
                                                                                                                                           ggtCCATGGCAAGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGa            >  2:4436965/3‑150 (MQ=255)
                                                                                                                                             tCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACa          <  1:1000629/150‑1 (MQ=255)
                                                                                                                                             tCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACa          >  2:3104091/1‑150 (MQ=255)
                                                                                                                                              caagggCCTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAg         >  2:2358930/5‑150 (MQ=255)
                                                                                                                                               cATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGc        >  1:5406202/1‑150 (MQ=255)
                                                                                                                                                    cTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTcggc   <  1:5470998/150‑1 (MQ=255)
                                                                                                                                                     ttGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTcggcg  <  1:1440989/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: