Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,685,902 G→C 12.3% P331A (CCC→GCC) ‡ 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,9020GC12.3% 198.5 / 15.4 72P331A (CCC→GCC) ‡PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base G (25/38);  new base C (7/2);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 6.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

TGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTC  >  NC_002947/5685754‑5686049
                                                                                                                                                    |                                                                                                                                                   
tGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCTACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAggg                                                                                                                                                     <  1:86179/149‑1 (MQ=255)
tGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAggg                                                                                                                                                     <  8:250557/149‑1 (MQ=255)
      aGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAg                                                                                                                                               >  2:62955/1‑149 (MQ=255)
      aGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAg                                                                                                                                               >  5:279727/1‑149 (MQ=255)
        cTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGc                                                                                                                                             <  8:276598/149‑1 (MQ=255)
          catcaATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGt                                                                                                                                           >  3:233950/1‑149 (MQ=255)
           atcaATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGtt                                                                                                                                          <  5:110620/149‑1 (MQ=255)
              aaTCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGcc                                                                                                                                       <  7:91235/149‑1 (MQ=255)
               aacggcggATCTTACTCCACGTGCTGAAACCAACGCTGCACATAGCTTTCGCAGCAGACAAACAGCTCGTACACGGCCGGCAGCATCAACATAGCCAATGCAACGCTCATCTGCTGCAACATCAACTTCAACGCCTTCAGGCGtttacg                                                                                                                                      >  2:171164/3‑145 (MQ=255)
                  ggcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCggg                                                                                                                                    <  6:279727/148‑1 (MQ=255)
                  ggcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCggg                                                                                                                                    >  3:295459/1‑148 (MQ=255)
                  ggcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCggg                                                                                                                                    <  4:295459/148‑1 (MQ=255)
                   gcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGc                                                                                                                                  >  2:301926/1‑149 (MQ=255)
                   gcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGc                                                                                                                                  <  5:246879/149‑1 (MQ=255)
                   gcgacTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGc                                                                                                                                  <  3:335459/149‑1 (MQ=255)
                     ggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCa                                                                                                                                <  5:167323/149‑1 (MQ=255)
                          ttCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCg                                                                                                                            >  8:319208/1‑148 (MQ=255)
                            cctccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCg                                                                                                                            <  6:148323/146‑1 (MQ=255)
                              tccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTc                                                                                                                       >  4:205449/1‑149 (MQ=255)
                              tccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCAGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTc                                                                                                                       >  4:56023/1‑149 (MQ=255)
                                cAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAACATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTc                                                                                                                       >  6:315920/1‑147 (MQ=255)
                                 aGGGGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCACGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg                                                                                                                    <  2:224227/149‑1 (MQ=255)
                                   gTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTc                                                                                                                       >  7:264080/1‑144 (MQ=255)
                                   gTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTc                                                                                                                       <  8:264080/144‑1 (MQ=255)
                                         caaccaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGc                                                                                                            <  5:247041/149‑1 (MQ=255)
                                         caaccaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGc                                                                                                            >  1:51106/1‑149 (MQ=255)
                                         caaccaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGc                                                                                                            >  1:114114/1‑149 (MQ=255)
                                          aaccaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGcc                                                                                                           >  2:382955/1‑149 (MQ=255)
                                            ccaacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGa                                                                                                         <  3:205449/149‑1 (MQ=255)
                                             caacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCTTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGAt                                                                                                        <  1:62955/149‑1 (MQ=255)
                                             caacGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGAt                                                                                                        <  6:206540/149‑1 (MQ=255)
                                              cacGCTGCACATAGCTTGCGCAGCCGACAACCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATTCCCCGCGTATCTGCTGCAACCGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATc                                                                                                       <  8:91150/148‑1 (MQ=255)
                                                 gCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGAt                                                                                                        >  1:157405/1‑145 (MQ=255)
                                                        aTAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTgg                                                                                             >  2:271283/1‑149 (MQ=255)
                                                         tAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGc                                                                                            >  1:54249/1‑149 (MQ=255)
                                                               gcgcAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCctgctg                                                                                      >  1:16505/1‑149 (MQ=255)
                                                                  cAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATAAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCaa                                                                                   <  8:255917/149‑1 (MQ=255)
                                                                   aGCCGACAAACATCTCGGCCTCGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCATGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAac                                                                                  <  8:60648/149‑1 (MQ=255)
                                                                   aGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCaa                                                                                   >  3:28767/1‑148 (MQ=255)
                                                                   aGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAac                                                                                  <  5:12107/149‑1 (MQ=255)
                                                                      cGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACg                                                                               <  2:157405/149‑1 (MQ=255)
                                                                       gACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCGGATCAGCTCCTTGGCCTGCTGCAACCCGt                                                                              >  4:133410/1‑149 (MQ=255)
                                                                         cAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTc                                                                            <  2:176800/149‑1 (MQ=255)
                                                                          aaaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCa                                                                           <  6:295027/149‑1 (MQ=255)
                                                                           aaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGTCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAg                                                                          <  7:316466/149‑1 (MQ=255)
                                                                           aaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAg                                                                          >  3:140128/1‑149 (MQ=255)
                                                                             cAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGc                                                                         >  6:144292/1‑148 (MQ=255)
                                                                                cTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGatt                                                                     <  2:114114/149‑1 (MQ=255)
                                                                                   ggCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGattcat                                                                  <  4:28767/149‑1 (MQ=255)
                                                                                   ggCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGattca                                                                   <  1:206951/148‑1 (MQ=255)
                                                                                        cggccggcAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGc                                                             <  7:319208/149‑1 (MQ=255)
                                                                                         ggccggcAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCa                                                            <  1:171164/149‑1 (MQ=255)
                                                                                                  caTCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGc                                                   <  2:16505/149‑1 (MQ=255)
                                                                                                         aTGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTAtctc                                            >  7:204413/1‑149 (MQ=255)
                                                                                                         aTGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTAtctc                                            <  1:271283/149‑1 (MQ=255)
                                                                                                          tGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCg                                           <  4:140128/149‑1 (MQ=255)
                                                                                                           ggCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCgg                                          >  4:164253/1‑149 (MQ=255)
                                                                                                             ccAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTGAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGcc                                        >  7:62828/1‑149 (MQ=255)
                                                                                                                 tGCCCGGCGGATCTGCTGCAACAGCCCCGGCAAGGGCGGCAGTCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGt                                    <  4:233950/149‑1 (MQ=255)
                                                                                                                   ccccGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTc                                   <  1:223610/148‑1 (MQ=255)
                                                                                                                                ctgcAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGc                     <  2:51106/149‑1 (MQ=255)
                                                                                                                                    aaCAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGcc                 >  2:191811/1‑149 (MQ=255)
                                                                                                                                         ccccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCgg                          >  7:33306/1‑135 (MQ=255)
                                                                                                                                         ccccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCgg                          <  8:33306/135‑1 (MQ=255)
                                                                                                                                         ccccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGc            >  4:315597/1‑149 (MQ=255)
                                                                                                                                          cccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGcc           <  1:301926/149‑1 (MQ=255)
                                                                                                                                          cccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGcc           >  1:274449/1‑149 (MQ=255)
                                                                                                                                          cccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGcc           >  8:281260/1‑149 (MQ=255)
                                                                                                                                          cccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGcc           >  8:286400/1‑149 (MQ=255)
                                                                                                                                           ccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCg          >  7:14161/1‑149 (MQ=255)
                                                                                                                                             ggCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGtt        <  5:315920/149‑1 (MQ=255)
                                                                                                                                              gCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGttt       <  2:54249/149‑1 (MQ=255)
                                                                                                                                                aaGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTtc     <  5:144292/149‑1 (MQ=255)
                                                                                                                                                aaGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTtc     <  1:382955/149‑1 (MQ=255)
                                                                                                                                                   ggCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTtcctc  <  3:56023/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTTTTCCTC  >  NC_002947/5685754‑5686049

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