Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,458,107 T→C 6.0% intergenic (+179/‑26) PP_2150 → / → sthA thiamine biosynthesis lipoprotein/pyridine nucleotide transhydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,458,1070TC 6.0% 198.6 / 2.4 67intergenic (+179/‑26)PP_2150/sthAthiamine biosynthesis lipoprotein/pyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (2/2);  total (32/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 = 2.51e-01

GGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCAGGTC  >  NC_002947/2457963‑2458252
                                                                                                                                                |                                                                                                                                                 
ggTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCgg                                                                                                                                               >  4:330523/1‑149 (MQ=255)
ggTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCg                                                                                                                                                <  2:110706/148‑1 (MQ=255)
 gTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCAAGCCACGGTCGCTCGCggg                                                                                                                                              <  8:178471/149‑1 (MQ=255)
  tGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGc                                                                                                                                             <  7:303846/149‑1 (MQ=255)
  tGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGc                                                                                                                                             <  5:227090/149‑1 (MQ=255)
      gCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGt                                                                                                                                         >  1:211445/1‑149 (MQ=255)
             tAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCgg                                                                                                                                               >  6:283034/1‑136 (MQ=255)
             tAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCgg                                                                                                                                               <  5:283034/136‑1 (MQ=255)
                    ggCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACg                                                                                                                           >  1:117517/1‑149 (MQ=255)
                         aaTTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATgg                                                                                                                      >  1:154118/1‑149 (MQ=255)
                             ggCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTCCTGAAGGAGTACGCATGGCTg                                                                                                                   <  6:258109/148‑1 (MQ=255)
                              gCCTACGACGCTACCAAGTGGTAATGTGCGAGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGGCCCGGCGGCATGATCGTGCCCCACGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTc                                                                                                                 <  8:351793/149‑1 (MQ=255)
                               ccTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTc                                                                                                                 >  1:115124/1‑148 (MQ=255)
                               ccTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCt                                                                                                                <  3:16166/149‑1 (MQ=255)
                                 tACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATCATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTa                                                                                                               >  2:191271/1‑148 (MQ=255)
                                    gacgCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAAc                                                                                                           >  7:70920/1‑149 (MQ=255)
                                       gcgACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTg                                                                                                                   >  2:169524/1‑138 (MQ=255)
                                       gcgACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTg                                                                                                                   <  1:169524/138‑1 (MQ=255)
                                         gACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTGTACAACTacga                                                                                                      >  4:201838/1‑149 (MQ=255)
                                                ggTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtg                                                                                               >  5:312885/1‑149 (MQ=255)
                                                ggTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtg                                                                                               <  1:217933/149‑1 (MQ=255)
                                                  ttAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtggt                                                                                             >  1:44084/1‑149 (MQ=255)
                                                     aTGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGc                                                                                           <  7:33448/148‑1 (MQ=255)
                                                     aTGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCt                                                                                          <  5:180344/149‑1 (MQ=255)
                                                            cggcgTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTcc                                                                                   <  2:154118/149‑1 (MQ=255)
                                                             ggcgTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccg                                                                                  >  1:221194/1‑149 (MQ=255)
                                                                 ttCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccggccc                                                                              <  2:226633/149‑1 (MQ=255)
                                                                          gATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTg                                                                      >  6:235065/1‑148 (MQ=255)
                                                                                 cTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGcgcg                                                              <  7:311520/149‑1 (MQ=255)
                                                                                    cACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCgg                                                             <  2:211445/147‑1 (MQ=255)
                                                                                     aCCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccggc                                                                                <  5:248977/127‑1 (MQ=255)
                                                                                     aCCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccggc                                                                                >  6:248977/1‑127 (MQ=255)
                                                                                       ccGATTTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAAt                                                         >  3:333958/1‑148 (MQ=255)
                                                                                       ccGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAAt                                                         >  4:38661/1‑148 (MQ=255)
                                                                                         gAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGaa                                                      >  5:329189/1‑149 (MQ=255)
                                                                                             tttCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAcgcc                                                  >  4:301070/1‑149 (MQ=255)
                                                                                             acgatcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAc                                                     <  5:282645/142‑1 (MQ=255)
                                                                                             acgatcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAc                                                     >  6:282645/5‑146 (MQ=255)
                                                                                               gatcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAc                                                     >  6:1671/3‑144 (MQ=255)
                                                                                               gatcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAc                                                     <  5:1671/142‑1 (MQ=255)
                                                                                                ctcATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGTAGTACGCATGGCTGTCTACAACTACGACTTAGTGGTGCTTGGTTCCGGCCCGGCCGGTGAATGCGCGGCAATGAAcgccgc                                                <  8:70920/148‑1 (MQ=255)
                                                                                                ctcATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAcgccgc                                                <  3:301070/148‑1 (MQ=255)
                                                                                                atcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAcgccgc                                                >  6:132197/2‑148 (MQ=255)
                                                                                                 tcATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACgccgcc                                               >  4:308684/1‑148 (MQ=255)
                                                                                                 tcATGGCGCCCAGGCGGCATGATCGTGCCCCGCGAGCGACCGAGGCTTGCGGGCAAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGACAGGCCGTTGAAGGCGCGGCAATGAACgccgcc                                               >  6:20081/1‑148 (MQ=255)
                                                                                                      gcgcCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCa                                         <  2:117517/149‑1 (MQ=255)
                                                                                                        gcCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAgg                                       <  2:115124/149‑1 (MQ=255)
                                                                                                        gcCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAgg                                       >  2:344400/1‑149 (MQ=255)
                                                                                                        gcCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAgg                                       <  7:241150/149‑1 (MQ=255)
                                                                                                         ccccGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAggg                                      <  2:221194/149‑1 (MQ=255)
                                                                                                          cccGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGgc                                     <  8:105662/149‑1 (MQ=255)
                                                                                                          cccGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGgc                                     <  3:330523/149‑1 (MQ=255)
                                                                                                           ccGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGgcg                                    >  4:181592/1‑149 (MQ=255)
                                                                                                                      gaactggcccGCAAGCCACGGTCGCTCGCGGGGCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGa                         >  1:44805/8‑149 (MQ=38)
                                                                                                                        cGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCAACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATg                       >  2:220062/1‑149 (MQ=255)
                                                                                                                         gTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTAGGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATg                       <  1:191271/148‑1 (MQ=255)
                                                                                                                         gTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATgg                      >  5:17141/1‑149 (MQ=255)
                                                                                                                           gCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTc                    >  2:189189/1‑149 (MQ=255)
                                                                                                                            gcccGCAAGCCACGGTCGCTCGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGc                                          <  7:252475/125‑1 (MQ=255)
                                                                                                                            gcccGCAAGCCACGGTCGCTCGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGc                                          >  8:252475/2‑126 (MQ=255)
                                                                                                                             cccGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGa                  <  5:132197/149‑1 (MQ=255)
                                                                                                                              ccGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGAt                 <  6:329189/149‑1 (MQ=255)
                                                                                                                              ccGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGAt                 <  3:38661/149‑1 (MQ=255)
                                                                                                                                       aCCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGc        <  3:201838/149‑1 (MQ=255)
                                                                                                                                       aCCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGACGGTCATCGACCATCGc        <  6:153440/149‑1 (MQ=255)
                                                                                                                                         cGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGACGGTCATCGACCATCGcta      <  6:312885/149‑3 (MQ=255)
                                                                                                                                             gCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCAGGTc  >  2:322969/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCAGGTC  >  NC_002947/2457963‑2458252

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: