Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,894 T→G 9.6% intergenic (‑200/‑52) tktA ← / → sahR transketolase A/methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,657,8940TG 9.6% 205.3 / 13.6 73intergenic (‑200/‑52)tktA/sahRtransketolase A/methionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base G (4/3);  total (38/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 = 8.37e-01

CGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCG  >  NC_002947/5657750‑5658042
                                                                                                                                                |                                                                                                                                                        
cGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACtttt                                                                                                                                                       <  2:56617/148‑1 (MQ=255)
 gCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                                      <  2:144058/148‑1 (MQ=255)
 gCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                                      >  5:46514/1‑148 (MQ=255)
  cccACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                                      <  5:212290/147‑1 (MQ=255)
  cccACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                                      >  6:212290/1‑147 (MQ=255)
         ggggCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTTGATATTgg                                                                                                                                              <  1:373637/148‑1 (MQ=255)
              tcTTGGGGAAAACATTTAAGCGCCTAGCGCCTGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGTGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTTCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGc                                                                                                                                        <  8:203987/148‑1 (MQ=255)
               cTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGc                                                                                                                                        >  5:412886/1‑148 (MQ=255)
                     gCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTCGATATTgg                                                                                                                                              >  7:354042/1‑136 (MQ=255)
                     gCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTCGATATTgg                                                                                                                                              <  8:354042/136‑1 (MQ=255)
                      cAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCAt                                                                                                                                <  6:46514/149‑1 (MQ=255)
                      cAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCAt                                                                                                                                >  4:9439/1‑149 (MQ=255)
                       aaaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACtttt                                                                                                                                                       >  7:179604/1‑125 (MQ=255)
                       aaaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACtttt                                                                                                                                                       <  8:179604/125‑1 (MQ=255)
                       aaaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATAACAAAACTTTTTGATATTGGCCTTGCTAGGGCATc                                                                                                                               >  7:64887/1‑149 (MQ=255)
                          cATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGAt                                                                                                                            <  1:28594/149‑1 (MQ=255)
                           aTTTAAGCGCCAAGCGCCTGCTTTAAGCCGCAACTGAACAAGCAGAGCGAGACTCGCGGTGCCCCAAAGCGTGTAACATGAGGCTTTCGGCTTGCAGCTCATCGCCAATAACAAAACATATTGATATTGGCCTTGCTAGGGCATCGAta                                                                                                                           >  7:90011/1‑148 (MQ=255)
                            tttAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGc                                                                                                                          <  5:15886/149‑1 (MQ=255)
                            tttAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGc                                                                                                                          <  6:412886/149‑1 (MQ=255)
                               aagcgccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCt                                                                                                                        <  2:218193/148‑1 (MQ=255)
                               aagcgccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTg                                                                                                                       >  5:66449/1‑149 (MQ=255)
                                   gccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTc                                                                                                                     >  5:167269/1‑147 (MQ=255)
                                   gccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTa                                                                                                                  <  8:180968/149‑1 (MQ=255)
                                          gccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTgcc                                                                                                          >  5:26090/1‑149 (MQ=255)
                                                tttAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCa                                                                                                  >  2:263812/1‑149 (MQ=255)
                                                  tAAGCAGCAACTGAAGAAGCAGAGCGCGACTCGGTGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAATCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCa                                                                                                  <  4:223086/147‑1 (MQ=255)
                                                      ccGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTc                                                                                                                     >  8:361852/1‑128 (MQ=255)
                                                      ccGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTc                                                                                                                     <  7:361852/128‑1 (MQ=255)
                                                      ccGCAACTGAAGAAGCAGACCGCGACTCGGGCCGACCCAAAGCGAGTAACTTGAGCCTTTCGGCTTGCAGCTCATCGCCAAAATCAAAACTTTTAGAAAATTGCCTACCAAGGGAATCCATGCCTGTCTAGACTGCCC‑‑CCCCATGCAt                                                                                             >  4:285654/1‑148 (MQ=255)
                                                            cTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCtgcgtg                                                                                      >  4:144232/1‑149 (MQ=255)
                                                             tGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCtgcgtgc                                                                                     <  6:66449/149‑1 (MQ=255)
                                                               aagaagCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCtgcgtgc                                                                                     <  1:263812/147‑1 (MQ=255)
                                                                 gaagCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCaa                                                                                 <  6:167269/149‑1 (MQ=255)
                                                                       gagCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATc                                                                           >  4:408776/1‑149 (MQ=255)
                                                                       gagCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTATTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAAACTGCGTTCGCAATCGAt                                                                            >  8:320338/1‑148 (MQ=255)
                                                                          cgcgACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATcccc                                                                        >  2:338362/1‑149 (MQ=255)
                                                                            cgACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGa                                                                      <  8:149947/149‑1 (MQ=255)
                                                                            cgACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGa                                                                      <  3:220480/149‑1 (MQ=255)
                                                                              aCTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAg                                                                     <  4:197062/148‑1 (MQ=255)
                                                                              aCTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAg                                                                     >  3:197062/1‑148 (MQ=255)
                                                                               cTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGc                                                                    >  3:144731/1‑148 (MQ=255)
                                                                                  ggggCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGCTGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAcg                                                                >  5:207325/1‑149 (MQ=255)
                                                                                    ggCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAg                                                                     >  4:373335/1‑142 (MQ=255)
                                                                                    ggCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAg                                                                     <  3:373335/142‑1 (MQ=255)
                                                                                       gcCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGaaa                                                           <  3:144232/149‑1 (MQ=255)
                                                                                         cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATg                                                                                                <  6:265587/110‑1 (MQ=255)
                                                                                         cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCt                                                                                           >  1:246659/1‑115 (MQ=255)
                                                                                         cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCt                                                                                           <  2:246659/115‑1 (MQ=255)
                                                                                         cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTGTTTGATATTTGCCTTTCTATGGCCTAGATTCCTTTCTAGCCTTGCC‑‑CCCCATGAATATTCGTTCGCAATAGATACCCCAGCAACGCGAAACa                                                         >  1:318307/1‑149 (MQ=255)
                                                                                                  tgtAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGcc                                                >  7:247166/1‑149 (MQ=255)
                                                                                                    tAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCt                                              <  4:144731/149‑1 (MQ=255)
                                                                                                       cTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgtg                                           <  3:408776/149‑1 (MQ=255)
                                                                                                          gAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCa                                         >  4:233635/1‑148 (MQ=255)
                                                                                                           aGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCATCCGGGGCACGGACTCGATNAAACTGCGTGCGCAATAGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgtg                                           >  1:337946/1‑149 (MQ=255)
                                                                                                                  cGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGt                                 <  6:26090/148‑1 (MQ=255)
                                                                                                                         cAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgc                         <  3:233635/149‑1 (MQ=255)
                                                                                                                           gCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgc                       >  3:327264/1‑149 (MQ=255)
                                                                                                                            cTCATCTCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgct                      <  4:2070/149‑1 (MQ=255)
                                                                                                                               aTCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgctgcg                   >  3:156449/1‑149 (MQ=255)
                                                                                                                                 ggCCAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGcc                 >  3:171007/2‑149 (MQ=255)
                                                                                                                                  gCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgt                                            <  6:250215/121‑1 (MQ=255)
                                                                                                                                  gCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgt                                            >  5:250215/1‑121 (MQ=255)
                                                                                                                                  gCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCt                <  1:338362/149‑1 (MQ=255)
                                                                                                                                   ccAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCt                >  2:13350/1‑148 (MQ=255)
                                                                                                                                   ccAATATCAAAAAGTTTTGATATTGGCCATGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCAGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGc                                     >  1:232788/1‑127 (MQ=255)
                                                                                                                                   acAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATTCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGc                                     <  2:232788/126‑1 (MQ=255)
                                                                                                                                     aaTATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGcc                 <  6:309325/145‑1 (MQ=255)
                                                                                                                                     aaTATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGcc                 >  5:309325/1‑145 (MQ=255)
                                                                                                                                      atatCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTgg                                                     <  2:243103/108‑1 (MQ=255)
                                                                                                                                      atatCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTgg                                                     >  1:243103/1‑108 (MQ=255)
                                                                                                                                       tatCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTg               >  2:30345/1‑145 (MQ=255)
                                                                                                                                       tatCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTg               <  1:30345/145‑1 (MQ=255)
                                                                                                                                         tCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTa         >  2:216384/1‑149 (MQ=255)
                                                                                                                                          cAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTg                                <  2:56949/125‑1 (MQ=255)
                                                                                                                                          cAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTg                                >  1:56949/1‑125 (MQ=255)
                                                                                                                                           aaaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGTACGTAtt       <  8:64887/149‑1 (MQ=255)
                                                                                                                                           aaaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTAtt       >  6:346351/1‑149 (MQ=255)
                                                                                                                                            aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTg      >  2:361333/1‑149 (MQ=255)
                                                                                                                                            aaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGACTGAACGTATTg      >  7:182214/1‑149 (MQ=255)
                                                                                                                                             aaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTgc     <  1:216384/149‑1 (MQ=255)
                                                                                                                                                tttttGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTgcgcg  >  8:292323/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                        
CGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCG‑‑CCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTGCGCG  >  NC_002947/5657750‑5658042

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: