Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,546,502 T→C 13.5% V10V (GTA→GTG rlmB ← 23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,546,5020TC13.5% 248.4 / 22.0 101V10V (GTA→GTGrlmB23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base T (36/51);  new base C (5/9);  total (41/60)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

TCTCACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACTGCTGCG  >  NC_002947/5546354‑5546644
                                                                                                                                                    |                                                                                                                                                
tctcACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTa                                                                                                                                                 <  2:1576541/150‑1 (MQ=255)
 ctcACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTAc                                                                                                                                                <  2:1329287/150‑1 (MQ=255)
  tcACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACg                                                                                                                                               >  1:832003/1‑150 (MQ=255)
         gcTCGGCCTGGCCGACCTGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgca                                                                                                                                        <  2:158168/150‑4 (MQ=255)
         gcTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgca                                                                                                                                        <  2:902066/150‑4 (MQ=255)
         gcTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgca                                                                                                                                        <  1:2316630/150‑4 (MQ=255)
          cTCGGCCTGGCCGGCCGGCACGCGGTTTTGAGCGCCATGAGCCGGAACAGTGTCGGGGCGTGCCGCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGAAGCAATGCCTGCACGTCTTACACGCCGtgcag                                                                                                                                       <  2:2366301/150‑5 (MQ=255)
          cTCGGCCTGGCCGACCGGCACGCGCTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgcag                                                                                                                                       <  2:1212083/150‑5 (MQ=255)
           tCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGGTTCGAGGCGTGGCGCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCGGCAGCAATGCCTGCACGGCGTGTACGCCGTAGGtt                                                                                                                                     >  2:853509/1‑150 (MQ=255)
               ccTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGAtttttt                                                                                                                                 <  2:939663/150‑1 (MQ=255)
               ccTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGAtttttt                                                                                                                                 >  1:939663/1‑150 (MQ=255)
               ccTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgcag                                                                                                                                       <  1:1513762/145‑5 (MQ=255)
               ccTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgcag                                                                                                                                       >  2:1513762/1‑141 (MQ=255)
                cTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTc                                                                                                                                <  2:569264/150‑1 (MQ=255)
                  ggCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCa                                                                                                                              <  2:2171691/150‑1 (MQ=255)
                   gCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAg                                                                                                                             <  1:1783841/150‑1 (MQ=255)
                    ccGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGc                                                                                                                            <  1:1859784/150‑1 (MQ=255)
                     cGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCt                                                                                                                           <  1:1822588/150‑1 (MQ=255)
                     cGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCt                                                                                                                           >  2:765147/1‑150 (MQ=255)
                       aCCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGa                                                                                                                         <  1:1673404/150‑1 (MQ=255)
                       aCCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGa                                                                                                                         <  1:795085/150‑1 (MQ=255)
                         cGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACt                                                                                                                       <  1:1521833/150‑1 (MQ=255)
                            cACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCAt                                                                                                                    >  2:427975/1‑150 (MQ=255)
                                cgGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAAt                                                                                                                <  1:1032043/150‑1 (MQ=255)
                                 ggTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCg                                                                                                                                            >  2:2099569/1‑122 (MQ=255)
                                 ggTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCg                                                                                                                                            <  1:2099569/122‑1 (MQ=255)
                                 ggTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAAtt                                                                                                               <  1:926463/150‑1 (MQ=255)
                                     tcccGCGGCCAGCGTCAGCAGGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTGGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGctc                                                                                                           <  1:853509/148‑1 (MQ=255)
                                        cgcgGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTg                                                                                                        >  1:2096621/1‑150 (MQ=255)
                                         gcgGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGt                                                                                                      <  1:1776277/150‑1 (MQ=255)
                                                gcgccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCg                                                                                               >  1:1448341/1‑150 (MQ=255)
                                                 cgccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGttttttgcg                                                                                              >  2:166436/1‑146 (MQ=255)
                                                  gccagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTg                                                                                             <  1:1214193/150‑1 (MQ=255)
                                                    cagcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGa                                                                                           <  1:561928/150‑1 (MQ=255)
                                                     agcAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGCGGAg                                                                                          >  2:1181521/1‑150 (MQ=255)
                                                       cAAGGTTTCGAGGGGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAgcg                                                                                        >  2:406155/1‑150 (MQ=255)
                                                       cAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAgcg                                                                                        >  2:1236965/1‑150 (MQ=255)
                                                            tttCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGa                                                                                   <  1:1270678/150‑1 (MQ=255)
                                                             ttCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGAc                                                                                  <  2:832003/150‑1 (MQ=255)
                                                             ttCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGAc                                                                                  >  2:2388182/1‑150 (MQ=255)
                                                                   gCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCCATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCCCTTGTTTTTTCGTGGAGCGCCGGACTGcggc                                                                            >  2:419356/1‑150 (MQ=255)
                                                                     gTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcggcg                                                                           >  1:1694271/1‑149 (MQ=255)
                                                                     gTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcggcg                                                                           <  2:1694271/149‑1 (MQ=255)
                                                                     gTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGCACGCCGTAGATTTTTTCTATACGTTTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcggcgg                                                                          <  1:1241375/150‑1 (MQ=255)
                                                                     gTGGCTCGCTGCGCCCTTCCGAAAGAAAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcggcgg                                                                          >  2:351061/1‑150 (MQ=255)
                                                                         cTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGa                                                                                           >  1:79625/1‑129 (MQ=255)
                                                                         cTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGa                                                                                           >  1:79240/1‑129 (MQ=255)
                                                                         cTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGa                                                                                           <  2:79240/129‑1 (MQ=255)
                                                                         cTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGa                                                                                           <  2:79625/129‑1 (MQ=255)
                                                                          tCGCTGCGCCCTTCCGACAGCCAGATGTTGTGGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCtt                                                                     >  1:2276818/1‑150 (MQ=255)
                                                                          tCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCtt                                                                     <  2:335820/150‑1 (MQ=255)
                                                                                gggccTTACGACGGCCAGTTCTGTTTGACCCGCTTTGGAGGTTGCTGGAGCAGTGCCTGCACGGTGTGTACGCCGTCGATTGCTTCCAGCTGACTCGTGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGTGCCTTTCCGGt                                                               <  2:87405/147‑1 (MQ=255)
                                                                                    cTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGGTGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCtt                                                           <  1:605195/150‑1 (MQ=255)
                                                                                     ttCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGAtt                                                                                                                                     <  1:2410057/76‑1 (MQ=255)
                                                                                     ttCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGAtt                                                                                                                                     >  2:2410057/1‑76 (MQ=255)
                                                                                      tCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTgg                                                         <  2:452101/150‑1 (MQ=255)
                                                                                        cGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTc                                                       <  1:1940681/150‑1 (MQ=255)
                                                                                        cGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGCGCTTGGTc                                                       <  2:1908283/150‑1 (MQ=255)
                                                                                           cAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGt                                                    >  1:2117355/1‑150 (MQ=255)
                                                                                             gCCAGATCTGCTTTACCCGCCTTGGATGGTGCTGCAGCAAAGCCTGCACGGCGTGTACCCCTTACATCTATGTCGAACATCTCATGGAGTGCTCTTGTTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGtt                                                   <  1:166436/150‑1 (MQ=255)
                                                                                             gCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGttt                                                  <  1:765147/150‑1 (MQ=255)
                                                                                              ccAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcggcgg                                                                          <  2:580916/125‑1 (MQ=255)
                                                                                              ccAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTcggcgg                                                                          >  1:580916/1‑125 (MQ=255)
                                                                                              ccAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGtt                                                   <  2:1648405/148‑1 (MQ=255)
                                                                                              ccAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGtt                                                   >  1:1648405/1‑148 (MQ=255)
                                                                                               cAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCTTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGc                                                >  2:1365185/1‑150 (MQ=255)
                                                                                                aGATCTGCTTGTCCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGcc                                               <  1:1642442/150‑1 (MQ=255)
                                                                                                aGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgc                                                                                                                                         >  2:2163971/1‑60 (MQ=38)
                                                                                                aGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTACACGCCGtgc                                                                                                                                         <  1:2163971/62‑3 (MQ=38)
                                                                                                     tGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGt                                                               <  1:2034904/129‑1 (MQ=255)
                                                                                                     tGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGt                                                               >  2:2034904/1‑129 (MQ=255)
                                                                                                     tGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCg                                              <  2:2268834/146‑1 (MQ=255)
                                                                                                     tGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCg                                              >  1:2268834/1‑146 (MQ=255)
                                                                                                      gCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGc                                                                                                                            >  2:1752441/1‑68 (MQ=255)
                                                                                                      gCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGc                                                                                                                            <  1:1752441/68‑1 (MQ=255)
                                                                                                       cTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGc                                        >  1:185923/1‑150 (MQ=255)
                                                                                                         tGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCtt                                      <  2:2276818/150‑1 (MQ=255)
                                                                                                              cGCTTTGGATGGTGCTGCAGCAATGCCTGCAGGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATGGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCAc                                 <  2:1257440/150‑1 (MQ=255)
                                                                                                              cGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCAc                                 <  2:1758329/150‑1 (MQ=255)
                                                                                                               gCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCAcc                                >  2:2304951/1‑150 (MQ=255)
                                                                                                                     gATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTTTTTTTTTTTCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTtctc                          >  2:480206/1‑150 (MQ=255)
                                                                                                                     gATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATCTTTTCCTGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTtctc                          >  2:1532753/1‑150 (MQ=255)
                                                                                                                        ggTGTTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGATGGCTTGCCACCTTTCTCCGa                       <  1:351061/150‑1 (MQ=255)
                                                                                                                         gTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGAt                      <  1:427975/150‑1 (MQ=255)
                                                                                                                          tgctgcAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTCCGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGAtt                     <  2:1663585/150‑1 (MQ=255)
                                                                                                                              gcagcaATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTtgct                 >  1:2249849/1‑150 (MQ=255)
                                                                                                                              gcagcaATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGAAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTtgct                 <  1:1236965/150‑1 (MQ=255)
                                                                                                                               cagcaATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTtgctg                <  1:2082818/150‑1 (MQ=255)
                                                                                                                                     tGCCTGCACGGAGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTAGATCTCCAGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCa          <  1:2388182/150‑1 (MQ=255)
                                                                                                                                       ccTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGAc                                           >  1:1134118/1‑115 (MQ=255)
                                                                                                                                       ccTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATGGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGAc                                           <  2:1134118/115‑1 (MQ=255)
                                                                                                                                        cTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCActg       <  1:419356/150‑1 (MQ=255)
                                                                                                                                         tGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCActgc      >  2:1766972/1‑150 (MQ=255)
                                                                                                                                           tatggCGTACACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCtt                                      <  2:189196/113‑1 (MQ=255)
                                                                                                                                           taCGGCGTACACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGTTTGTCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCActgctg    >  2:2275859/2‑150 (MQ=255)
                                                                                                                                           taCGGCGTACACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCtt                                      >  1:189196/2‑116 (MQ=255)
                                                                                                                                            aCGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACtgctgc   <  2:262827/150‑1 (MQ=255)
                                                                                                                                             cGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACTGCTgcg  >  2:2359731/1‑150 (MQ=255)
                                                                                                                                                 gtgtACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGttt                                                  <  2:487796/98‑1 (MQ=255)
                                                                                                                                                 gtgtACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGttt                                                  >  1:487796/1‑98 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCTCACGGCGCTCGGCCTGGCCGACCGGCACGCGGTTTTCCGCGGCCAGCGCCAGCAAGGTTTCGAGGCGTGGCTCGCTGCGCCCTTCCGACAGCCAGATCTGCTTGACCCGCTTTGGATGGTGCTGCAGCAATGCCTGCACGGCGTGTACGCCGTAGATTTTTTCCAGCTGACTCATGAATTGCTCTTGTTTTTGCGTGGAGCGCCGGACTTCGGCGGGCCTTTACGGTGCTTGGTCGGTTTGCCGGACGGCTTGCCACCTTTCTCCGATTTGCTGCCAGCACTGCTGCG  >  NC_002947/5546354‑5546644

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