Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,852,188 G→T 5.9% intergenic (+99/+106) PP_5531 → / ← PP_2505 hypothetical protein/GAF domain/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,852,1880GT 5.9% 215.3 / 3.5 71intergenic (+99/+106)PP_5531/PP_2505hypothetical protein/GAF domain/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base T (2/2);  total (33/38)
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.06e-01

GGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCC  >  NC_002947/2852042‑2852332
                                                                                                                                                  |                                                                                                                                                
ggAAAGGCCGATTACGTGCAGCACATGATCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTgg                                                                                                                                                 <  6:350918/148‑1 (MQ=255)
  aaaGGCCGATTACGTGCAGCACATGTTAGCCTATCTGCGACAGTGAATCCTGCACGAGGACCGCGCTAAGCAACTGAGCCGCGCTCTCCACTCGGCGCGGTCGCAGTTTCATTGTCCAGCCATGTGGACGCCATTCACGCCCCTCGtt                                                                                                                                               >  3:135846/1‑148 (MQ=255)
      gCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCaa                                                                                                                                           >  5:23420/1‑148 (MQ=255)
         gATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAg                                                                                                                                       >  5:129944/1‑149 (MQ=255)
           ttACGTGAATCACATGTTCGCCTATCTGGGAGAGTGAATCCTGGACGATGTGGGCGCTAAGCAAGTGTGCCGAGCTCGAGTCTCGGGGGGGTCGCTGTTTGTTTGTCCAGCGATGTGGACTCCATTCCCGCCTCTGGTTGCAATAAggg                                                                                                                                     <  8:164565/149‑1 (MQ=255)
               gTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTTATTGCAACCAGGggc                                                                                                                                   <  3:72492/147‑1 (MQ=255)
               gTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTTATTGCAACCAGGggc                                                                                                                                   >  4:72492/1‑147 (MQ=255)
               gTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAg                                                                                                                                       <  3:79603/143‑1 (MQ=255)
               gTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAg                                                                                                                                       >  4:79603/1‑143 (MQ=255)
                tGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTACCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGggcggc                                                                                                                                >  1:314204/1‑149 (MQ=255)
                 gcagcaCATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTCGACGCCATTCCCGCCCCTGGATGCAATAAGGggcggc                                                                                                                                >  8:304836/1‑148 (MQ=255)
                  cagcaCATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTTATTGCAACCAGGggcggga                                                                                                                               <  7:286355/148‑3 (MQ=255)
                   agcaCATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGAt                                                                                                                             <  8:151951/149‑1 (MQ=255)
                     cacaTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGAt                                                                                                                             >  4:367661/1‑147 (MQ=255)
                     cacaTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGAt                                                                                                                             >  2:331561/1‑147 (MQ=255)
                      acaTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAg                                                                                                                          >  8:354734/1‑149 (MQ=255)
                       caTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGt                                                                                                                         <  8:140146/149‑1 (MQ=255)
                        aTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAg                                                                                                                                       <  4:338762/134‑1 (MQ=255)
                        aTGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAg                                                                                                                                       >  3:338762/1‑134 (MQ=255)
                          gTTCGCCTACCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTg                                                                                                                       >  7:8770/1‑148 (MQ=255)
                           ttCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGa                                                                                                                      <  8:128167/148‑1 (MQ=255)
                             cGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGaaa                                                                                                                    <  5:295807/148‑1 (MQ=255)
                              gcgtATCTGCGAGAGTGAATCCTGCTCGTGGTTGGCGCTAAGCATCTGTGCCGCGCTCTGCTCTGGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCGATTCCGGCCCGTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGc                                                                                                                  <  4:135846/146‑1 (MQ=255)
                              gCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGc                                                                                                                  >  8:77509/1‑149 (MQ=255)
                                cTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGa                                                                                                                <  6:291650/149‑1 (MQ=255)
                                   tCTGCGACAGTGACTCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAg                                                                                                               <  5:329613/147‑1 (MQ=255)
                                       cGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGAt                                                                                                          <  8:8770/148‑1 (MQ=255)
                                        gACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCa                                                                                                        <  2:166064/149‑1 (MQ=255)
                                               aaTCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCaaa                                                                                                  <  1:189855/148‑1 (MQ=255)
                                               aaTCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGATTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAAt                                                                                                 >  8:369999/1‑149 (MQ=255)
                                                aTCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAAt                                                                                                 <  8:252778/148‑1 (MQ=255)
                                                      cACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTaa                                                                                           <  6:285617/148‑1 (MQ=255)
                                                      cACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTaa                                                                                           <  1:331561/148‑1 (MQ=255)
                                                      cACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTa                                                                                            <  8:280803/147‑1 (MQ=255)
                                                          aGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAg                                                                                      <  7:77509/149‑1 (MQ=255)
                                                               ggCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGaa                                                                                 >  7:340340/1‑149 (MQ=255)
                                                                 cgcTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAAtg                                                                               >  1:120475/1‑149 (MQ=255)
                                                                    tAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAAt                                                                                <  5:181742/145‑1 (MQ=255)
                                                                    tAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAAt                                                                                >  6:181742/1‑145 (MQ=255)
                                                                      aGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGAt                                                                           <  6:23420/148‑1 (MQ=255)
                                                                      aGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATg                                                                          <  7:354734/149‑1 (MQ=255)
                                                                      aGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATg                                                                          >  1:349165/1‑149 (MQ=255)
                                                                         aaCTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCGCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGt                                                                       <  5:110715/149‑1 (MQ=255)
                                                                              tgCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCt                                                                  >  8:90889/1‑149 (MQ=255)
                                                                                ccGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTcc                                                                >  3:91027/1‑149 (MQ=255)
                                                                                     cTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGAt                                                           >  1:179537/1‑149 (MQ=255)
                                                                                      tCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAACCGCGCTGATCAGGCAAATTGGTAACTCGGCCGATTGTGATGGGTCCCCTCCCGGagc                                                          >  8:79913/1‑147 (MQ=255)
                                                                                         ccACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCg                                                        >  5:79947/1‑148 (MQ=255)
                                                                                          cACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGga                                                      <  3:388328/149‑1 (MQ=255)
                                                                                                gggCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCtt                                                <  5:219799/149‑1 (MQ=255)
                                                                                                      tCGCAGTTTGTTTGTCCAGCGATGTGGACGCAATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTc                                          <  7:125584/149‑1 (MQ=255)
                                                                                                        gCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCt                                        >  6:167461/1‑149 (MQ=255)
                                                                                                            tttgtttgTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAg                                     >  8:328725/1‑148 (MQ=255)
                                                                                                              tgtttgtCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGa                                   <  7:369999/148‑1 (MQ=255)
                                                                                                                tttgtCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCg                                                              <  7:339478/119‑1 (MQ=255)
                                                                                                                tttgtCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCg                                                              >  8:339478/1‑119 (MQ=255)
                                                                                                                   gtCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAATGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACg                              <  3:34195/148‑1 (MQ=255)
                                                                                                                    tCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCt                            >  3:108775/1‑149 (MQ=255)
                                                                                                                     ccAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTg                           >  4:116860/1‑149 (MQ=255)
                                                                                                                      cAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTgg                          >  5:207769/1‑149 (MQ=255)
                                                                                                                      cAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTgg                          >  1:330485/1‑149 (MQ=255)
                                                                                                                      cAGCGATGTGCACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAATGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCggg                          <  1:230710/149‑1 (MQ=255)
                                                                                                                         cGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCgg                                                       <  2:19102/117‑1 (MQ=255)
                                                                                                                         cGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCgg                                                       >  1:19102/1‑117 (MQ=255)
                                                                                                                           aTGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACAt                     <  3:367661/149‑1 (MQ=255)
                                                                                                                                  cGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCCTAACTTCCTCAGGATAACGCTGGCACATAGACAAc              >  8:325034/1‑149 (MQ=255)
                                                                                                                                     cATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTg           <  2:314204/149‑1 (MQ=255)
                                                                                                                                     cATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTg           <  6:79947/149‑1 (MQ=255)
                                                                                                                                       ttCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGt         <  4:91027/149‑1 (MQ=255)
                                                                                                                                         gccGCCCCTTATTGCAACCAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCa       >  3:175460/2‑149 (MQ=255)
                                                                                                                                         cccGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTc        <  7:304836/148‑1 (MQ=255)
                                                                                                                                              cccTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTcc  <  2:179537/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGAAAGGCCGATTACGTGCAGCACATGTTCGCCTATCTGCGACAGTGAATCCTGCACGAGGTGGGCGCTAAGCAACTGTGCCGCGCTCGCCACTCGGGGCGGTCGCAGTTTGTTTGTCCAGCGATGTGGACGCCATTCCCGCCCCTGGTTGCAATAAGGGGCGGCGATCAGTTGAAAGCGAGCTGATCAGGCAAATTGGTAACTCAGCCGAATGTGATGCGTACCCTCCCGGATCCGGAGATCTTAACTTCCTCAGGATAACGCTGGCACATAGACAACTTGGTCACGTCC  >  NC_002947/2852042‑2852332

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