Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,045,688 C→T 6.2% T301M (ACG→ATG)  PP_1819 → methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,045,6880CT 6.2% 194.4 / 2.6 64T301M (ACG→ATG) PP_1819methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base T (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GGCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGGCGCCGCGCCGGGATGAGG  >  NC_002947/2045548‑2045834
                                                                                                                                             |                                                                                                                                                  
ggCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTgg                                                                                                                                            >  5:88332/1‑149 (MQ=255)
 gCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTg                                                                                                                                             <  1:181651/147‑1 (MQ=255)
  cccGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCg                                                                                                                                          <  6:88332/149‑1 (MQ=255)
      cgAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                                      <  4:310146/149‑1 (MQ=255)
      cgAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                                      <  4:32314/149‑1 (MQ=255)
      cgAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                                      <  4:340172/149‑1 (MQ=255)
             gggCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCgg                                                                                                                               >  1:200142/1‑149 (MQ=255)
             gggCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAACGCCAGCGCCGT‑GGCCAGGCCGATGGCGATGGCATTg                                                                                                                                       <  5:359313/141‑1 (MQ=37)
             gggCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAACGCCAGCGCCGT‑GGCCAGGCCGATGGCGATGGCATTg                                                                                                                                       >  6:359313/1‑141 (MQ=37)
              ggCTGCAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGc                                                                                                                              >  4:278451/1‑149 (MQ=255)
              ggCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGc                                                                                                                              <  7:97516/149‑1 (MQ=255)
                        gcAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTggcgg                                                                                                                    >  7:15342/1‑149 (MQ=255)
                            ttGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTg                                                                                                                >  6:86639/1‑149 (MQ=255)
                                   ccAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCAGCGCCGT‑GGCCAGGCCGATGGCGATGGCATTGcg                                                                                                                                     <  7:325860/121‑2 (MQ=25)
                                   ccAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCAGCGCCGT‑GGCCAGGCCGATGGCGATGGCATTGcg                                                                                                                                     >  8:325860/1‑120 (MQ=25)
                                     aGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTAc                                                                                                       >  6:231731/1‑149 (MQ=255)
                                       cgcCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTAccc                                                                                                     <  2:307884/149‑1 (MQ=255)
                                             gCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCtg                                                                                               <  5:231731/149‑1 (MQ=255)
                                                    aCCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCg                                                                                         <  2:200142/148‑1 (MQ=255)
                                                     ccGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCgg                                                                                        <  1:6086/148‑1 (MQ=255)
                                                     ccGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGc                                                                                       >  8:238061/1‑149 (MQ=255)
                                                          cGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGtt                                                                                  <  2:173925/149‑1 (MQ=255)
                                                             gtCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCg                                                                                         >  1:63503/1‑139 (MQ=255)
                                                             gtCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCg                                                                                         <  2:63503/139‑1 (MQ=255)
                                                              tCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                                      >  2:178272/1‑93 (MQ=255)
                                                              tCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                                      <  1:178272/93‑1 (MQ=255)
                                                              tCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGaa                                                                               <  7:70251/148‑1 (MQ=255)
                                                                gACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCt                                                                                                                 >  2:183338/1‑112 (MQ=255)
                                                                 aCCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCAt                                                                           >  5:328250/1‑149 (MQ=255)
                                                                      gCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGAt                                                                      <  6:66322/149‑1 (MQ=255)
                                                                      gCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGAt                                                                      <  6:328250/149‑1 (MQ=255)
                                                                          aTCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGGTTGTCGGCGTGTTGGCGGCCTGGGTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCTGGTGATCGcc                                                                  >  6:63399/1‑149 (MQ=255)
                                                                                 gCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGccc                                                           <  7:239732/149‑1 (MQ=255)
                                                                                 gCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGcc                                                            >  1:321583/1‑148 (MQ=255)
                                                                                    ggAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAgc                                                        <  6:209266/149‑1 (MQ=255)
                                                                                        cAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCAt                                                                           <  4:144544/126‑1 (MQ=255)
                                                                                        cAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCAt                                                                           >  3:144544/1‑126 (MQ=255)
                                                                                        cAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCAt                                                    <  3:278451/149‑1 (MQ=255)
                                                                                           gCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATc                                                                     >  1:181438/1‑129 (MQ=255)
                                                                                           gCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATc                                                                     <  2:181438/129‑1 (MQ=255)
                                                                                           ctccTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCAtt                                                   <  5:268246/145‑1 (MQ=255)
                                                                                                gCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCTCCCGCGCCCAGCGCATTGcggcgg                                            <  7:52665/149‑1 (MQ=255)
                                                                                                gCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGcggcgg                                            <  5:63399/149‑1 (MQ=255)
                                                                                                     gCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGAGTACCCGCGCTGTGGGCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCCGCGCATTGCGGCGGGTGaa                                       >  4:208373/1‑149 (MQ=255)
                                                                                                            gCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTggc                                                                                                                      >  8:159830/1‑63 (MQ=255)
                                                                                                            gCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTggc                                                                                                                      <  7:159830/63‑1 (MQ=255)
                                                                                                             cAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCtt                               >  8:161496/1‑149 (MQ=255)
                                                                                                              aGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTc                              >  3:256623/1‑149 (MQ=255)
                                                                                                               gcgcAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTc                              >  2:50978/1‑148 (MQ=255)
                                                                                                               gcgcAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCg                             <  2:321583/149‑1 (MQ=255)
                                                                                                                     tGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGccc                       <  7:161496/149‑1 (MQ=255)
                                                                                                                      gccatcgccatcGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGcc                        >  5:29763/1‑147 (MQ=255)
                                                                                                                         atcgccatcgGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTg                                    <  2:201161/132‑1 (MQ=255)
                                                                                                                         atcgccatcgGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTg                                    >  1:201161/1‑132 (MQ=255)
                                                                                                                          tcgccatcgGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGgcg                  >  3:37525/1‑149 (MQ=255)
                                                                                                                              catcgGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGgcgccgc              >  3:220862/1‑149 (MQ=255)
                                                                                                                                   gCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCTCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGGCGCCGCGCCgg         <  8:15342/149‑1 (MQ=255)
                                                                                                                                   gCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGACCTGGGCTTCGACGCCCAGGCGCCGCGCCgg         >  4:273308/1‑149 (MQ=255)
                                                                                                                                   gCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGGCGCCGCGCCgg         >  4:77278/1‑149 (MQ=255)
                                                                                                                                     cTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGAGTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGGCGCCGCGCCGGGa       <  3:208373/149‑1 (MQ=255)
                                                                                                                                        gCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCt                                <  7:163666/121‑1 (MQ=255)
                                                                                                                                        gCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCt                                >  8:163666/1‑121 (MQ=255)
                                                                                                                                         cc‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGGCGCCGCGCCGGGATGAg   >  5:336962/1‑149 (MQ=255)
                                                                                                                                          c‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGAATACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGGCGCCGCGCCGGGATGAgg  >  6:311623/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GGCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGGCGCCGCGCCGGGATGAGG  >  NC_002947/2045548‑2045834

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: ‑