Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,562,924 T→G 6.4% intergenic (‑295/+36) 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,9240TG 6.4% 771.2 / 31.0 269intergenic (‑295/+36)PP_4888/purAmethyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base T (133/116);  new base G (9/9);  total (143/126)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

Only 100 of 280 total aligned reads displayed.

GCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCGCAGC  >  NC_002947/5562776‑5563073
                                                                                                                                                    |                                                                                                                                                     
 cTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATTGGAGGGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCa                                                                                                                                                     <  2:4773505/150‑1 (MQ=255)
 cTCCTACAACGGCCTGCGCCAAACCATGGTTGTGTGCCCACTTAGGCGGGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCa                                                                                                                                                     <  2:6613413/150‑1 (MQ=38)
 cTCCTAAAACGCCCTGCGCAAAACAATGTTTTTGTGCACAAATCGGAGGGTTTTCACCCGCCAATACGGTAGCGGCGGCACCGGTGGATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCa                                                                                                                                                     <  1:5348276/150‑1 (MQ=255)
  tCCTACAACGGCCTGCGCCCACCCATCGTTGTGTGCCCACATCGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAg                                                                                                                                                    <  2:6209198/150‑1 (MQ=255)
   ccTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAgg                                                                                                                                                   <  2:2827923/150‑1 (MQ=255)
    cTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAggg                                                                                                                                                  >  1:7924770/1‑150 (MQ=255)
      aCAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGcc                                                                                                                                                >  1:6076085/1‑150 (MQ=255)
       cAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcg                                                                                                                                               >  2:5549927/1‑150 (MQ=255)
        aaCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgc                                                                                                                                              >  2:5173086/1‑150 (MQ=255)
        aaCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgc                                                                                                                                              >  1:4672692/1‑150 (MQ=255)
        aaCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgc                                                                                                                                              >  1:7124884/1‑150 (MQ=255)
         ccGGCCTGCGCCCAACCATGGTTGTGTGCCCACATAGGAGGGGTTTTACCGGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCcgcg                                                                                                                                             <  2:2159695/149‑1 (MQ=38)
          cGGCCTGCGCCCAACCTTCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCgcgg                                                                                                                                            <  2:4035155/150‑1 (MQ=255)
          cGGCCTGCGCCCAACCAGCGTGGGGTGCCCACATTGGCGGGGGTTTACCGGCGAATACGGTAGCGGCGGCACCGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCCgcgg                                                                                                                                            <  2:2650915/150‑1 (MQ=37)
           ggCCTGCGCCCAACCATGGTTGTGTGCCCACTTAGGAGCGGGTTTACCCGCGAATACGGTAGCGTCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCCgcggc                                                                                                                                           <  2:5843952/150‑1 (MQ=38)
           ggCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCgcggc                                                                                                                                           <  2:2757262/150‑1 (MQ=255)
              cTGCGCCCAACCATGGTTGTGTGCCCACTTAGGAGCGGGTTTACCCGCTAATACGGTAGCGGCGGCAAGGGTGGATTACCGGTGGTTATTGGCCAGCACAGACCTGCCCACAAATGAAAACGCCGCGGCCCTTCTCCGGGCCGCGGCGTc                                                                                                                                        <  2:13650/150‑1 (MQ=255)
              cTGCGCCCAACCATGGTTGTGTGCCCACATTGGTGGGGTTTTACCCGCGAATACGGTAGCGGCGGAAACGGTGGATTACGGGTGGATATTGGCCAGCACAGACATGCCCAAAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTc                                                                                                                                        <  1:7104760/150‑1 (MQ=38)
               tGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCa                                                                                                                                       >  2:3293907/1‑150 (MQ=255)
                 cgcCAACCCATGGTTGTGTGCCCACATAGGAGGGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGa                                                                                                                                     <  2:1848823/150‑1 (MQ=38)
                    ccAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACCGACATGCCCACAAATGAAAACGCCGCGGCCCGTCTCAGGGCCGCGGCGTCAGATCa                                                                                                                                  >  1:1213811/1‑150 (MQ=255)
                       aCCATCGTTGTGTGCCCACTTCGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGcc                                                                                                                               <  2:8428528/150‑1 (MQ=255)
                       aCCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACACGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGcc                                                                                                                               >  1:4280338/1‑150 (MQ=255)
                          aTCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTgg                                                                                                                            >  2:341980/1‑150 (MQ=255)
                           tCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGa                                                                                                                           >  1:4954770/1‑150 (MQ=255)
                            ggTTGTGGGCCCACATCGGCGCGGGTTTACCCGCGAATACGGTCGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                          <  1:627842/149‑1 (MQ=255)
                            cGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACGGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                          <  1:4338130/150‑1 (MQ=255)
                              ttGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAgg                                                                                                                        >  2:6400472/1‑150 (MQ=255)
                              ttGTGGGCCCACATCGGCGCGGGTTTACCCGCGAATACGGGAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAAGGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAgg                                                                                                                        <  2:494372/150‑1 (MQ=38)
                               tgtgtgCCCACTTTGGAGCGGGTTTACCGGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGc                                                                                                                       <  2:6021194/150‑1 (MQ=255)
                                gtgtgCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCa                                                                                                                      <  2:7201651/150‑1 (MQ=255)
                                      ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAg                                                                                                                >  1:4545676/1‑150 (MQ=255)
                                       cacaTCGGAGGGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCTGCCTTCCAGGCTGATctga                                                                                                               <  2:5439278/150‑4 (MQ=38)
                                       cacaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAgg                                                                                                               >  2:8770586/1‑150 (MQ=255)
                                        acaTCGGCGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGc                                                                                                              <  2:3969049/150‑1 (MQ=255)
                                              gAGGGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACgg                                                                                                        <  2:5332795/150‑1 (MQ=255)
                                                gCGGGTTGACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATGGGCCAGCACAGACATGCCCACAAATGTAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGt                                                                                                      <  2:5258039/150‑1 (MQ=255)
                                                   ggTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGAc                                                                                                   >  2:3703093/1‑150 (MQ=255)
                                                        aCCCGCGAATAGGGTAGCGGGGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACTGACATGCCCACAATTGAAAACGCCGCGGCCCTTCTCCGGGCGGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcagc                                                                                              <  1:7006649/150‑1 (MQ=255)
                                                          cggcgATTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                            <  1:7229482/148‑1 (MQ=255)
                                                          ccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                            >  2:4914102/1‑150 (MQ=255)
                                                          ccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                            >  2:3255733/1‑150 (MQ=255)
                                                          ccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                            >  1:1992691/1‑150 (MQ=255)
                                                          ccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCGGGAAGGCAGATCAGGCGAACGGGTGACGCAGCAc                                                                                            >  1:2659544/1‑150 (MQ=255)
                                                             ccaaTAGGGTAGCGGCGGAAAGGGTGTATTACGGGTGTATATTGGCCACCACAGACTTGCCAACAAATCAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGAt                                                                                         <  1:3322254/148‑1 (MQ=255)
                                                               aaTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCTGGCGAACGGGTGACGCAGCACGATgg                                                                                       >  2:1046417/1‑150 (MQ=255)
                                                               aaTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCATCTCAGGGCCGCGGCGTCAGATCAGCCAGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATgg                                                                                       >  1:2430228/1‑150 (MQ=255)
                                                                  aCGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGtct                                                                                    <  2:150348/150‑1 (MQ=255)
                                                                   cGGTCGCGGCGGCAACGGTGGATTACGGGTGTATATTGGCCAGCACAGACATGCCCACAAATCAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGtctc                                                                                   <  1:1911472/150‑1 (MQ=255)
                                                                   cGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGtctc                                                                                   <  2:6985244/150‑1 (MQ=255)
                                                                         cggcggCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCg                                                                             >  1:7183272/1‑150 (MQ=255)
                                                                         cggcggCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCg                                                                             <  1:3563651/150‑1 (MQ=255)
                                                                           gcggcAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATAAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGt                                                                           >  1:2367195/1‑150 (MQ=255)
                                                                                         ttACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATTAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGaa                                                             <  2:3876664/150‑1 (MQ=255)
                                                                                           aCCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAat                                                           <  2:8266930/150‑1 (MQ=255)
                                                                                             cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatga                                                         >  1:6775683/1‑150 (MQ=255)
                                                                                             cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatga                                                         <  1:3516804/150‑1 (MQ=255)
                                                                                              gttgtatatGGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgat                                                        <  1:3211429/145‑1 (MQ=255)
                                                                                               gTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgatg                                                       <  2:7138924/150‑1 (MQ=255)
                                                                                                gggATATTGGCCAGCACAGACCGGCCCACACCTGAACACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCCGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCCCGTTGCGGTCCGGGCCGGTCGAAATGATGt                                                      <  2:752422/149‑1 (MQ=255)
                                                                                                   atatTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGa                                                   >  2:249378/1‑150 (MQ=255)
                                                                                                   atatTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGa                                                   >  2:426775/1‑150 (MQ=255)
                                                                                                   atatTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCGTCTCAGGGCCGCGGCGTCAGATAAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGa                                                   >  1:3797073/1‑150 (MQ=255)
                                                                                                    tttgggCCACCACAGACATGCCCACAAATCAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGATCGACATCATTTCGac                                                  <  1:5148142/146‑2 (MQ=255)
                                                                                                    tatTGGCCAGCACAGACATGCCCACAATTGAAAACGCCGCGGCCCGGATAAGGGCCGCGGCGTCAGATCAGCCTTGAAGGCAGATCAGGCGAACGGGTGACGCAGCCCCCGGGGGACGTTGCGGTCCGGGCCGGTCGTGATGATGTCGAt                                                  >  1:5925315/1‑150 (MQ=255)
                                                                                                    tatTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATAAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGAt                                                  >  1:1501696/1‑150 (MQ=255)
                                                                                                      ttGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCg                                                >  1:7665802/1‑150 (MQ=255)
                                                                                                      ttGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCATCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCg                                                >  1:2550516/1‑150 (MQ=255)
                                                                                                        ggCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgc                                              >  1:1418485/1‑150 (MQ=255)
                                                                                                        ggCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgc                                              <  2:6773370/150‑1 (MQ=255)
                                                                                                         gCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcg                                             >  2:5129985/1‑150 (MQ=255)
                                                                                                          ccAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcgc                                            <  1:6512481/150‑1 (MQ=255)
                                                                                                           cAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcgcg                                           >  2:6874053/1‑150 (MQ=255)
                                                                                                             gCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGcc                                         <  2:399096/150‑1 (MQ=255)
                                                                                                                caGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAAc                                      <  2:7786087/150‑1 (MQ=255)
                                                                                                                     aTGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGct                                 >  2:7249463/1‑150 (MQ=255)
                                                                                                                         ccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTc                             <  1:2838896/150‑1 (MQ=255)
                                                                                                                         ccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCCCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACCATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCACCTCGTc                             >  1:8377463/1‑150 (MQ=255)
                                                                                                                          cacaAATGAAAACGCCGCGGCCCGTCTCAGGGCCGCGGCGGCAGATCAGCCGGGAAGGCAGATCAGGCGAACCGGGGACGCAGCACGATGGTCTCGGTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCCACCCGCGCTTCg                            >  1:4411670/1‑150 (MQ=255)
                                                                                                                           acaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGa                           >  1:1892210/1‑150 (MQ=255)
                                                                                                                             aaaTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATg                         >  2:1353038/1‑150 (MQ=255)
                                                                                                                              aaTGAAAACGCCGCGGCCCGTCTCAGGGCCGCGGCGTCCGATCAGCCTGGAATGCAGATCAGGCGAACGGGGGACGCAGCCCGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCACCGCTTCGATgc                        >  1:3974911/1‑150 (MQ=255)
                                                                                                                                tGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATgcgc                      <  2:658299/150‑1 (MQ=255)
                                                                                                                                  aaacCGCCGCGGCCCTCCTCAGGGCCGCGGCGTCATATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATTATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCtt                    <  1:5173086/150‑1 (MQ=255)
                                                                                                                                    aaCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa                  <  2:5728326/150‑1 (MQ=255)
                                                                                                                                     aCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGAt                 <  2:7924770/150‑1 (MQ=255)
                                                                                                                                       gccgcGGCCCTTCTCAGGGCCGCGGCGTCATAAAAGCCAGGAAGGCAGATCAGGCGAACGGGTGACGCACCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCATCGATGCCCATGATGt               >  2:299357/1‑150 (MQ=255)
                                                                                                                                       gccgcGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGt               >  1:2947677/1‑150 (MQ=255)
                                                                                                                                        ccgcGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTa              <  2:3571140/150‑1 (MQ=255)
                                                                                                                                         cgcgGCCCGTCACAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAg             <  1:8247234/150‑1 (MQ=255)
                                                                                                                                          gcgGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGt            >  2:1720400/1‑150 (MQ=255)
                                                                                                                                          gcgGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGt            >  1:7848042/1‑150 (MQ=255)
                                                                                                                                             gCCCTTCTCAGGGCCGCGGCGTCAGTTCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc         >  2:3486326/1‑150 (MQ=255)
                                                                                                                                              cccTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcg        >  2:5056951/1‑150 (MQ=255)
                                                                                                                                                cTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcgcg      >  1:8019554/1‑150 (MQ=255)
                                                                                                                                                cTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcgcg      >  1:5660565/1‑150 (MQ=255)
                                                                                                                                                    tcAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCGCAGc  >  2:7400787/1‑150 (MQ=255)
                                                                                                                                                    tcAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCGCAGc  >  1:875847/1‑150 (MQ=255)
                                                                                                                                                    tcAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCGCAGc  >  2:8549509/1‑150 (MQ=255)
                                                                                                                                                    tcAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCGCAGc  >  1:1013448/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCGCAGC  >  NC_002947/5562776‑5563073

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: