Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 895,210 G→C 6.0% intergenic (+50/+14) PP_0777 → / ← PP_0778 glutathione peroxidase/group 1 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947895,2100GC 6.0% 272.7 / 4.8 83intergenic (+50/+14)PP_0777/PP_0778glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (35/43);  new base C (2/3);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAG  >  NC_002947/895063‑895358
                                                                                                                                                     |                                                                                                                                                    
tggtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCgcggg                                                                                                                                                     >  5:41410/1‑149 (MQ=255)
  gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGTTGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCgcgggc                                                                                                                                                    <  5:15095/148‑1 (MQ=255)
  gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCgcgggcg                                                                                                                                                   >  5:386113/1‑149 (MQ=255)
  gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCgcgggcg                                                                                                                                                   <  6:200288/149‑1 (MQ=255)
      gCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGgcgc                                                                                                                                               >  2:210379/1‑149 (MQ=255)
       cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCgcgggcgg                                                                                                                                                  >  2:247293/1‑145 (MQ=255)
       cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCgcgggcgg                                                                                                                                                  <  1:247293/145‑1 (MQ=255)
       cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                              <  8:100697/149‑1 (MQ=255)
        aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTg                                                                                                                                             >  2:209375/1‑149 (MQ=255)
        aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTg                                                                                                                                             >  8:11430/1‑149 (MQ=255)
          ggATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGa                                                                                                                                           >  6:55222/1‑149 (MQ=255)
           gATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGAt                                                                                                                                          <  1:53876/149‑1 (MQ=255)
            aTGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTGCAGCGCCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGAt                                                                                                                                          <  1:257052/149‑1 (MQ=25)
               gCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                      <  4:11791/149‑1 (MQ=255)
                 cGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                      >  1:19842/1‑147 (MQ=255)
                 cGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGGGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTcaca                                                                                                                                    <  7:347968/149‑1 (MQ=255)
                       cTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCt                                                                                                                              >  7:291130/1‑149 (MQ=255)
                           cgcgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGAt                                                                                                                                          <  5:178334/133‑1 (MQ=38)
                           cgcgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGAt                                                                                                                                          >  6:178334/1‑133 (MQ=38)
                              gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCTGGCGGCGCTGGATCTCACAGGCGCTGCAAACa                                                                                                                       <  8:117162/149‑1 (MQ=255)
                              gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACa                                                                                                                       <  6:386113/149‑1 (MQ=255)
                                    cgccgcGCACGGCGCCGGATGACCCAGCGGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                                  >  1:23001/1‑149 (MQ=255)
                                    cgccgcGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                                  <  2:23001/148‑1 (MQ=255)
                                     gccgcGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgc                                                                                                                   <  7:345766/146‑1 (MQ=255)
                                      ccgcGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcgc                                                                                                                 >  4:204603/1‑147 (MQ=255)
                                      ccgcGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcgc                                                                                                                 <  3:204603/147‑1 (MQ=255)
                                              ggCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCa                                                                                                       <  4:176094/149‑1 (MQ=255)
                                                 gcCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGcc                                                                                                    >  5:305872/1‑149 (MQ=255)
                                                   cGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCg                                                                                                   >  8:286639/1‑148 (MQ=255)
                                                   cGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGt                                                                                                  <  1:13827/149‑1 (MQ=255)
                                                   cGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGt                                                                                                  <  6:41410/149‑1 (MQ=255)
                                                   cGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGt                                                                                                  <  8:291130/149‑1 (MQ=255)
                                                          cccAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGgcg                                                                                                                                >  3:312944/1‑112 (MQ=255)
                                                          cccAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGgcg                                                                                                                                <  4:312944/112‑1 (MQ=255)
                                                          cccAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCAc                                                                                             >  1:125305/1‑147 (MQ=255)
                                                          cccAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCAc                                                                                             >  7:243449/1‑147 (MQ=255)
                                                           ccAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCGCCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAggtg                                                                                                                                <  5:308051/111‑3 (MQ=21)
                                                           ccAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCGCCTGTGAGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAggtg                                                                                                                                >  6:308051/1‑109 (MQ=21)
                                                           ccAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                            >  2:40619/1‑115 (MQ=255)
                                                           ccAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                            <  1:40619/115‑1 (MQ=255)
                                                              gC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGtt                                                                                       >  3:209740/1‑149 (MQ=255)
                                                                  gtggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGcc                                                                                    >  6:357022/1‑149 (MQ=255)
                                                                  gtggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGcc                                                                                    <  7:11430/149‑1 (MQ=255)
                                                                  gtggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGcc                                                                                    <  5:103208/149‑1 (MQ=255)
                                                                   tggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCg                                                                                   <  5:55222/149‑1 (MQ=255)
                                                                     gtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCgg                                                                                  >  6:337827/1‑148 (MQ=255)
                                                                     gtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGc                                                                                 <  7:286639/149‑1 (MQ=255)
                                                                       gCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGTAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCAt                                                                               >  1:157517/1‑149 (MQ=255)
                                                                        cAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATg                                                                              >  4:57267/1‑149 (MQ=255)
                                                                          ggCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATgcg                                                                            >  6:348010/1‑149 (MQ=255)
                                                                           gCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                                  >  1:39201/1‑110 (MQ=255)
                                                                           gCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                                  <  2:39201/110‑1 (MQ=255)
                                                                           gCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATgcgc                                                                           <  1:209375/149‑1 (MQ=255)
                                                                             cATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTg                                                                         >  1:184042/1‑149 (MQ=255)
                                                                                 gagaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCa                                                                     <  4:27937/149‑1 (MQ=255)
                                                                                   gaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGTAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGa                                                                   <  2:157517/149‑1 (MQ=255)
                                                                                     aaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATa                                                                 >  4:59326/1‑149 (MQ=255)
                                                                                      aGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAc                                                                <  6:305872/149‑1 (MQ=255)
                                                                                      aGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAc                                                                <  2:19842/149‑1 (MQ=255)
                                                                                       gCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAcc                                                               <  5:337827/149‑1 (MQ=255)
                                                                                        cccTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAccc                                                              >  4:197830/1‑149 (MQ=255)
                                                                                              ccTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGaa                                                        <  2:184042/149‑1 (MQ=255)
                                                                                               cTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAc                                                                <  5:26745/140‑1 (MQ=255)
                                                                                               cTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAc                                                                >  6:26745/1‑140 (MQ=255)
                                                                                               cTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAAc                                                       >  7:195628/1‑149 (MQ=255)
                                                                                                tGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCTCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACt                                                      >  2:7364/1‑149 (MQ=255)
                                                                                                 gAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTg                                                     >  1:256611/1‑149 (MQ=255)
                                                                                                  aaGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGt                                                    <  3:57267/149‑1 (MQ=255)
                                                                                                   agagTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGtt                                                   <  2:125305/149‑1 (MQ=255)
                                                                                                    gagTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTc                                                  >  5:186266/1‑149 (MQ=255)
                                                                                                          cTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAAt                                            <  1:210379/149‑1 (MQ=255)
                                                                                                             aCGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATc                                           <  7:215937/147‑1 (MQ=255)
                                                                                                             aCGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATc                                           <  4:209740/147‑1 (MQ=255)
                                                                                                              cGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCTCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGcc                                        <  1:7364/149‑1 (MQ=255)
                                                                                                                    gTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAg                                  >  5:106582/1‑149 (MQ=255)
                                                                                                                          ggCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTg                             >  4:350933/1‑149 (MQ=255)
                                                                                                                                  tgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCt                       <  5:357022/147‑1 (MQ=255)
                                                                                                                                  tgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTgg                     >  1:188726/1‑149 (MQ=255)
                                                                                                                                    agaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCa                   >  5:228746/1‑149 (MQ=255)
                                                                                                                                    agaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCa                   <  3:59326/149‑1 (MQ=255)
                                                                                                                                       tCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGc                <  3:350933/149‑1 (MQ=255)
                                                                                                                                         gAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCt               <  8:195628/148‑1 (MQ=255)
                                                                                                                                             gCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTcgc          <  2:188726/149‑1 (MQ=255)
                                                                                                                                              ccGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGc                                     <  3:113727/121‑1 (MQ=255)
                                                                                                                                              ccGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTcgc          <  2:256611/148‑1 (MQ=255)
                                                                                                                                              ccGCGCGGGCGGCGCTGGATCTCACAGGAGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGc                                     >  4:113727/1‑121 (MQ=255)
                                                                                                                                               cgcgcgGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGGAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAcc                                                               <  4:23681/94‑1 (MQ=255)
                                                                                                                                               cgcgcgGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGGAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAcc                                                               >  3:23681/1‑94 (MQ=255)
                                                                                                                                                gcgcgGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCa       <  5:348010/149‑1 (MQ=255)
                                                                                                                                                     ggcggcGCTGGATCTCACAGGCGCTGCAAACATTGCCCAGAACGCCCAGCCGAAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCAATCAg  >  4:321720/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGC‑GGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTT‑TGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAG  >  NC_002947/895063‑895358

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