Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,557,923 C→A 10.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,557,9230CA10.1% 158.5 / 6.9 66intergenic (‑294/+37)PP_4888/purAmethyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base A (3/4);  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 = 6.11e-01

GCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCG  >  NC_002947/5557776‑5558069
                                                                                                                                                   |                                                                                                                                                  
gCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGCAATGCCCACAAATGAAAACGCCGCGGCCCAtc                                                                                                                                                    >  4:6042/1‑148 (MQ=255)
    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                                                                                                                   <  6:82274/143‑1 (MQ=255)
                                      ccACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGc                                                                                                                   >  5:82274/1‑143 (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                                                                                                                 <  7:106301/134‑1 (MQ=255)
                                                 cGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAg                                                                                                                 >  8:106301/1‑134 (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)
                                                                                                                                                   |                                                                                                                                                  
GCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGCG  >  NC_002947/5557776‑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: ‑