Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,557,924 T→G 7.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,557,9240TG 7.4% 171.0 / 4.2 66intergenic (‑295/+36)PP_4888/purAmethyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base G (2/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

CTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCG  >  NC_002947/5557780‑5558069
                                                                                                                                                |                                                                                                                                                 
cTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAgg                                                                                                                                               >  3:216046/1‑149 (MQ=255)
  aCAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGATGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGCTCAGGGc                                                                                                                                             >  4:116918/1‑149 (MQ=255)
    aaCGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcg                                                                                                                                           >  8:331683/1‑149 (MQ=255)
        gCCTGCTCCCAACCATCGTTGGGGGCCCACATCGGAGCGGTTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCCgcggc                                                                                                                                       <  6:81266/149‑1 (MQ=37)
           ggcgcCAAACCATCTTTGTGTCCCCAAATCGGCGGGTTTTTACCGGCTAATACGGTAGCGGCGGCAACGGTGTATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTc                                                                                                                                    <  4:203501/148‑1 (MQ=255)
           ggcgcCAAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCCGCGGCGtt                                                                                                                                    <  8:277279/148‑2 (MQ=38)
             cgcCCAACCACCGTTGTGGGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCa                                                                                                                                   <  5:188453/148‑1 (MQ=255)
                ccAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGAt                                                                                                                                >  2:102473/1‑148 (MQ=255)
                   aCCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGc                                                                                                                            >  4:368450/1‑149 (MQ=255)
                         gTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                      <  1:34956/149‑1 (MQ=255)
                         gTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                      <  7:331683/149‑1 (MQ=255)
                          ttGTTTGCCAAAATAGGTGGGGTTTAACCCGCAAATACGGTAGCGGCGGAAACGGTGGATTACCGGTGTATTTTGGCCAGCACAGACATCCCCACAAATCAAAACGCCGCGGCCCTTCTCTGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                      <  2:200665/148‑1 (MQ=255)
                              gtgCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAgg                                                                                                                    >  4:31232/1‑146 (MQ=255)
                                  ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGc                                                                                                                   >  5:82274/1‑143 (MQ=255)
                                  ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGc                                                                                                                   <  6:82274/143‑1 (MQ=255)
                                  ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGAt                                                                                                               >  8:113151/1‑147 (MQ=255)
                                  ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGAt                                                                                                               >  5:1232/1‑147 (MQ=255)
                                   cacaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGACCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAg                                                                                                            >  5:30291/1‑149 (MQ=255)
                                    acaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAgg                                                                                                           <  4:216046/149‑1 (MQ=255)
                                           cggggttttACCCGCTAATAGGGTAGCGGCGGCAACGGTGTATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACgg                                                                                                    <  6:322334/143‑1 (MQ=255)
                                            gCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTATCACGGCCGCGGAGTAACATAAACCAGGAAGGCAGCTAAGGCGAACggg                                                                                                   >  6:59171/1‑149 (MQ=37)
                                             cGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAg                                                                                                                 >  8:106301/1‑134 (MQ=255)
                                             cGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAg                                                                                                                 <  7:106301/134‑1 (MQ=255)
                                                gTTTACCCTCTAATACGGTAGCTGCGGCAACGGTTGATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGATAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCTGCCTTCCAGGCTGATct                                                                                                             <  8:292448/135‑2 (MQ=255)
                                                gTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCTGCCTTCCAGGCTGATct                                                                                                             >  7:292448/1‑134 (MQ=38)
                                                     cccGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcag                                                                                           >  7:78393/1‑148 (MQ=255)
                                                            aTACGGTAGGGGCGGAAACGGTGGATTACCGGTGTATATTGGCCACCACATACATGCCCACAAATCAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGAt                                                                                                                                <  6:10359/104‑1 (MQ=255)
                                                            aTACGGTAGCGGCGGCAACGGTGGATTACCGGTGTATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATg                                                                                    <  2:273163/148‑1 (MQ=255)
                                                            aTACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGAt                                                                                                                                >  5:10359/1‑104 (MQ=255)
                                                                ggTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCACGGCCGCGGCGTCAGATCCCCCAGGAAGGCAGCTCAGGCGAACGGGTGACGCAGCACGAAGGtcac                                                                               >  3:367819/1‑147 (MQ=255)
                                                                    gcggcgGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTg                                                                           <  1:102473/149‑1 (MQ=255)
                                                                     cggcggCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGc                                                                          >  1:21779/1‑149 (MQ=255)
                                                                      ggcggcAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCg                                                                         <  5:127770/149‑1 (MQ=255)
                                                                       gcggaAACGGTGTATACCCGGTGTATATTGCCCAGAACTGACATGCCAACAAATGAAAACGCCGCGGCCCTTATCCGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCgg                                                                        <  3:6042/149‑1 (MQ=255)
                                                                          gcAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTcc                                                                     <  7:113151/149‑1 (MQ=255)
                                                                           cAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCg                                                                    <  6:165305/149‑1 (MQ=255)
                                                                           cAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGAGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCg                                                                    <  5:59171/149‑1 (MQ=255)
                                                                                   gATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCTGCCTTCCAGGCTGATct                                                                                                             >  3:52905/1‑99 (MQ=255)
                                                                                     ttACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCg                                                           >  5:297895/1‑148 (MQ=255)
                                                                                        ccGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCGGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAat                                                       >  1:266718/1‑149 (MQ=255)
                                                                                           gTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAAGGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgat                                                    <  3:31232/149‑1 (MQ=255)
                                                                                           gTGGATATGGCCCAGCACAGACATCCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatgat                                                    <  3:38667/149‑1 (MQ=255)
                                                                                              gATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTc                                                 >  1:165778/1‑149 (MQ=255)
                                                                                                     gCCAGCACAGACATGCCCACAAATGAAAACGCCCCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGATCGACATCATTTCGACCGgc                                          <  3:230880/149‑1 (MQ=255)
                                                                                                      ccAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCAGGGCCGGTCGAAATGATGTCGATCGgcg                                         >  4:106196/1‑149 (MQ=255)
                                                                                                       cAGCACAGACATGCCCACAAATGAAAACGCCGCGACCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcgc                                        <  6:30291/149‑1 (MQ=255)
                                                                                                             aGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAAc                                  <  7:352521/149‑1 (MQ=255)
                                                                                                               aCATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCa                                <  2:280991/149‑1 (MQ=255)
                                                                                                               aCATGCCAAAAATTCATAACGCCGCGGCCCTTCTCAGGGCGGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGTGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAAcc                                 <  4:367819/148‑1 (MQ=255)
                                                                                                                 aTGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGc                              >  1:329934/1‑149 (MQ=255)
                                                                                                                    cccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGGCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCCCGATGGGCTCGTTGCGGGCCGGGGCGGTCGAAATGATGTCGCTCGGCGCGGCAACCAGctct                           >  1:361834/1‑149 (MQ=255)
                                                                                                                    cccACAAATGAAAACGCCGCGGCCCTTCCCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGGCTCGTTGCGGTCCGGGCCGGGCGAAATGATGTCGATCGGCGCGCCAACCAGctct                           <  3:116918/149‑1 (MQ=255)
                                                                                                                     ccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCa                                    <  3:333947/139‑1 (MQ=255)
                                                                                                                     ccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCa                                    >  4:333947/1‑139 (MQ=255)
                                                                                                                       acaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTc                         >  4:134817/1‑148 (MQ=255)
                                                                                                                         aaaTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGAt                      <  3:134817/149‑1 (MQ=255)
                                                                                                                         aaaTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGAt                      <  3:106196/149‑1 (MQ=255)
                                                                                                                                 aCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa              <  3:368450/149‑1 (MQ=255)
                                                                                                                                 aCGCCGCGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa              >  4:123906/1‑149 (MQ=255)
                                                                                                                                   gccgcGGCCCTTCTCAGGGCCGCGGCGTAAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATg            >  8:268145/1‑149 (MQ=255)
                                                                                                                                   gccgcGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCtt                          <  8:349864/135‑1 (MQ=255)
                                                                                                                                   gccgcGGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCtt                          >  7:349864/1‑135 (MQ=255)
                                                                                                                                    ccgcGGCCCTTCTCAGGACCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGt           >  7:275802/1‑149 (MQ=255)
                                                                                                                                        ggCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCGTGactacatc       <  8:78393/149‑8 (MQ=255)
                                                                                                                                        ggCCCGTCTCAGGGCCGCGGCGTCAGATCAGCCGGGAAGGCAGATCAGGCGAACGGGTGACGCACCACGATGGTATCGTTGCGGTCCGGGGCGGGCGCAATGATGGCGAGCGGCGCGCCCACCCGCTCGTCGATGCGGGTGATGGAGt        >  8:296326/1‑148 (MQ=255)
                                                                                                                                          cccTTCTCAGGGCCGCGGCGTAAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc     <  7:268145/149‑1 (MQ=255)
                                                                                                                                               ctcAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcgcg  >  1:179786/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCG  >  NC_002947/5557780‑5558069

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: