Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,698,755 A→G 7.7% E168G (GAA→GGA) ‡ cheBC → chemotaxis response regulator protein‑glutamate methylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,698,7550AG 7.7% 182.9 / 5.3 65E168G (GAA→GGA) ‡cheBCchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base G (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

GCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACC  >  NC_002947/1698609‑1698899
                                                                                                                                                     |                                                                                                                                                
gctgcGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGaa                                                                                                                                                 <  3:21699/148‑1 (MQ=255)
     gcAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGtgctg                                                                                                                                           <  4:371729/149‑1 (MQ=255)
      cAAGATCCTAAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGtgctgc                                                                                                                                          >  6:229365/1‑149 (MQ=255)
       aaGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCGGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTTCGAGtgctgcc                                                                                                                                         <  2:357440/149‑2 (MQ=255)
       aaGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTgctgct                                                                                                                                         >  6:105084/1‑149 (MQ=255)
                aaTATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGGTGGGGGTTGGGTGGGGAAAACGGGTCTACGGCGGGTGGCCCGCCAGCCGACGAAGTGGTGATCAAggg                                                                                                                                   >  1:312054/1‑149 (MQ=255)
                 atatCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGtgctgc                                                                                                                                          >  6:350527/1‑138 (MQ=255)
                 atatCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGtgctgc                                                                                                                                          <  5:350527/138‑1 (MQ=255)
                 atatCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGAT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTg                                                                                                                               >  7:261342/1‑149 (MQ=255)
                  tatCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGc                                                                                                                              <  4:197526/149‑1 (MQ=255)
                                   gCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCcggcg                                                                                                             >  6:258650/1‑149 (MQ=255)
                                        cagcGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCg                                                                                                                            >  6:108475/1‑129 (MQ=255)
                                        cagcGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCg                                                                                                                            <  5:108475/129‑1 (MQ=255)
                                           cGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGt                                                                                                     <  4:231239/149‑1 (MQ=255)
                                           cGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGt                                                                                                     <  2:87810/149‑1 (MQ=255)
                                                 ccTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAg                                                                                                                                                <  4:177875/100‑1 (MQ=255)
                                                 ccTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAg                                                                                                                                                >  3:177875/1‑100 (MQ=255)
                                                  cTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTc                                                                                              >  4:364245/1‑149 (MQ=255)
                                                  cTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTc                                                                                              >  7:211629/1‑149 (MQ=255)
                                                         aaGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAtg                                                                                       <  5:105084/149‑1 (MQ=255)
                                                         aaGGGCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAtg                                                                                       <  3:235226/149‑1 (MQ=255)
                                                             tCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGCTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGa                                                                                   >  2:67520/1‑149 (MQ=255)
                                                               ggTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGGTCAAGGGTTTGCCGGGAGCGTTTCGGGCGGCCTGGGTGCTGGTCCGGCATGTgggg                                                                                  >  8:313320/1‑146 (MQ=255)
                                                               ggTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGAc                                                                                  <  6:290454/148‑1 (MQ=255)
                                                                   gCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgt                                                                             <  2:312054/149‑1 (MQ=255)
                                                                    cTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgtg                                                                            >  5:314421/1‑149 (MQ=255)
                                                                       cACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgtgt                                                                           <  3:364245/147‑1 (MQ=255)
                                                                       cACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGttt                                                                         <  5:258650/149‑1 (MQ=255)
                                                                        aCCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTg                                                                        <  1:14786/149‑1 (MQ=255)
                                                                            cagcgcGAGGCTTCCCAGCGCCCCGGGTT‑GGT‑‑GGGAATCGGGTCTTAGGCGGGTGGCCCCGCAGCACACGAATTGCTGCTAAAGGGTTTTCCGGCACCGTATGCGGCCGCCCTCGTGCTGGTCCAGCCTGTGCACCACGTGTTTgaccc                                                                    >  2:110371/3‑145 (MQ=255)
                                                                               cgcgAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCa                                                                   <  1:23308/147‑1 (MQ=255)
                                                                                gcgAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                  <  6:144795/147‑1 (MQ=255)
                                                                                gcgAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                  >  5:144795/1‑147 (MQ=255)
                                                                                     gCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGc                                                           <  8:211629/149‑1 (MQ=255)
                                                                                     gCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGCGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGc                                                           >  7:39049/1‑149 (MQ=255)
                                                                                          ctcgcgcAGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                      <  6:174005/146‑1 (MQ=255)
                                                                                          ccAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                      >  4:92200/1‑149 (MQ=255)
                                                                                                     ggTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCagcagc                                           >  1:315005/1‑149 (MQ=255)
                                                                                                           gt‑‑ggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAgcagcgcagc                                      <  8:261342/149‑1 (MQ=255)
                                                                                                            t‑‑ggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGttt                                                                         <  4:329552/113‑1 (MQ=255)
                                                                                                            t‑‑ggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGttt                                                                         >  3:329552/1‑113 (MQ=255)
                                                                                                                  aaTCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGcc                                 >  4:259486/1‑149 (MQ=255)
                                                                                                                  aaTCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGcc                                 >  6:68606/1‑149 (MQ=255)
                                                                                                                       ggTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCg                            >  2:185861/1‑149 (MQ=255)
                                                                                                                          cTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGc                                                   <  5:375108/123‑1 (MQ=255)
                                                                                                                          cTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGc                                                   >  6:375108/1‑123 (MQ=255)
                                                                                                                           ttCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCAGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTgc                        >  1:281244/1‑149 (MQ=255)
                                                                                                                             cggcggGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTTCCGGTg                         <  3:92200/146‑1 (MQ=255)
                                                                                                                             cggcggGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTgc                        <  6:314421/147‑1 (MQ=255)
                                                                                                                                   gTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTgc                        <  4:340509/141‑1 (MQ=255)
                                                                                                                                   gTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTgc                        >  3:340509/1‑141 (MQ=255)
                                                                                                                                   gTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGcc                >  8:13522/1‑149 (MQ=255)
                                                                                                                                   gTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGcc                >  8:19069/1‑149 (MQ=255)
                                                                                                                                     ggCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCg              >  1:228665/1‑149 (MQ=255)
                                                                                                                                       ttgagcagcaCTTCGAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                  <  1:191356/91‑1 (MQ=255)
                                                                                                                                       ttgagcagcaCTTCGAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                  >  2:191356/5‑95 (MQ=255)
                                                                                                                                       cccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTg            <  5:229365/149‑1 (MQ=255)
                                                                                                                                        ccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAgcagcgc                                         <  5:126431/119‑1 (MQ=255)
                                                                                                                                        ccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAgcagcgc                                         >  6:126431/1‑119 (MQ=255)
                                                                                                                                         gagcagcaCTTCGAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTgcgc                      <  3:281556/135‑1 (MQ=255)
                                                                                                                                         gagcagcaCTTCGAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTgcgc                      >  4:281556/3‑137 (MQ=255)
                                                                                                                                           gcagcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAgg        <  3:259486/149‑1 (MQ=255)
                                                                                                                                              gcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCg      >  7:110848/1‑148 (MQ=255)
                                                                                                                                                 cTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGccc               <  1:157500/136‑1 (MQ=255)
                                                                                                                                                 cTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGccc               >  2:157500/1‑136 (MQ=255)
                                                                                                                                                 cTCGAAGTGCTGCTCAAGGGGTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAcc  <  6:349029/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTT‑GGT‑‑GGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACC  >  NC_002947/1698609‑1698899

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