Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,685,956 G→C 6.8% P331P (CCC→CCG) ‡ PP_4988 ← signal transduction histidine kinase CheA/CheY‑like receiver

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,685,9560GC 6.8% 416.8 / 17.9 133P331P (CCC→CCG) ‡PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base G (58/66);  new base C (5/4);  total (63/70)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 139 total aligned reads displayed.

AAGTGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCT  >  NC_002947/5685807‑5686104
                                                                                                                                                     |                                                                                                                                                     
aaGTGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAg                                                                                                                                                       >  2:3063576/1‑150 (MQ=255)
  gTGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAggg                                                                                                                                                     <  2:2704564/150‑1 (MQ=255)
       tCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCa                                                                                                                                                <  2:445528/150‑1 (MQ=255)
        cAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAg                                                                                                                                               >  2:818694/1‑150 (MQ=255)
          gCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACATCTCGGCCACGGCCGGCAGCATCAACATGGCCAATTCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGc                                                                                                                                             <  2:1722252/150‑1 (MQ=255)
          gCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCGTTAACATGGGCCATGCCCCGCTTATCGGCGGCAACAGCCCCGGGAAGGGCTggcggg                                                                                                                                             >  2:2337905/1‑145 (MQ=255)
               tcaATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAgcccttgc                                                                                                                                        <  1:1255037/150‑8 (MQ=255)
                    cggcggCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgggg                                                                                                                                   <  2:3053774/150‑1 (MQ=255)
                    cggcggCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgggg                                                                                                                                   <  1:2611087/150‑1 (MQ=255)
                    cggcggCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgggg                                                                                                                                   <  1:67400/150‑1 (MQ=255)
                     ggcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGc                                                                                                                                  <  1:818694/150‑1 (MQ=255)
                     ggcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGc                                                                                                                                        >  1:2793784/1‑144 (MQ=255)
                     ggcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGc                                                                                                                                        <  2:2793784/144‑1 (MQ=255)
                          ctctTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTc                                                                                                                             <  1:2948662/150‑1 (MQ=255)
                             ttCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGa                                                                                                                          <  1:988081/150‑1 (MQ=255)
                              tcctccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGAc                                                                                                                        >  2:2225309/1‑150 (MQ=255)
                                ctccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTc                                                                                                                      <  1:1510845/150‑1 (MQ=255)
                                   cAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                   <  1:331133/150‑1 (MQ=255)
                                   cAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                   >  1:1678024/1‑150 (MQ=255)
                                   cAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATCGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGCCTCCAg                                                                                                                   >  1:2974226/1‑150 (MQ=255)
                                    aGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAgg                                                                                                                  >  1:2198331/1‑150 (MQ=255)
                                         ctgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGc                                                                                                                                             >  1:2162820/1‑119 (MQ=255)
                                         ctgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATc                                                                                                             >  1:2169450/1‑150 (MQ=255)
                                         ctgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATc                                                                                                             <  2:2370910/150‑1 (MQ=255)
                                          tgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCGGGGct                                                                                                                                 <  1:616767/130‑2 (MQ=39)
                                          tgcAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCGGGGct                                                                                                                                 >  2:616767/1‑129 (MQ=39)
                                              accaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGa                                                                                                        <  2:340775/150‑1 (MQ=255)
                                               ccaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGAt                                                                                                       >  1:2531572/1‑150 (MQ=255)
                                                caacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAg                                                                                                                                                       >  2:1179201/1‑102 (MQ=255)
                                                caacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATc                                                                                                      >  1:2403848/1‑150 (MQ=255)
                                                 aacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCa                                                                                                     <  1:1383391/150‑1 (MQ=255)
                                                      tGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGc                                                                                                                                                    >  1:202496/1‑99 (MQ=255)
                                                      tGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGc                                                                                                                                                    <  2:202496/99‑1 (MQ=255)
                                                      tGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCGGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATcc                                                                                                <  1:2382834/150‑1 (MQ=255)
                                                         acaTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                   >  1:1238672/1‑128 (MQ=255)
                                                               cTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCctgc                                                                                       <  1:165455/150‑1 (MQ=255)
                                                                 tGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCctgctg                                                                                     <  2:934053/150‑1 (MQ=255)
                                                                      aGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAaca                                                                                <  2:2198331/150‑1 (MQ=255)
                                                                          gACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGtt                                                                            <  2:560069/150‑1 (MQ=255)
                                                                           aCAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTc                                                                           >  2:2248301/1‑150 (MQ=255)
                                                                           aCAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTc                                                                           <  1:1154908/150‑1 (MQ=255)
                                                                            cAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCa                                                                          >  2:414792/1‑150 (MQ=255)
                                                                             aaaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAg                                                                         >  2:322302/1‑150 (MQ=255)
                                                                             aaaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAg                                                                         >  1:1128342/1‑150 (MQ=255)
                                                                             aaaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAg                                                                         >  1:509932/1‑150 (MQ=255)
                                                                              aaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                        >  1:2812304/1‑150 (MQ=255)
                                                                               aCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                   >  1:2563363/1‑106 (MQ=255)
                                                                               aCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                   <  2:2563363/106‑1 (MQ=255)
                                                                                 aGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgggg                                                                                                                                   >  1:1150347/1‑89 (MQ=255)
                                                                                 aGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCgggg                                                                                                                                   <  2:1150347/89‑1 (MQ=255)
                                                                                 aGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                        <  2:2842721/147‑1 (MQ=255)
                                                                                 aGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                        >  1:2842721/1‑147 (MQ=255)
                                                                                 aGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGat                                                                     >  1:552384/1‑150 (MQ=255)
                                                                                    tCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGattca                                                                  <  2:1128342/150‑1 (MQ=255)
                                                                                     cGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGattcat                                                                 <  2:216687/150‑1 (MQ=255)
                                                                                       gCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCa                                                                          <  2:2814899/139‑1 (MQ=255)
                                                                                       gCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGattcattc                                                               >  2:2866408/1‑150 (MQ=255)
                                                                                        ccACGGCCGGCAGCATCAACATGGCCAATGCAAAGCTATCCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGAttcattca                                                              >  1:1491964/1‑150 (MQ=255)
                                                                                         gaccgccggcAGCATCAACGTGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAg                                                             <  2:493293/146‑1 (MQ=255)
                                                                                          aCGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAggg                                                                                                                                                     <  2:878834/62‑1 (MQ=255)
                                                                                          aCGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAggg                                                                                                                                                     >  1:878834/1‑62 (MQ=255)
                                                                                              ccggcAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCa                                                        <  1:1087520/150‑1 (MQ=255)
                                                                                              ccgccAGCCTCAACAGGGCCACTGCCCCGCGGATCTGCTCCAACAGCCCCGGCAAGGGCGGCAGGCGTTGCCGGGGGCAGTCGGACTCCAGGCAAGCGCCGTGCAGATCCTTGGGCTGCTGCAACACGTTCAGAGATTCCCTCCGCACCa                                                        <  2:2974226/150‑1 (MQ=255)
                                                                                                    gcaTCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGc                                                  <  2:2169450/150‑1 (MQ=255)
                                                                                                     caTCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCg                                                 <  1:2866408/150‑1 (MQ=255)
                                                                                                           cATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGc                                                                                                                                  <  1:768107/64‑1 (MQ=255)
                                                                                                           cATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGc                                                                                                                                  >  2:768107/1‑64 (MQ=255)
                                                                                                                ccAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCATGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCa                                      >  1:857990/1‑150 (MQ=255)
                                                                                                                ccAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCa                                      >  2:2084856/1‑150 (MQ=255)
                                                                                                                  aaTGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGAt                                                                                                       <  2:1963622/83‑1 (MQ=255)
                                                                                                                     gCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTcc                                 <  2:2531572/150‑1 (MQ=255)
                                                                                                                      ccccGCGCATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGCCACCGGGCTGGGGGCAGGTGGAATGTAGCTGCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCg                                <  1:361851/150‑1 (MQ=255)
                                                                                                                       cccGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGTCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGGTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGa                               >  1:1616492/1‑150 (MQ=255)
                                                                                                                       cccGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCctgctg                                                                                     >  2:3025667/1‑96 (MQ=255)
                                                                                                                        ccGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAg                              >  1:1298278/1‑150 (MQ=255)
                                                                                                                         cgcgGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTc                                  >  1:447403/1‑145 (MQ=255)
                                                                                                                            ggATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCg                          <  2:1298278/150‑1 (MQ=255)
                                                                                                                               tCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCTGGGCCAGTCGTACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACCTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCa                       <  1:1174824/150‑1 (MQ=255)
                                                                                                                                  gctgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCg                          <  2:1261467/144‑1 (MQ=255)
                                                                                                                                  gctgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGc                    >  2:1183797/1‑150 (MQ=255)
                                                                                                                                    tgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTg                  >  1:679771/1‑150 (MQ=255)
                                                                                                                                     gcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGc                 <  1:322302/150‑1 (MQ=255)
                                                                                                                                       aaCAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCt               <  2:2812304/150‑1 (MQ=255)
                                                                                                                                       aaCAGCCACGGCAAGGCCTACAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCt               >  2:2492334/1‑150 (MQ=255)
                                                                                                                                        ctaGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCtt              <  1:2038558/148‑1 (MQ=255)
                                                                                                                                         cAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGattcatt                                                                <  1:395534/99‑1 (MQ=255)
                                                                                                                                         cAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCCTg             <  1:3063576/150‑1 (MQ=255)
                                                                                                                                           gCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGc           <  2:1678746/150‑1 (MQ=255)
                                                                                                                                           gCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTAATTCAGCAACAGCTCGCGTATCTCGTCGATGTCCGAGCCCGGCAGGCTGCCTTgag           >  2:2038558/1‑148 (MQ=255)
                                                                                                                                            ccccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGcc                                                                                          <  2:2901604/70‑1 (MQ=255)
                                                                                                                                            ccccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGACGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGcc          >  1:1745832/1‑150 (MQ=255)
                                                                                                                                            ccccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGAACCAGCTGGCGTATCTCGGCCAGGGCCGAGCCCGGCAAGCCGcagtgccc          >  1:2867223/1‑143 (MQ=255)
                                                                                                                                             cccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCg         >  1:1806926/1‑150 (MQ=255)
                                                                                                                                             cccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCg         >  1:1380835/1‑150 (MQ=255)
                                                                                                                                             cccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACACCAGGTAATCGCCGATCAGAACCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCAGGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCg         >  1:1013333/1‑150 (MQ=255)
                                                                                                                                              ccGGCGACGCCTGCAGCCCTTGCCGGGGCTAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGt        >  2:767345/1‑150 (MQ=255)
                                                                                                                                              ccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGTCCTGCTCCAACACGTTCAGCTATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGt        <  2:1185442/150‑1 (MQ=255)
                                                                                                                                              ccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTAGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGt        >  1:2320788/1‑150 (MQ=255)
                                                                                                                                                  cAAGTGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGAGCAGATCCTTGGCCTGTTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTtc    <  2:592536/150‑1 (MQ=255)
                                                                                                                                                    aGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTtcct  <  2:509932/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGTGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCT  >  NC_002947/5685807‑5686104

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: