Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 840,650 T→A 5.8% I18F (ATC→TTC)  ispE ← 4‑diphosphocytidyl‑2C‑methyl‑D‑erythritol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947840,6500TA 5.8% 201.1 / 2.4 69I18F (ATC→TTC) ispE4‑diphosphocytidyl‑2C‑methyl‑D‑erythritol kinase
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base A (2/2);  total (32/37)
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.06e-01

AGGTTGCTGTCGTGGGGGACGGCCTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCAGCTTCAGGCGCTCGGGCAACC  >  NC_002947/840502‑840794
                                                                                                                                                    |                                                                                                                                                 
aGGTTGCTGTCGTGGGGGACGGCCTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCgatgat                                                                                                                                                   <  8:32192/149‑1 (MQ=255)
 ggTTGCTGTCGTGGGGGACGGCCTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCACTAAGTGGTAGCCGTCTGCGCGGCGGCCGatgatg                                                                                                                                                  >  4:259629/1‑149 (MQ=255)
      cTGTCGTGGGGGACGGCCTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCa                                                                                                                                              <  4:294580/148‑1 (MQ=255)
         tCGTGGGGGACGGCCTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGcc                                                                                                                                           <  3:259629/148‑1 (MQ=255)
         tCGTGGGGGACGGCCTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCca                                                                                                                                          >  2:215102/1‑149 (MQ=255)
                   cGGCCTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGc                                                                                                                                <  3:378546/149‑1 (MQ=255)
                       cTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTgg                                                                                                                            >  8:77475/1‑149 (MQ=255)
                       cTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTgg                                                                                                                            >  4:170634/1‑149 (MQ=255)
                       cTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTgg                                                                                                                            >  1:65106/1‑149 (MQ=255)
                       cTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGATTg                                                                                                                             >  1:348968/1‑148 (MQ=255)
                       cTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCAAAGGTTGAGCTTgg                                                                                                                            >  3:99354/1‑149 (MQ=255)
                                ggTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCggg                                                                                                                       <  2:65106/145‑1 (MQ=255)
                                  tgtgCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGAtgtg                                                                                                                                                <  4:306387/118‑1 (MQ=255)
                                  tgtgCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGAtgtg                                                                                                                                                >  3:306387/1‑118 (MQ=255)
                                     gCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGAAGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAgg                                                                                                              >  5:90805/1‑149 (MQ=255)
                                     gCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGAAGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAgg                                                                                                              >  5:96227/1‑149 (MQ=255)
                                      cAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAgg                                                                                                              <  2:149450/148‑1 (MQ=255)
                                      cAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAgg                                                                                                              <  3:170634/148‑1 (MQ=255)
                                       aaGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGt                                                                                                            >  6:324847/1‑149 (MQ=255)
                                          cGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAg                                                                                                         <  1:212039/149‑1 (MQ=255)
                                            gATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCt                                                                                                       >  8:29425/1‑149 (MQ=255)
                                             aTCACGCCATCGTCACGCCGGGCAAAACTCAGGTCGTCGCCATGGTCTAGACACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGc                                                                                                        <  6:202185/147‑1 (MQ=255)
                                                aCGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTg                                                                                                   >  3:352721/1‑149 (MQ=255)
                                                  gCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCa                                                                                                 <  5:324847/149‑1 (MQ=255)
                                                   ccATCGTCCCGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGGTGTAGCCGTCGGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCAt                                                                                                <  3:28534/149‑1 (MQ=255)
                                                     aTCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGc                                                                                              <  5:28532/149‑1 (MQ=255)
                                                        gTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGAt                                                                                           >  5:119642/1‑149 (MQ=255)
                                                          cACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATc                                                                                          >  8:79273/1‑148 (MQ=255)
                                                            cGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGAAGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAAt                                                                                       <  6:96227/149‑1 (MQ=255)
                                                            cGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGAAGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAAt                                                                                       <  6:90805/149‑1 (MQ=255)
                                                             gCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATg                                                                                      >  2:68440/1‑149 (MQ=255)
                                                                 gtcAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGATTCCAGTTCGTGGAAGCCGTCTGCGCGGCGGCCGATGATCTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAACTTTTGCATGCGATCAATGccc                                                                                   >  5:304896/3‑148 (MQ=255)
                                                                   cAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCGATGCCCCAg                                                                                >  2:297469/1‑149 (MQ=255)
                                                                      aaCTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTg                                                                             >  6:248690/1‑149 (MQ=255)
                                                                         tCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGgcg                                                                          >  8:233779/1‑149 (MQ=255)
                                                                                gtcgCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGc                                                                   <  2:189465/149‑1 (MQ=255)
                                                                                  cgCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCa                                                                 <  1:215102/149‑1 (MQ=255)
                                                                                   gCCATGGTCTAGAAACTGTAAGACCGTTTCCAGTTCGTGGTAGCCTTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGTGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCa                                                                 <  6:304896/148‑1 (MQ=255)
                                                                                      aTGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTcc                                                             <  7:79273/149‑1 (MQ=255)
                                                                                      aTGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCAACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTcc                                                             >  5:379576/1‑149 (MQ=255)
                                                                                          tCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGa                                                         <  4:372192/149‑1 (MQ=255)
                                                                                            tAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGac                                                        >  6:80990/1‑148 (MQ=255)
                                                                                               aaaCTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGaccac                                                     >  7:229262/1‑148 (MQ=255)
                                                                                                aaCTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGAccacca                                                   <  4:99354/149‑1 (MQ=255)
                                                                                                  cTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAg‑g                                                <  5:137473/149‑1 (MQ=255)
                                                                                                   tGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAg‑gg                                               >  4:30763/1‑149 (MQ=255)
                                                                                                   tGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAg‑gg                                               >  1:4659/1‑149 (MQ=255)
                                                                                                   tGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAg‑gg                                               >  1:60602/1‑149 (MQ=255)
                                                                                                         acCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGt                                         >  4:33541/1‑149 (MQ=255)
                                                                                                           cGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTc                                        >  3:195463/1‑148 (MQ=255)
                                                                                                            gtatCCCGTTCGTGGTAGCCGGCTGCGCGGCGGCCGATGATGTGCACCCACAGGGTCCGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCa                                      <  8:169772/146‑1 (MQ=255)
                                                                                                             tttCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAg                                     <  5:248690/149‑1 (MQ=255)
                                                                                                             tttCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAg                                     <  7:77475/149‑1 (MQ=255)
                                                                                                             tttCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGAATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAg                                     >  2:64559/1‑149 (MQ=255)
                                                                                                                 cAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTcc                                 <  5:80990/149‑1 (MQ=255)
                                                                                                                   gTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGc                                                                   >  7:68314/1‑114 (MQ=255)
                                                                                                                   gTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGc                                                                   <  8:68314/114‑1 (MQ=255)
                                                                                                                    ttCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTg                              <  6:119642/149‑1 (MQ=255)
                                                                                                                       gTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGc                                                                   <  3:370950/110‑1 (MQ=255)
                                                                                                                       gTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGc                                                                   >  4:370950/1‑110 (MQ=255)
                                                                                                                       gTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAg                                     <  2:115726/139‑1 (MQ=255)
                                                                                                                       gTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAg                                     >  1:115726/1‑139 (MQ=255)
                                                                                                                        tGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGt                          <  7:210127/149‑1 (MQ=255)
                                                                                                                        tGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGt                          >  7:330883/1‑149 (MQ=255)
                                                                                                                          gTAGCCGTCTGCGCGGCGGCCTATGATGGGCATCCACAGGTTGAGCTTGGTCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGc                        <  2:348968/149‑1 (MQ=255)
                                                                                                                          gTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGc                        <  2:4659/149‑1 (MQ=255)
                                                                                                                           tAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCa                       <  4:83374/149‑1 (MQ=255)
                                                                                                                              ccGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCAGCt                    >  6:76165/1‑149 (MQ=255)
                                                                                                                               cGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAGAGGTCACGTCCAGGTCCTTGCCGTGCAGCt                    <  1:226506/149‑1 (MQ=255)
                                                                                                                                gtctcggcggcgGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCAGCTTc                  <  2:328833/143‑1 (MQ=255)
                                                                                                                                     cgcgGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCAGCTTCAGgc              <  7:29425/148‑1 (MQ=255)
                                                                                                                                          cggcCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGCGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCAGCTTCAGGCGCTCgg        <  3:30763/149‑1 (MQ=255)
                                                                                                                                          cggcCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCAGCTTCAGGCGCTCgg        <  8:229262/149‑1 (MQ=255)
                                                                                                                                               gatgatGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCa                       <  1:85451/129‑1 (MQ=255)
                                                                                                                                               gatgatGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCa                       >  2:85451/1‑129 (MQ=255)
                                                                                                                                               gatgatGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCAGCTTCAGGCGCTCGGGCAAc   >  4:200779/1‑149 (MQ=255)
                                                                                                                                                atgatgTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCAGCTTCAGGCGCTCGGGCAAcc  <  4:352721/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGGTTGCTGTCGTGGGGGACGGCCTCGATCTCGGTGTGCAAGCGGATCACGCCATCGTCACGCAGGGCAAAACTCAGTTCGTCGCCATGGTCTAGAAACTGAAACACCGTTTCCAGTTCGTGGTAGCCGTCTGCGCGGCGGCCGATGATGTGCAGCCACAGGTTGAGCTTGGCCGGGGCGGGCAGGGTGAGCTTTTGCATGCGATCAATGCCCCAGCTGGCGCGGCTGCCAGTCCTTGACCACCAG‑GGTCACGTCCAGGTCCTTGCCGTGCAGCTTCAGGCGCTCGGGCAACC  >  NC_002947/840502‑840794

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: ‑