Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,472,205 T→A 7.8% K95* (AAG→TAG) ‡ 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,472,2050TA 7.8% 274.9 / 13.1 90K95* (AAG→TAG) ‡rlmH23S rRNA (pseudouridine(1915)‑N(3))‑methyltransferase RlmH
Reads supporting (aligned to +/- strand):  ref base T (42/41);  new base A (3/4);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

CAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTCGGC  >  NC_002947/5472059‑5472352
                                                                                                                                                  |                                                                                                                                                   
cAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGcc                                                                                                                                                  >  2:75533/1‑150 (MQ=255)
  gcggcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCAt                                                                                                                                                <  2:1457895/150‑1 (MQ=255)
     gcgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGa                                                                                                                                             <  2:866458/150‑1 (MQ=255)
      cgACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGAc                                                                                                                                            <  2:657599/150‑1 (MQ=255)
            cGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAg                                                                                                                                      <  2:42063/150‑1 (MQ=255)
             gACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                     <  2:992254/150‑1 (MQ=255)
             gACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                     <  2:35034/150‑1 (MQ=255)
             gACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                     <  1:569579/150‑1 (MQ=255)
                 aaCGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGa                                                                                                                                 >  2:207038/1‑150 (MQ=255)
                    gctGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACaa                                                                                                                              >  2:42870/3‑150 (MQ=255)
                     ttGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                     >  2:894111/1‑142 (MQ=255)
                     ttGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                     <  1:894111/142‑1 (MQ=255)
                     ttGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAAt                                                                                                                             >  2:1656436/1‑150 (MQ=255)
                     ttGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAAt                                                                                                                             <  2:44096/150‑1 (MQ=255)
                      tGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATg                                                                                                                            <  2:443132/150‑1 (MQ=255)
                         tCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGt                                                                                                                         <  1:207038/150‑1 (MQ=255)
                         tCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGt                                                                                                                         >  2:1177188/1‑150 (MQ=255)
                          cGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGtt                                                                                                                        <  1:1177188/150‑1 (MQ=255)
                           ggCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTc                                                                                                                       <  1:167212/150‑1 (MQ=255)
                               cgcgcgcAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCa                                                                                                                   <  1:1927008/150‑1 (MQ=255)
                                     cAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGc                                                                                                             <  1:52760/150‑1 (MQ=255)
                                      aaaCCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCACGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCa                                                                                                            >  2:2178784/1‑150 (MQ=255)
                                      aaaCCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGagta                                                                                                                                          >  1:1439302/1‑117 (MQ=255)
                                      aaaCCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGagta                                                                                                                                          <  2:1439302/120‑4 (MQ=255)
                                      aaaCCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGag                                                                                                                                            >  1:1640032/1‑117 (MQ=255)
                                      aaaCCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGag                                                                                                                                            <  2:1640032/118‑2 (MQ=255)
                                          cTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACttt                                                                                                        >  1:1829010/1‑150 (MQ=255)
                                          cTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGATCGGTACTACATGGATTGCCATGGACGTCGAGGGGGACAAGGCGTTCCCCAGGCTGCACttt                                                                                                        >  2:844364/1‑150 (MQ=255)
                                           tCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTg                                                                                                       >  2:1693842/1‑150 (MQ=255)
                                           tCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTg                                                                                                       <  1:1545378/150‑1 (MQ=255)
                                                     ggCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGACAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATgg                                                                                             >  1:296270/1‑150 (MQ=255)
                                                         cTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCa                                                                                                                   >  1:375452/1‑124 (MQ=255)
                                                         cTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCa                                                                                                                   <  2:375452/124‑1 (MQ=255)
                                                         cTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTc                                                                                         <  1:1693842/150‑1 (MQ=255)
                                                          tCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCg                                                                                        <  1:1128041/150‑1 (MQ=255)
                                                          tCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGAAAATGCGTTGCCCAGGCTGCACTTTGCTCAGCATGGCTTCg                                                                                        >  1:1101381/1‑150 (MQ=255)
                                                                   cccACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGa                                                                               >  1:2314696/1‑150 (MQ=255)
                                                                    ccaccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGCACTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGAc                                                                              >  2:437532/1‑150 (MQ=255)
                                                                    ccaccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGAc                                                                              <  2:404853/150‑1 (MQ=255)
                                                                     caccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAAATCGGGCTCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCTTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACg                                                                             >  1:1813394/1‑150 (MQ=255)
                                                                     caccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCAAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACg                                                                             >  1:695690/1‑150 (MQ=255)
                                                                       ccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCa                                                                                                                                                 >  1:1098275/1‑80 (MQ=255)
                                                                       ccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCa                                                                                                                                                 <  2:1098275/80‑1 (MQ=255)
                                                                       ccaTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGa                                                                           >  1:596902/1‑150 (MQ=255)
                                                                        caTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCAATCGCAACGCCGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGAt                                                                          >  1:1210761/1‑150 (MQ=255)
                                                                         aTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGag                                                                                                                                            >  2:2059441/1‑82 (MQ=255)
                                                                         aTCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCAAGCCATGGag                                                                                                                                            <  1:2059441/83‑2 (MQ=255)
                                                                          tCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCa                                                                        >  1:1091331/1‑150 (MQ=255)
                                                                            aaaTTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAgg                                                                      <  2:917224/150‑1 (MQ=255)
                                                                            aaaTTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAgg                                                                      <  2:2314696/150‑1 (MQ=255)
                                                                            aaaTTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAgg                                                                      >  2:1439518/1‑150 (MQ=255)
                                                                              aTTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggcg                                                                    >  1:1957183/1‑150 (MQ=255)
                                                                               ttCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCa                                                                                                                   <  1:1390592/102‑1 (MQ=255)
                                                                               ttCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCa                                                                                                                   >  2:1390592/1‑102 (MQ=255)
                                                                                  aCCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGc                                                                                                                           <  2:32715/91‑1 (MQ=255)
                                                                                  aCCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGc                                                                                                                           >  1:32715/1‑91 (MQ=255)
                                                                                  aCCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggcgggcg                                                                >  1:1957449/1‑150 (MQ=255)
                                                                                      tGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGt                                                            <  2:1957183/150‑1 (MQ=255)
                                                                                         gcgcATCCAGGCGCCAGCGGTCCAGCTCGGTCTCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTc                                                           <  1:1088597/148‑1 (MQ=255)
                                                                                         gcgcATCCAGGCGCCAGCGGTCCAGCTCGGTCTCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTc                                                           >  2:1088597/1‑148 (MQ=255)
                                                                                          cgcATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                     <  2:1498002/73‑1 (MQ=255)
                                                                                          cgcATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAgg                                                                                                                                     >  1:1498002/1‑73 (MQ=255)
                                                                                             aTCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCAtt                                                     >  1:102406/1‑150 (MQ=255)
                                                                                                cAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCACGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTCACGGATCAGGCGGGCGACGTCGGCATTCtt                                                  <  1:2178784/150‑1 (MQ=255)
                                                                                                  ggggcCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGc                                                <  2:387182/147‑1 (MQ=255)
                                                                                                       cAGCGTGCCCGCTCGGTCCCCAGCTGCTCGGTACTCCATGTCTTGCCATGGACCTCGAGTGTGACAATGCGTTCCCCAGCCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACgg                                           <  1:844364/150‑1 (MQ=255)
                                                                                                       cAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACgg                                           <  1:1584289/150‑1 (MQ=255)
                                                                                                       cAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACgg                                           <  2:102406/150‑1 (MQ=255)
                                                                                                        aGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACggg                                          >  2:1303934/1‑150 (MQ=255)
                                                                                                           ggTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTAt                                       <  2:1781633/150‑1 (MQ=255)
                                                                                                               cAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAg                                   >  1:1475793/1‑150 (MQ=255)
                                                                                                                 gcgcGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGAc                                                                              >  1:2330583/4‑105 (MQ=255)
                                                                                                                 gCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGAc                                                                              <  2:2330583/105‑1 (MQ=255)
                                                                                                                  cTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCgg                                <  2:2164868/150‑1 (MQ=255)
                                                                                                                   tCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCggg                               >  2:908538/1‑150 (MQ=255)
                                                                                                                       tCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGAttt                           <  2:1475793/150‑1 (MQ=255)
                                                                                                                         gCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCgg                                                         <  2:92353/118‑1 (MQ=255)
                                                                                                                         gCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCgg                                                         >  1:92353/1‑118 (MQ=255)
                                                                                                                              ctgctCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGGGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAg                    >  1:2045909/1‑150 (MQ=255)
                                                                                                                                    ggTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggcggg                                                                  >  1:2009311/1‑98 (MQ=255)
                                                                                                                                    ggTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAggcggg                                                                  <  2:2009311/98‑1 (MQ=255)
                                                                                                                                    ggTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAg              >  2:521797/1‑150 (MQ=255)
                                                                                                                                     gTACTCCCTGGCCTTTTATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGc             <  2:1813394/150‑1 (MQ=255)
                                                                                                                                      tACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTAtt                                      >  1:250662/1‑124 (MQ=255)
                                                                                                                                      tACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTAtt                                      <  2:250662/124‑1 (MQ=255)
                                                                                                                                         tCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAgatcgg         >  1:511800/1‑145 (MQ=255)
                                                                                                                                          ccATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAg        >  1:1938056/1‑150 (MQ=255)
                                                                                                                                           cATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGc       >  1:926423/1‑150 (MQ=255)
                                                                                                                                            aTGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCt      >  1:1867815/1‑150 (MQ=255)
                                                                                                                                                cTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTcggc  <  1:1656436/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAGCGGCGACAGCGACCAACGTTGCTCGGCGCGCGCGCAAACCTCAGGCGCCAGGCCCTCCGGGCCGCCCACCATCAAATTCACCGTGCGCGCATCCAGGCGCCAGCGGTCCAGCTCGGTCGCCAGCTGCTCGGTACTCCATGGCTTGCCATGGACCTCGAGGGTGACAATGCGTTCCCCAGGCTGCACTTTGCTCAGCATGGCTTCGCCCTCCTGACGGATCAGGCGGGCGACGTCGGCATTCTTGCCACGGGTATTCAGCGGGATTTCCACCAGCTCCAGCGACAGCTCGGC  >  NC_002947/5472059‑5472352

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