Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,562,923 C→A 6.1% intergenic (‑294/+37) PP_4888 ← / ← purA methyl‑accepting chemotaxis transducer/adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,562,9230CA 6.1% 870.3 / 30.6 269intergenic (‑294/+37)PP_4888/purAmethyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (134/118);  new base A (9/8);  total (143/126)
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.99e-01

Only 100 of 281 total aligned reads displayed.

CCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCGCA  >  NC_002947/5562774‑5563071
                                                                                                                                                     |                                                                                                                                                    
ccGCTCCTACACCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtc                                                                                                                                                      <  2:6697991/150‑1 (MQ=255)
ccGCTCCTACACCGGCCGGCGCCCAACCTTGGTTTTGTGCCCACATCGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACATACATGCCCACAAATGAAAACGCCGCGGCCCTtc                                                                                                                                                      <  2:3735204/150‑1 (MQ=38)
ccGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtc                                                                                                                                                      >  1:1171736/1‑150 (MQ=255)
ccGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtc                                                                                                                                                      >  1:8205058/1‑150 (MQ=255)
ccGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtc                                                                                                                                                      >  2:5525938/1‑150 (MQ=255)
  gCTCATACACCGGCCTGCGCCCAACAATCTTTGTGTGCCCACTTTGGAGGGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtctc                                                                                                                                                    <  1:3436427/150‑1 (MQ=38)
   cTCCTAAAACGGCCTGCGCCAAAAAATCTTTGTGTGCCCACATAGGAGCGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCa                                                                                                                                                   <  1:1427135/150‑1 (MQ=38)
   cTCCTAAAACGCCCTGCGCAAAACAATGTTTTTGTGCACAAATCGGAGGGTTTTCACCCGCCAATACGGTAGCGGCGGCACCGGTGGATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCa                                                                                                                                                   <  1:5348276/150‑1 (MQ=255)
    accTACACCGGCGTGCGCCAAACCTTGTTTGTGTGCCAAATTTGTTGGGGTTTTCCCCGCTAATACGGTAGCGGCGGCAAGGGTGTATTACCGGTGTATTTTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAg                                                                                                                                                  <  1:4451471/149‑1 (MQ=255)
        aCAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGcc                                                                                                                                              >  1:6076085/1‑150 (MQ=255)
        aCAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGcc                                                                                                                                              >  1:8276455/1‑150 (MQ=255)
         cAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcg                                                                                                                                             >  2:8488618/1‑150 (MQ=255)
          aaCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATCGGCGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgc                                                                                                                                            <  2:2654175/150‑1 (MQ=255)
          aaCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgc                                                                                                                                            >  1:4672692/1‑150 (MQ=255)
          aaCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgc                                                                                                                                            >  2:5137358/1‑150 (MQ=255)
           ccggccggcGCCCAACCAGGGTGGTTTGCCCACATTGGTGGGGGTTTACCCGCGAATAGGGTAGCGGCGGCAACGGTGGATTACGGGTGGATTTTGGCCTGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCcgcg                                                                                                                                           <  2:829296/149‑1 (MQ=25)
           aCGGCGGGCGCCAAACCTTGGTTGTGTGCCAACTTTGGAGGGGTTTTACCGGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCcgcg                                                                                                                                           <  1:8512970/150‑1 (MQ=25)
            cGGCCTGCGCCCAACCTTCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCgcgg                                                                                                                                          <  2:4035155/150‑1 (MQ=255)
             ggCCTGCGCCCAACCATGGTTGTGTGCCCACTTAGGAGCGGGTTTACCCGCGATTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCgcggc                                                                                                                                         <  2:6559411/150‑1 (MQ=255)
             ggCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCgcggc                                                                                                                                         <  1:2132519/150‑1 (MQ=255)
                cTGCGCCCAACCATGGTTGTGTGCCCACTTAGGAGCGGGTTTACCCGCTAATACGGTAGCGGCGGCAAGGGTGGATTACCGGTGGTTATTGGCCAGCACAGACCTGCCCACAAATGAAAACGCCGCGGCCCTTCTCCGGGCCGCGGCGTc                                                                                                                                      <  2:13650/150‑1 (MQ=255)
                cTGCGCCCAACCATGGTTGTGTGCCCACATTGGTGGGGTTTTACCCGCGAATACGGTAGCGGCGGAAACGGTGGATTACGGGTGGATATTGGCCAGCACAGACATGCCCAAAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTc                                                                                                                                      <  1:7104760/150‑1 (MQ=38)
                cTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTc                                                                                                                                      >  2:5615419/1‑150 (MQ=255)
                 tGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCa                                                                                                                                     >  2:7736855/1‑150 (MQ=255)
                   cgcCAAACCATTTTGGTGTGCCAAATTTGGTGGGGGTTTACCCGGGAATACGGTAGCGGGGGCAACGGTGTTTTACGGGTGTTTATTGGCCAGCACAGACATGCCCAAAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGa                                                                                                                                   <  2:751513/150‑1 (MQ=255)
                       cAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAg                                                                                                                               >  2:4841910/1‑150 (MQ=255)
                       cAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAg                                                                                                                               >  2:3615188/1‑150 (MQ=255)
                        aaCCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGc                                                                                                                              >  2:2943043/1‑150 (MQ=255)
                            aTCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTgg                                                                                                                          >  2:5880425/1‑150 (MQ=255)
                            aTCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTgg                                                                                                                          >  1:3238864/1‑150 (MQ=255)
                             tCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGa                                                                                                                         >  1:4954770/1‑150 (MQ=255)
                              cGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACGGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                        <  1:4338130/150‑1 (MQ=255)
                               gTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAg                                                                                                                       >  1:2849615/1‑150 (MQ=255)
                                ttGTGGGCCCACATCGGCGCGGGTTTACCCGCGAATACGGGAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAAGGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAgg                                                                                                                      <  2:494372/150‑1 (MQ=38)
                                  gtgtgCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCa                                                                                                                    <  2:7201651/150‑1 (MQ=255)
                                     tgCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGAt                                                                                                                 >  2:451667/1‑150 (MQ=255)
                                       cccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCa                                                                                                               >  1:2370353/1‑150 (MQ=255)
                                        ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAg                                                                                                              >  1:4545676/1‑150 (MQ=255)
                                         cacaTCGGAGGGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCTGCCTTCCAGGCTGATctga                                                                                                             <  2:5439278/150‑4 (MQ=38)
                                         cacaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAgg                                                                                                             >  2:8770586/1‑150 (MQ=255)
                                          acaTCGGCGGGGTTTTCCCCGCTAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGTATATTGGCCTGCACAGACATGCCCACAAATGAAACCGCCGCGGCCCTTCTCCGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGc                                                                                                            <  2:787258/150‑1 (MQ=255)
                                          acaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGc                                                                                                            <  2:7947568/150‑1 (MQ=255)
                                               ttcgCGGGTAAACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACACATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGAGCAGGCGAACg                                                                                                       >  2:5212429/4‑150 (MQ=39)
                                                gAGGGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACgg                                                                                                      <  2:5332795/150‑1 (MQ=255)
                                                    gggTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGa                                                                                                  >  2:2817103/1‑150 (MQ=255)
                                                       tttACCCGCGACTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgc                                                                                               <  2:3354919/150‑1 (MQ=255)
                                                          ccccGCGAATACGGTCGCGGCGGCAAGGGTGGATTACGGTTGTATTTTGGCCAGAACAGACTTGCCCACAAAGCAAAACGCCGCGGCCCTTCTCCGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcagc                                                                                            <  1:3065276/149‑1 (MQ=255)
                                                          aCCCGCGAATAGGGTAGCGGGGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACTGACATGCCCACAATTGAAAACGCCGCGGCCCTTCTCCGGGCGGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcagc                                                                                            <  1:7006649/150‑1 (MQ=255)
                                                            cggcgATTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                          <  1:7229482/148‑1 (MQ=255)
                                                            ccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                          >  1:3528635/1‑150 (MQ=255)
                                                            ccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                          >  2:6119017/1‑150 (MQ=255)
                                                            ccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCGGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                          >  1:2659544/1‑150 (MQ=255)
                                                            ccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATAAGCCAGGAAGGCAGATCAGGCGAAAGGGTGACGCAGCAc                                                                                          >  1:395874/1‑150 (MQ=255)
                                                                 aaTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCCCGGCGTCAGATCAGCCTGGAAGGCACATCAGGCGAACGGGTGACGCAGCACCATgg                                                                                     >  1:2430402/1‑150 (MQ=255)
                                                                     cGGTCGCGGCGGCAACGGTGGATTACGGGTGTATATTGGCCAGCACAGACATGCCCACAAATCAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGtctc                                                                                 <  1:1911472/150‑1 (MQ=255)
                                                                     cGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGtctc                                                                                 <  2:2588201/150‑1 (MQ=255)
                                                                        tAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGtt                                                                              <  2:150555/150‑1 (MQ=255)
                                                                           cggcggCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCg                                                                           >  1:7183272/1‑150 (MQ=255)
                                                                           cggcggCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCg                                                                           <  2:3264527/150‑1 (MQ=255)
                                                                             gcggcAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGt                                                                         >  1:4259608/1‑150 (MQ=255)
                                                                              cggcACGGGTGTTTTACGGGTGTATTTGGCCCAGCCCCGACCGCCCCACAAAGGAAAACGCCGCGGCCCTTCTCTGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTc                                                                        <  1:3690640/150‑1 (MQ=255)
                                                                                            tACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGaaa                                                          <  2:2140559/150‑1 (MQ=255)
                                                                                            tACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGaaa                                                          <  1:8048449/150‑1 (MQ=255)
                                                                                               cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatga                                                       >  1:6775683/1‑150 (MQ=255)
                                                                                                ggTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgat                                                      <  2:5736872/150‑1 (MQ=255)
                                                                                                 gTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgatg                                                     <  2:7138924/150‑1 (MQ=255)
                                                                                                     atatTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGa                                                 >  2:426775/1‑150 (MQ=255)
                                                                                                      tttgggCCACCACAGACATGCCCACAAATCAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGATCGACATCATTTCGac                                                <  1:5148142/146‑2 (MQ=255)
                                                                                                      tatTGGCCAGCACAGACATGCCCACAATTGAAAACGCCGCGGCCCGGATAAGGGCCGCGGCGTCAGATCAGCCTTGAAGGCAGATCAGGCGAACGGGTGACGCAGCCCCCGGGGGACGTTGCGGTCCGGGCCGGTCGTGATGATGTCGAt                                                >  1:5925315/1‑150 (MQ=255)
                                                                                                      tatTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGAt                                                >  1:7971687/1‑150 (MQ=255)
                                                                                                      tatTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCGTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGAt                                                >  1:1717445/1‑150 (MQ=255)
                                                                                                        ttGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCg                                              <  2:7005858/150‑1 (MQ=255)
                                                                                                        ttGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCg                                              >  1:7665802/1‑150 (MQ=255)
                                                                                                          ggCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGCACGGGTGACGCAGCACGATGGTCTCGTTGCGGGCCGGGCCGGTCGAAATGATGTCGATCGgc                                            >  1:190029/1‑150 (MQ=255)
                                                                                                           gCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcg                                           >  2:3904966/1‑150 (MQ=255)
                                                                                                           gCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcg                                           >  2:5129985/1‑150 (MQ=255)
                                                                                                           gCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcg                                           >  2:2477764/1‑150 (MQ=255)
                                                                                                            ccAGCACAGACATGCCCACAACTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcgc                                          <  2:6584903/150‑1 (MQ=255)
                                                                                                            ccAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcgc                                          <  1:6512481/150‑1 (MQ=255)
                                                                                                               gCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGcc                                       <  2:399096/150‑1 (MQ=255)
                                                                                                                     aCATGCCCACAAATGAAAACGCCGCGGGCCATCACAGGGCCGCGGCGTCAGAACAGCCAGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAAATATGTCGATCGGCGCGCCAACCAg                                 >  1:5271706/1‑150 (MQ=255)
                                                                                                                       aTGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGct                               >  2:1777162/1‑150 (MQ=255)
                                                                                                                       aTGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGct                               >  2:7249758/1‑150 (MQ=255)
                                                                                                                               aaaTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATg                       <  2:6256906/150‑1 (MQ=255)
                                                                                                                               aaaTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATg                       >  2:1353156/1‑150 (MQ=255)
                                                                                                                                      gaCGCCGCGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa                >  2:3760649/2‑150 (MQ=255)
                                                                                                                                      caccccgcGGCCCGTCTCCGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa                <  1:6408780/146‑1 (MQ=255)
                                                                                                                                       aCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGAt               <  1:1064607/150‑1 (MQ=255)
                                                                                                                                       aCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGAt               <  2:6864085/150‑1 (MQ=255)
                                                                                                                                         gccgcGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGt             >  2:6908020/1‑150 (MQ=255)
                                                                                                                                          ccgcGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTa            >  2:1334789/1‑150 (MQ=255)
                                                                                                                                           cgcgGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAg           >  1:2570066/1‑150 (MQ=255)
                                                                                                                                            gcgGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGt          >  1:7848042/1‑150 (MQ=255)
                                                                                                                                             cgGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCATTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTc         >  1:1401076/1‑150 (MQ=255)
                                                                                                                                               gCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc       >  2:3485678/1‑150 (MQ=255)
                                                                                                                                                cccTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcg      >  2:5056951/1‑150 (MQ=255)
                                                                                                                                                 ccTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcgc     <  2:3238864/150‑1 (MQ=255)
                                                                                                                                                 ccTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcgc     >  1:7141008/1‑150 (MQ=255)
                                                                                                                                                 ccTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcgc     >  2:3520055/1‑150 (MQ=255)
                                                                                                                                                  cTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcgcg    >  1:5660565/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCGCA  >  NC_002947/5562774‑5563071

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: