Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,676,956 A→G 5.9% intergenic (‑111/‑42) PP_3241 ← / → PP_3242 hypothetical protein/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,676,9560AG 5.9% 201.1 / 2.1 68intergenic (‑111/‑42)PP_3241/PP_3242hypothetical protein/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base G (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

ATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCATGCTGGTGGTGAT  >  NC_002947/3676808‑3677104
                                                                                                                                                     |                                                                                                                                                    
aTGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGaa                                                                                                                                                      <  1:101235/149‑1 (MQ=255)
     aaaaCCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACtata                                                                                                                                                 <  4:151027/149‑1 (MQ=255)
     aaaaCCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACtata                                                                                                                                                 <  6:12445/149‑1 (MQ=255)
     aaaaCCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACtata                                                                                                                                                 <  3:178509/149‑1 (MQ=255)
      aaaCCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGCTTCCGGCCGTGCGACGCTGAACtata                                                                                                                                                 <  1:4003/149‑1 (MQ=255)
             cgcgTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGAc                                                                                                                                         <  2:157096/149‑1 (MQ=255)
              gcgTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACg                                                                                                                                        >  3:346269/1‑149 (MQ=255)
                gTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGc                                                                                                                                      >  8:172745/1‑149 (MQ=255)
                          gagaCTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCg                                                                                                                              >  1:162300/1‑147 (MQ=255)
                             aCTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAggg                                                                                                                         >  8:306418/1‑149 (MQ=255)
                                 ttCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAg                                                                                                                     <  7:172745/149‑1 (MQ=255)
                                  tCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGTGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAg                                                                                                                     <  8:309108/148‑1 (MQ=255)
                                     tGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGAt                                                                                                                 >  4:251198/1‑149 (MQ=255)
                                      gCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGAtt                                                                                                                >  1:246315/1‑149 (MQ=255)
                                        ttCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTc                                                                                                               >  1:178001/1‑148 (MQ=255)
                                        ttCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTc                                                                                                               >  4:157804/1‑148 (MQ=255)
                                           aCTCGTTTAATCCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGATGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAAt                                                                                                           <  6:320086/149‑1 (MQ=255)
                                           aCTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAAt                                                                                                           <  3:248710/149‑1 (MQ=255)
                                           aCTCCTTTAAACCGTGGCGCGACCGCAATGTCGCGGACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTACTCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAAt                                                                                                           >  6:374416/1‑149 (MQ=255)
                                             tCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAAtgt                                                                                                         >  8:351219/1‑149 (MQ=255)
                                              ccTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGa                                                                                                                                             <  7:231711/112‑1 (MQ=255)
                                              ccTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGa                                                                                                                                             >  8:231711/1‑112 (MQ=255)
                                               cTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAAtgt                                                                                                         <  2:162300/147‑1 (MQ=255)
                                               cTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTgc                                                                                                       >  8:177628/1‑149 (MQ=255)
                                                tttAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTgcg                                                                                                      >  8:230485/1‑149 (MQ=255)
                                                tttAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTgcg                                                                                                      >  6:330484/1‑149 (MQ=255)
                                                  tAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTgcg                                                                                                      >  2:339579/1‑147 (MQ=255)
                                                        gTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAg                                                                                               >  2:260332/1‑148 (MQ=255)
                                                         tGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGa                                                                                             <  5:354927/149‑1 (MQ=255)
                                                             gcgACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCACGGATTCCAATGTGCGCAGCCAGGaaa                                                                                           >  6:379331/1‑147 (MQ=255)
                                                                aCCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGga                                                                                      <  7:230485/149‑1 (MQ=255)
                                                                aCCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGATTAATCTTTTTTGCCGCCGGTTACGGCCGTGCGACGCTGAACTAAAGTGAATAAGGCGAACGCCGTAGGGCCAGGGATTGCAATGTGCGCAGCCAGGAAACAGga                                                                                      >  1:72579/1‑149 (MQ=255)
                                                                 ccTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGgac                                                                                     >  6:66348/1‑149 (MQ=255)
                                                                   tCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGgacga                                                                                   >  3:359173/1‑149 (MQ=255)
                                                                     aaTGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGgacgac                                                                                  >  1:14526/1‑148 (MQ=255)
                                                                         tCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCACGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATg                                                                              <  5:379331/148‑1 (MQ=255)
                                                                             gTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAGCTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCtt                                                                         >  1:278308/1‑149 (MQ=255)
                                                                             gTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAGCTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCtt                                                                         >  1:283018/1‑149 (MQ=255)
                                                                                cGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTc                                                                      <  4:59307/149‑1 (MQ=255)
                                                                                cGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTc                                                                      <  5:66348/149‑1 (MQ=255)
                                                                                   ttGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTc                                                                      >  4:289293/1‑146 (MQ=255)
                                                                                   ttGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTc                                                                      <  3:289293/146‑1 (MQ=255)
                                                                                   ttGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGcc                                                                   <  2:178001/149‑1 (MQ=255)
                                                                                   ttGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGc                                                                    >  4:12173/1‑148 (MQ=255)
                                                                                      aGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCAt                                                                 <  7:177628/148‑1 (MQ=255)
                                                                                        ggTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAg                                                              >  4:360785/1‑149 (MQ=255)
                                                                                              ggCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAGCTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGa                                                        <  2:283018/149‑1 (MQ=255)
                                                                                              ggCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAGCTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGa                                                        <  2:278308/149‑1 (MQ=255)
                                                                                                  gggAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGc                                                    >  7:123262/1‑149 (MQ=255)
                                                                                                     aaaaGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACAATAGTGAAGACTGCGAACGCCGTAGGGCCAGGGATTCCAAAGTGCGAAGCCAGGAAACAGGACGACGATGGGCTTATCGCCATCAGTGGCGAGGGCcgc                                                 >  7:343365/1‑149 (MQ=255)
                                                                                                       aaGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCgcgc                                               >  4:93844/1‑149 (MQ=255)
                                                                                                           tttAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTg                                           >  7:34924/1‑149 (MQ=255)
                                                                                                              aaTCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGgcg                                        >  1:95138/1‑149 (MQ=255)
                                                                                                               aTCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGgcg                                        <  2:95138/148‑1 (MQ=255)
                                                                                                                tCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCt                                      >  4:297428/1‑149 (MQ=255)
                                                                                                                  ttttttGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGATACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGt                                    >  2:160232/1‑149 (MQ=255)
                                                                                                                    ttttGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCt                                      <  3:12173/145‑1 (MQ=255)
                                                                                                                         ccgccgGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCt                             <  3:251198/149‑1 (MQ=255)
                                                                                                                         ccgccgGTGCCGGCCGTGCGACTCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTGCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCt                             <  1:260332/149‑1 (MQ=255)
                                                                                                                          cgcggGTTCCGGCCGGGCGACGGTGAACTATAGTGAAGACTGCGAACGCCGTAGGGCCAGGTATTCCAATGTGCTCAGCCAGGAAACAGGACGTCGATGGTCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCt                             <  4:74039/148‑1 (MQ=255)
                                                                                                                          cgccgGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTg                            >  1:44620/1‑149 (MQ=255)
                                                                                                                              ggTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTAc                        <  5:374416/149‑1 (MQ=255)
                                                                                                                                 ttCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTg                      <  2:14526/149‑1 (MQ=255)
                                                                                                                                      gCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCAtt                 <  4:359173/149‑1 (MQ=255)
                                                                                                                                       ccGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTc                <  5:330484/149‑1 (MQ=255)
                                                                                                                                         gTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCAt              <  7:306418/149‑1 (MQ=255)
                                                                                                                                         gTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCAt              <  7:351219/149‑1 (MQ=255)
                                                                                                                                            cGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCATGCt           <  2:44620/149‑1 (MQ=255)
                                                                                                                                                 cTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCATGCtggtgg      <  8:34924/149‑1 (MQ=255)
                                                                                                                                                    aaCTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGCTTCCAATGTGCGCAGCCAGGAAACAGGACGCCGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCATGCGGGTGGTGa   >  7:75761/1‑149 (MQ=255)
                                                                                                                                                     aCTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCATGCTGGTGGTGa   <  1:339579/148‑1 (MQ=255)
                                                                                                                                                     aCTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCATGCTGGTGGTGAt  <  4:316399/149‑1 (MQ=255)
                                                                                                                                                     aCTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCATGCTGGTGGTGAt  >  6:358282/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGGAAAAACCAACGCGTCTGATGAAGAGACTATTCATGCTTCACTCCTTTAAACCGTGGCGCGACCTCAATGTCGCGTACGGTTGAGGGTGAAGGCGGGGAAAAGGTTTAATCTTTTTTGCCGCCGGTTCCGGCCGTGCGACGCTGAACTATAGTGAAGACGGCGAACGCCGTAGGGCCAGGGATTCCAATGTGCGCAGCCAGGAAACAGGACGACGATGGGCTTTTCGCCATCAGTGGCGAGGGCCGCGCCGTGGCGCTGTTCCGCCTGGTACTGGCATTCATGCTGGTGGTGAT  >  NC_002947/3676808‑3677104

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: