Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,458,100 C→G 6.8% intergenic (+172/‑33) 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,1000CG 6.8% 229.9 / 6.1 73intergenic (+172/‑33)PP_2150/sthAthiamine biosynthesis lipoprotein/pyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base C (33/35);  new base G (2/3);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

CCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCA  >  NC_002947/2457953‑2458248
                                                                                                                                                   |                                                                                                                                                    
cgaGGGTATGTGTGAATGTTATGTTGTGCGGGCATATTTGGCGTACGACTCGACCAAGGTTTAATGTGCGCGGCGTTCATGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCTAGCGACCg                                                                                                                                                     <  6:13737/147‑1 (MQ=255)
ccAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCg                                                                                                                                                     >  5:153440/1‑149 (MQ=255)
ccAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCg                                                                                                                                                     >  1:226633/1‑149 (MQ=255)
  aGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTg                                                                                                                                                   >  8:311520/1‑149 (MQ=255)
  aGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTg                                                                                                                                                   <  4:51345/149‑1 (MQ=255)
      aaTGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCAAGCCACGGTCGCtc                                                                                                                                               <  7:315047/149‑2 (MQ=255)
          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)
                gtatatGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGt                                                                                                                                                    <  6:185829/132‑1 (MQ=255)
                gtatatGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGt                                                                                                                                                    >  5:185829/3‑134 (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                                                                                           <  1:217933/149‑1 (MQ=255)
                                                          ggTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtg                                                                                           >  5:312885/1‑149 (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                                                                            >  6:248977/1‑127 (MQ=255)
                                                                                               aCCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccggc                                                                            <  5:248977/127‑1 (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                                                 >  6:282645/5‑146 (MQ=255)
                                                                                                       acgatcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAc                                                 <  5:282645/142‑1 (MQ=255)
                                                                                                         gatcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAc                                                 <  5:1671/142‑1 (MQ=255)
                                                                                                         gatcATGCCGCCGGGGCGCCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAc                                                 >  6:1671/3‑144 (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             <  3:38661/149‑1 (MQ=255)
                                                                                                                                        ccGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGAt             <  6:329189/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)
                                                                                                                                                   |                                                                                                                                                    
CCAGGGAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCA  >  NC_002947/2457953‑2458248

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: