Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,850,353 G→A 6.2% V284V (GTG→GTA)  rlmD → 23S rRNA (uracil(1939)‑C(5))‑methyltransferase RlmD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,850,3530GA 6.2% 283.2 / 5.6 97V284V (GTG→GTA) rlmD23S rRNA (uracil(1939)‑C(5))‑methyltransferase RlmD
Reads supporting (aligned to +/- strand):  ref base G (47/44);  new base A (3/3);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GCCAGGTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAGGTGGTGGCAGT  >  NC_002947/1850204‑1850502
                                                                                                                                                     |                                                                                                                                                      
gccaggTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTg                                                                                                                                                        <  1:825157/150‑1 (MQ=255)
 ccaggTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGc                                                                                                                                                       <  1:2337060/150‑1 (MQ=255)
  caggTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCa                                                                                                                                                      <  1:1788564/150‑1 (MQ=255)
   aggTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAg                                                                                                                                                     <  1:2329097/150‑1 (MQ=255)
     gTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGt                                                                                                                                                   >  1:2346214/1‑150 (MQ=255)
      ttGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTg                                                                                                                                                  >  2:1104925/1‑150 (MQ=255)
        gCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGaa                                                                                                                                                >  2:1686628/1‑150 (MQ=255)
         cAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAAt                                                                                                                                               >  2:482677/1‑150 (MQ=255)
         cAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAAt                                                                                                                                               <  2:899069/150‑1 (MQ=255)
          aaGCGTTCTGCGCGCAGGCCAACGCCCAGCTGTGGCTGCAAGTCGACGGCGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTACAGGTGAATg                                                                                                                                              >  2:719753/1‑150 (MQ=255)
          aaGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATg                                                                                                                                              >  2:2315280/1‑150 (MQ=255)
          aaGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATg                                                                                                                                              >  2:2315500/1‑150 (MQ=255)
          aaGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATg                                                                                                                                              >  2:1143367/1‑150 (MQ=255)
          aaGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTACAGGTGAATg                                                                                                                                              >  2:717955/1‑150 (MQ=255)
          aaGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTACAGGTGAATg                                                                                                                                              >  2:717936/1‑150 (MQ=255)
            gCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGcc                                                                                                                                            <  2:972784/150‑1 (MQ=255)
               ttgtGCGAGCAGGGCCGCGCCCAGCTGTGGCGGCAGGGCGAAGGTGGTCCGGGGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCGGGCAGCTGGAGCTGGCCTGGCGCGCGGGCGATTTCGTGCAGGTGAATGCCCAg                                                                                                                                         <  2:980706/147‑1 (MQ=255)
               ttCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGccctg                                                                                                                                         <  2:1648589/150‑3 (MQ=255)
               ttCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAg                                                                                                                                         <  2:849517/150‑1 (MQ=255)
               ttCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGAATATATTTTTTCGTGCAGGTGAATGCCCAg                                                                                                                                         >  1:1181905/1‑150 (MQ=255)
                 cTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGt                                                                                                                                       <  1:726430/150‑1 (MQ=255)
                 cTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGCGCAGGTGACTGCACCGGt                                                                                                                                       <  2:693416/150‑1 (MQ=255)
                      aGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATa                                                                                                                                  <  1:426627/150‑1 (MQ=255)
                       gCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGACCTGGTCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATAc                                                                                                                                 >  1:1222834/1‑150 (MQ=255)
                          ggCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGc                                                                                                                              >  2:2156014/1‑150 (MQ=255)
                          ggCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGc                                                                                                                              <  2:2121367/150‑1 (MQ=255)
                             cAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCgat                                                                                                                           <  2:270527/150‑1 (MQ=255)
                             cAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCACTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCgat                                                                                                                           >  1:424219/1‑150 (MQ=255)
                                cGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCgatgat                                                                                                                        <  1:1417099/150‑1 (MQ=255)
                                   ccAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGa                                                                                                                     >  1:1291032/1‑150 (MQ=255)
                                    cAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAg                                                                                                                    >  1:420356/1‑150 (MQ=255)
                                    cAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAg                                                                                                                    >  1:1812746/1‑150 (MQ=255)
                                     aGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGc                                                                                                                   >  1:486123/1‑150 (MQ=255)
                                            gCTGCAAGGCGAAGGTGAACCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGCGCAGGCttt                                                                                                            >  2:27052/1‑150 (MQ=255)
                                             cTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTg                                                                                                           <  2:1848137/150‑1 (MQ=255)
                                             cTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTCCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTg                                                                                                           <  2:1848928/150‑1 (MQ=255)
                                              tGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTgg                                                                                                          >  2:898845/1‑150 (MQ=255)
                                               gCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGc                                                                                                         <  1:239732/150‑1 (MQ=255)
                                               gCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGc                                                                                                         <  1:898845/150‑1 (MQ=255)
                                                cAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGcc                                                                                                        >  1:1565193/1‑150 (MQ=255)
                                                    gCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgc                                                                                                    <  1:2315500/150‑1 (MQ=255)
                                                    gCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgc                                                                                                    <  1:2315280/150‑1 (MQ=255)
                                                    gCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgc                                                                                                    >  1:947105/1‑150 (MQ=255)
                                                      gAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCt                                                                                                       >  2:2056487/1‑145 (MQ=255)
                                                      gAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCt                                                                                                       <  1:2056487/145‑1 (MQ=255)
                                                      gAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgctc                                                                                                  >  2:936994/1‑150 (MQ=255)
                                                       aaGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgctcg                                                                                                 <  1:406177/150‑1 (MQ=255)
                                                          gTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGgctcgctc                                                                                              >  2:782618/1‑150 (MQ=255)
                                                           tGAGCCGGCGCCCGTGGACCCAGCGGGGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTcc                                                                                             >  1:322479/1‑150 (MQ=255)
                                                              gCCGGCGCCCGTGGACCCAGCGGCGCAGTTGAGATTTGCCCTGGCCCCCTGACAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCcaca                                                                                          >  1:1566718/1‑150 (MQ=255)
                                                                cGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAg                                                                                        >  1:1440843/1‑150 (MQ=255)
                                                                    gcCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCcga                                                                                    <  1:2278816/150‑1 (MQ=255)
                                                                       cGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCcgncga                                                                                 <  2:1497175/150‑1 (MQ=255)
                                                                       cGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCcgacga                                                                                 >  2:1746621/1‑150 (MQ=255)
                                                                       cGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCcgacga                                                                                 <  1:1686628/150‑1 (MQ=255)
                                                                        gTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAg                                                                                >  2:2135628/1‑150 (MQ=255)
                                                                         tGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGc                                                                               <  2:2346214/150‑1 (MQ=255)
                                                                              ccAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAgg                                                                                                                                        >  2:1765876/1‑88 (MQ=255)
                                                                              ccAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAgg                                                                                                                                        <  1:1765876/88‑1 (MQ=255)
                                                                              ccAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGtgtg                                                                          <  1:1104925/150‑1 (MQ=255)
                                                                              ccAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGtgtg                                                                          >  2:744660/1‑150 (MQ=255)
                                                                              ccAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGtgtg                                                                          >  1:180869/1‑150 (MQ=255)
                                                                                   ggcgCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTACAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGa                                                                     <  1:717955/150‑1 (MQ=255)
                                                                                   ggcgCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTACAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGa                                                                     <  1:717936/150‑1 (MQ=255)
                                                                                   ggcgCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTACAGGGGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGa                                                                     <  1:719753/150‑1 (MQ=255)
                                                                                      gcAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACtt                                                                  >  2:1587499/1‑150 (MQ=255)
                                                                                            gggCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGc                                                                               <  2:1136641/131‑1 (MQ=255)
                                                                                            gggCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGc                                                                               >  1:1136641/1‑131 (MQ=255)
                                                                                              gCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCg                                                          <  2:582784/150‑1 (MQ=255)
                                                                                                tttGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACtt                                                                  <  2:752744/140‑1 (MQ=255)
                                                                                                tttGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACtt                                                                  >  1:752744/1‑140 (MQ=255)
                                                                                                    cccTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTgg                                                    <  1:482677/150‑1 (MQ=255)
                                                                                                    cccTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTgg                                                    >  2:1284307/1‑150 (MQ=255)
                                                                                                    cccTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTgg                                                    <  2:1384382/150‑1 (MQ=255)
                                                                                                    cccTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATGTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTTGCCTGGCTCGCTCCACCAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTgg                                                    >  2:824494/1‑150 (MQ=255)
                                                                                                     ccTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTggg                                                   <  2:1565193/150‑1 (MQ=255)
                                                                                                             ccTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTg                                           >  1:1756379/1‑150 (MQ=255)
                                                                                                                 gCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑t                                      >  2:1767440/1‑150 (MQ=255)
                                                                                                                   agctggaactggCCTGACGCCCGGACGATTTCGTGCAGGAGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGACGACGAGCGTGTGATGGACTTGTTCTGCTGTCTGGGCAACTTTGCCC‑TGc                                    >  2:1011934/1‑150 (MQ=255)
                                                                                                                        gagctggCCTGGCGCCCGGGCGATTTCGTGCAGGAGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGATTTGGCCTGGCTCGCGCCACAAGCCGACGAGGGGGTGCTGGCCTTGTACTGCGGTCTGGGCAACTTTGCCT‑TGCCGcag                               >  2:832436/1‑148 (MQ=255)
                                                                                                                            tggcctggcGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGccc                           <  1:178795/150‑1 (MQ=255)
                                                                                                                              gcctggcGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGc                         >  2:1794662/1‑150 (MQ=255)
                                                                                                                                ctggcGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCa                       >  2:1062746/1‑150 (MQ=255)
                                                                                                                                   gcgcCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGgc                    <  2:812835/150‑1 (MQ=255)
                                                                                                                                     gcCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGgc                    >  2:377400/1‑148 (MQ=255)
                                                                                                                                     gcCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGgc                    <  1:377400/148‑1 (MQ=255)
                                                                                                                                        cGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTg               <  2:1222834/150‑1 (MQ=255)
                                                                                                                                         gggCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGa              >  1:1602701/1‑150 (MQ=255)
                                                                                                                                           gCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCCTTGCCGCTGGCCCGCCAGGCGCGTGAg             <  2:2112841/150‑1 (MQ=255)
                                                                                                                                            cGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAggt           <  2:947105/150‑1 (MQ=255)
                                                                                                                                             gATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAggtg          >  1:2086471/1‑150 (MQ=255)
                                                                                                                                              aTTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAggtgg         <  1:1794662/150‑1 (MQ=255)
                                                                                                                                               tttCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTggtc                                                                    >  2:1666536/1‑89 (MQ=255)
                                                                                                                                               tttCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTggtc                                                                    <  1:1666536/91‑3 (MQ=255)
                                                                                                                                                 tCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAGGtggtgg      <  2:1743301/150‑1 (MQ=255)
                                                                                                                                                   gTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAGGTGGTGGCa    >  1:1769248/1‑150 (MQ=255)
                                                                                                                                                    tGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAGGTGGTGGCAg   <  1:1767440/150‑1 (MQ=255)
                                                                                                                                                    tGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGCGAGGTGGTGGCAg   >  2:230597/1‑150 (MQ=255)
                                                                                                                                                     gCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAGGTGGTGGCAGt  <  2:418831/150‑1 (MQ=255)
                                                                                                                                                     gCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAGGTGGTGGCAGt  <  1:1587499/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCCAGGTTGCAAGCGTTCTGCGAGCAGGCCAACGCCCAGCTGTGGCTGCAAGGCGAAGGTGAGCCGGCGCCCGTGGACCCAGCGGCGCAGTTGGGCTTTGCCCTGGCCCCCTGGCAGCTGGAGCTGGCCTGGCGCCCGGGCGATTTCGTGCAGGTGAATGCCCAGGTCAATACGGCGATGATCGAGCAGGCTTTGGCCTGGCTCGCTCCACAAGCCGACGAGCGTGTGCTGGACTTGTTCTGCGGTCTGGGCAACTTTGCCC‑TGCCGCTGGCCCGCCAGGCGCGTGAGGTGGTGGCAGT  >  NC_002947/1850204‑1850502

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