Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,698,797 Δ1 bp 6.6% coding (503/1014 nt) 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,7970A. 6.6% 233.4 / 10.4 76coding (503/1014 nt)cheBCchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base A (25/46);  new base . (3/2);  total (28/48)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

TGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCC  >  NC_002947/1698650‑1698944
                                                                                                                                                    |                                                                                                                                                    
tgctgcGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAg                                                                                                                                                    <  1:957346/150‑1 (MQ=255)
 gctgcGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGt                                                                                                                                                   >  2:1181818/1‑150 (MQ=255)
  ctgcGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGtg                                                                                                                                                  <  2:533286/150‑1 (MQ=255)
   tgcGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGtgc                                                                                                                                                 <  1:206298/150‑1 (MQ=255)
      gcAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGtgctgc                                                                                                                                              <  1:1113057/150‑1 (MQ=255)
            tCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGTCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAgg                                                                                                                                        <  1:274705/150‑1 (MQ=255)
            tCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAgg                                                                                                                                        <  2:215107/150‑1 (MQ=255)
            tCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAgg                                                                                                                                        <  1:190044/150‑1 (MQ=255)
            tCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAgg                                                                                                                                        <  1:1286603/150‑1 (MQ=255)
              cTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGt                                                                                                                                      >  2:345323/1‑150 (MQ=255)
              cTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGt                                                                                                                                      >  2:389082/1‑150 (MQ=255)
                          tggctggTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCg                                                                                                                          <  1:159539/150‑1 (MQ=255)
                             ctggTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGTGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGttt                                                                                                                       <  2:820586/150‑1 (MQ=255)
                              tggTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTc                                                                                                                      >  2:429407/1‑150 (MQ=255)
                                     ccAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCggcggc                                                                                                               >  2:1011594/1‑150 (MQ=255)
                                       agcagcGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCa                                                                                                             >  1:1239504/1‑150 (MQ=255)
                                         cagcGTGCACCTGCCGCAAGGTTGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATc                                                                                                           <  2:192056/150‑1 (MQ=255)
                                            cGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTg                                                                                                        <  1:1101950/150‑1 (MQ=255)
                                             gTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGc                                                                                                       <  2:116512/150‑1 (MQ=255)
                                               gCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGCCTTGTTGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTg                                                                                                     <  1:978847/150‑1 (MQ=255)
                                                 aCCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGt                                                                                                   >  2:1147720/1‑150 (MQ=255)
                                                  ccTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTc                                                                                                  >  1:893852/1‑150 (MQ=255)
                                                   cTGCCGCAAGGTCGGTGGCTGCGCCGATGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGGAGAAGTCGAAGGGCTGCGCAAGGGTTTGTCGGCAGCGGTTCAGGCGGCGATATGGCTGGTcg                                                                                                 >  2:542013/1‑149 (MQ=255)
                                                    tGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCa                                                                                                <  1:429407/150‑1 (MQ=255)
                                                        gCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAt                                                                                            <  2:869545/150‑1 (MQ=255)
                                                        gCAAGGTCGGTGGCTGCACAGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAt                                                                                            >  2:86506/1‑150 (MQ=255)
                                                            ggTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCAGCAGCACTTCGAGtgctgccgg                                                                                                                                           >  1:856987/1‑96 (MQ=255)
                                                            ggTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCAGCAGCACTTCGAGtgctgccgg                                                                                                                                           <  2:856987/99‑4 (MQ=255)
                                                                 gTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTg                                                                                                     <  2:54300/132‑1 (MQ=255)
                                                                 gTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTg                                                                                                     >  1:54300/1‑132 (MQ=255)
                                                                   ggCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgt                                                                                <  2:1147468/150‑1 (MQ=255)
                                                                      tGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGtt                                                                             >  1:192868/1‑150 (MQ=255)
                                                                      tGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGtt                                                                             >  1:993775/1‑150 (MQ=255)
                                                                      tGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCAgcacga                                                                             <  1:1147720/150‑6 (MQ=255)
                                                                         aCCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCCGCGt                                                                                                                         <  2:342967/104‑1 (MQ=255)
                                                                         aCCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGt                                                                                                                         >  1:342967/1‑104 (MQ=255)
                                                                          ccGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGcc                                                                         <  2:668646/150‑1 (MQ=255)
                                                                              tGCGCGAGGCGTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                     <  1:1019626/150‑1 (MQ=255)
                                                                                 gcgAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCa                                                                  <  1:1011594/150‑1 (MQ=255)
                                                                                   gAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATg                                                                <  1:441755/150‑1 (MQ=255)
                                                                                    aGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATgg                                                               >  1:1142158/1‑150 (MQ=255)
                                                                                     ggCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGc                                                              <  1:185922/150‑1 (MQ=255)
                                                                                     cgatctccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGTTAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGTTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGc                                                              <  1:1042768/144‑1 (MQ=255)
                                                                                      gCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGcc                                                             <  2:192868/150‑1 (MQ=255)
                                                                                        ttCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAtg                                                                                           >  2:428593/1‑119 (MQ=255)
                                                                                        ttCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAtg                                                                                           <  1:428593/119‑1 (MQ=255)
                                                                                          cccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                         <  2:1211706/150‑1 (MQ=255)
                                                                                           ccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATg                                                        <  2:1189443/150‑1 (MQ=255)
                                                                                           ccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATg                                                        >  2:322396/1‑150 (MQ=255)
                                                                                            cAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCCCTCGAAGTGCTGCTCAAAGGTTTGCCGGCAGCGTTTCCGGCGGCCCTCGTGCTGGTCCAGCATGTGGACGAGGTGTTTGCCGCAGGCATGGCCGAAtcg                                                       <  2:771105/150‑3 (MQ=255)
                                                                                            cAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCCTCGTGCTGGTCCAGCATGTTGACCAGGTGTTTGCCGCAGGCATGGccagatcc                                                       >  2:1042768/1‑144 (MQ=255)
                                                                                              gcgcCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCt                                                     <  2:1025948/150‑1 (MQ=255)
                                                                                              gcgcCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCt                                                     <  2:993775/150‑1 (MQ=255)
                                                                                               cgccgcGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGTCATGGCCGAATGGCTc                                                    <  1:869274/150‑1 (MQ=255)
                                                                                                ggcgcgGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTca                                                   <  1:903048/148‑1 (MQ=255)
                                                                                                      ggTTGGTGGCCATCGGGTCTTCGGCAGGTGGCCCGGCATCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAgcagcg                                             <  2:1239504/150‑1 (MQ=255)
                                                                                                             ggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAgg                                      <  1:345323/150‑1 (MQ=255)
                                                                                                             ggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAgg                                      <  2:1142158/150‑1 (MQ=255)
                                                                                                              gCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGc                                     >  1:924459/1‑150 (MQ=255)
                                                                                                              gCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGc                                     <  1:1181818/150‑1 (MQ=255)
                                                                                                                  tCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGc                                 >  1:170261/1‑150 (MQ=255)
                                                                                                                  tCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGc                                 >  2:69014/1‑150 (MQ=255)
                                                                                                                    gggTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCg                               <  1:980356/150‑1 (MQ=255)
                                                                                                                      gTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGt                             <  1:69014/150‑1 (MQ=255)
                                                                                                                      gTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGt                             >  2:187644/1‑150 (MQ=255)
                                                                                                                               gggTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGc                    <  2:548309/150‑1 (MQ=255)
                                                                                                                                  tGGCCCGACAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCCGCGTTTGCCGCGGCCATCGGGCTGGCCCAGCATGTGGACCAGGTGTTTGGCGCAGACATGGCCGAATGGCTCAGGAGCGCAGCAGGCCTGCCGGGGCGCCTGGCCCg                 >  2:1215438/1‑150 (MQ=255)
                                                                                                                                   ccttgagcagcaCTTCGAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGCCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGt                >  2:501719/7‑150 (MQ=255)
                                                                                                                                    cttgagcagcaCTTCGAGTGCTGCCCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                         >  2:322405/6‑108 (MQ=255)
                                                                                                                                    cttgagcagcaCTTCGAGTGCTGCCCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                         <  1:322405/103‑1 (MQ=255)
                                                                                                                                     cccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGa              >  2:634442/1‑150 (MQ=255)
                                                                                                                                      ccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGaa             <  1:86506/150‑1 (MQ=255)
                                                                                                                                       cGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAg            <  2:999552/150‑1 (MQ=255)
                                                                                                                                        ggCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAgg           >  1:930817/1‑150 (MQ=255)
                                                                                                                                        ggCAGCACTCGAAGTGCTGCTCAAGGGGTGCCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAgg           <  1:324935/150‑1 (MQ=255)
                                                                                                                                                  cGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAccgcc  <  2:930817/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCTGCGCAAGATCCTCAATATCGGCTGGCTGGTCGGCCAGCAGCGTGCACCTGCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCC  >  NC_002947/1698650‑1698944

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: