Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,460,072 T→C 5.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,460,0720TC 5.0% 448.9 / 7.8 140intergenic (+179/‑26)PP_2150/sthAthiamine biosynthesis lipoprotein/pyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base T (52/80);  new base C (2/5);  total (55/85)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 144 total aligned reads displayed.

GAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCAGGTCGGTG  >  NC_002947/2459923‑2460221
                                                                                                                                                     |                                                                                                                                                     
gAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCtt                                                                                                                                                       <  1:2344445/150‑1 (MQ=255)
  aTGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGc                                                                                                                                                     <  2:876905/150‑1 (MQ=255)
   tGGGTGAAGCATATGTAGTGCGGGCAACATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCg                                                                                                                                                    <  2:1161586/150‑1 (MQ=255)
    gggTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCAAGCCACGGTCGCTCGCgg                                                                                                                                                   <  1:901837/150‑1 (MQ=255)
    gggTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCAAGCCACGGTCGCTCGCgg                                                                                                                                                   <  1:901712/150‑1 (MQ=255)
      gTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGc                                                                                                                                                 >  1:3192105/1‑150 (MQ=255)
      gTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCAAGCCACGGTCGCTCGCgggg                                                                                                                                                 <  1:3433389/150‑2 (MQ=255)
       tgcaGCATATGTAGTGCGGGCAACAGTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGGTCAGGGGTGATGAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCAAGCCACGGTCGCTCGCggggc                                                                                                                                                <  2:3095513/147‑3 (MQ=255)
           gCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGtt                                                                                                                                            >  2:2632462/1‑150 (MQ=255)
            cATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTc                                                                                                                                           <  1:1222097/150‑1 (MQ=255)
                tGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCTGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCTTGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGaa                                                                                                                                       <  2:2224299/150‑1 (MQ=255)
                tGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGTCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGaa                                                                                                                                       <  2:2275072/150‑1 (MQ=255)
                tGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGaa                                                                                                                                       <  1:736515/150‑1 (MQ=255)
                 gTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAg                                                                                                                                      >  2:4333382/1‑150 (MQ=255)
                   aGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGa                                                                                                                                    <  1:3171083/150‑1 (MQ=255)
                   aGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGa                                                                                                                                    <  1:2467732/150‑1 (MQ=255)
                   aGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGa                                                                                                                                    <  2:907961/150‑1 (MQ=255)
                   aGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGa                                                                                                                                    <  2:1365647/150‑1 (MQ=255)
                   aGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGa                                                                                                                                    <  1:2061733/150‑1 (MQ=255)
                   aGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGa                                                                                                                                    <  1:180738/150‑1 (MQ=255)
                   aGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGa                                                                                                                                    <  1:1358836/150‑1 (MQ=255)
                      gCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTa                                                                                                                                 <  1:4415440/150‑1 (MQ=255)
                       cGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTAc                                                                                                                                <  1:865010/150‑1 (MQ=255)
                          gCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCa                                                                                                                             >  1:266900/1‑150 (MQ=255)
                          gCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCa                                                                                                                             >  1:368475/1‑150 (MQ=255)
                            aaaaTTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATg                                                                                                                           >  1:3760556/1‑150 (MQ=255)
                             aaaTTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGTGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATgg                                                                                                                          <  2:2129168/150‑1 (MQ=255)
                             aaaTTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGTGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATgg                                                                                                                          <  2:2129297/150‑1 (MQ=255)
                             aaaTTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCa                                                                                                                             <  2:1041861/147‑1 (MQ=255)
                             aaaTTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCa                                                                                                                             >  1:1041861/1‑147 (MQ=255)
                                ttGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTg                                                                                                                       <  1:761327/150‑1 (MQ=255)
                                  ggCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTc                                                                                                                     <  2:106993/150‑1 (MQ=255)
                                   gCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCt                                                                                                                    >  2:2040002/1‑150 (MQ=255)
                                    ccTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTa                                                                                                                   <  2:4956501/150‑1 (MQ=255)
                                      tACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACa                                                                                                                 >  2:733544/1‑150 (MQ=255)
                                      tACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACa                                                                                                                 >  1:2884848/1‑150 (MQ=255)
                                           cgcgACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTac                                                                                                            <  1:1998057/150‑1 (MQ=255)
                                              caacAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTacgac                                                                                                         <  2:1414151/147‑1 (MQ=255)
                                                ccAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGt                                                                                                       <  2:3480964/150‑1 (MQ=255)
                                                 cAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTa                                                                                                      >  1:3523642/1‑150 (MQ=255)
                                                 cAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTa                                                                                                      >  2:2747889/1‑150 (MQ=255)
                                                  aaGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAg                                                                                                     >  1:3469049/1‑150 (MQ=255)
                                                   aGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgt                                                                                                    >  1:3511110/1‑150 (MQ=255)
                                                     ggTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtgg                                                                                                  >  1:1711399/1‑150 (MQ=255)
                                                       ttAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtggtg                                                                                                >  2:1543116/1‑150 (MQ=255)
                                                       ttAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAgtggtg                                                                                                >  2:1543185/1‑150 (MQ=255)
                                                         aaTGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCt                                                                                              <  2:1110910/150‑1 (MQ=255)
                                                               cgcgGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTc                                                                                        >  1:2466796/1‑150 (MQ=255)
                                                                 cggcgTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccg                                                                                      <  2:4495237/150‑1 (MQ=255)
                                                                    cgTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccggcc                                                                                   <  2:2963482/150‑1 (MQ=255)
                                                                     gTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccggccc                                                                                  >  1:5114718/1‑150 (MQ=255)
                                                                     gTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTccggccc                                                                                  <  1:1722093/150‑1 (MQ=255)
                                                                           ggTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCccggccgg                                                                            <  2:368475/150‑1 (MQ=255)
                                                                            gTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTATGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGt                                                                           <  1:2894770/150‑1 (MQ=255)
                                                                            gTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGt                                                                           <  2:2201235/150‑1 (MQ=255)
                                                                                aTTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAg                                                                       <  1:2224274/150‑1 (MQ=255)
                                                                                 ttAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAgg                                                                      >  2:2565386/1‑150 (MQ=255)
                                                                                    gACTGCACCCGTATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACGACGACGTAGTGGTGCTGGGGTCCGGCCCGGCCGGTGAAGgcgc                                                                   <  1:2646608/150‑1 (MQ=255)
                                                                                       tGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCt                                                                                                                        <  2:3109848/94‑1 (MQ=255)
                                                                                       tGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGc                                                                <  1:566475/150‑1 (MQ=255)
                                                                                          aCCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAAt                                                             <  2:2089219/150‑1 (MQ=255)
                                                                                          aCCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAAt                                                             <  2:3435176/150‑1 (MQ=255)
                                                                                            ccGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGa                                                           >  2:4876136/1‑150 (MQ=255)
                                                                                            ccGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGa                                                           >  1:3374266/1‑150 (MQ=255)
                                                                                             cGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGaa                                                          >  2:1275362/1‑150 (MQ=255)
                                                                                             cGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGaa                                                          <  1:3958407/150‑1 (MQ=255)
                                                                                                   ttCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAcgccgc                                                    >  1:2194464/1‑150 (MQ=255)
                                                                                                   ttCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAAcgccgc                                                    <  2:709798/150‑1 (MQ=255)
                                                                                                      tcATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCaa                                                 >  2:4278230/1‑150 (MQ=255)
                                                                                                       cATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCaaa                                                <  1:605135/150‑1 (MQ=255)
                                                                                                        aTGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCACTGAACGCCGCCAAAg                                               >  1:3352059/1‑150 (MQ=255)
                                                                                                            cgcCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAgg                                           <  1:4437551/150‑1 (MQ=255)
                                                                                                            cgcCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAgg                                           <  1:4672958/150‑1 (MQ=255)
                                                                                                            cgcCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAgg                                           <  1:4673081/150‑1 (MQ=255)
                                                                                                              ccccGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGgc                                         <  1:1001167/150‑1 (MQ=255)
                                                                                                               cccGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGgcg                                        >  1:2455014/1‑150 (MQ=255)
                                                                                                                ccGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGgcgc                                       >  2:3332085/1‑150 (MQ=255)
                                                                                                                 cggcggCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCa                                      >  2:4491795/1‑150 (MQ=255)
                                                                                                                  gggggcATGATCGGGCCCCGCGAGCGACCTTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCaa                                     <  1:1206450/147‑1 (MQ=255)
                                                                                                                     ggcATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGt                                  >  2:1089404/1‑150 (MQ=255)
                                                                                                                      gcATGATCGTGCCCCGCGAGCGACCGTGACTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTa                                 <  1:848749/150‑1 (MQ=255)
                                                                                                                        aagATCGAGCCCCGCGAGCGACCGTGGCTTGCGGCCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTCGc                               >  2:2306874/3‑150 (MQ=255)
                                                                                                                          gATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGa                             >  1:3674395/1‑150 (MQ=255)
                                                                                                                            tCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATg                           >  1:113557/1‑150 (MQ=255)
                                                                                                                             cGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATgg                          <  2:4300657/150‑1 (MQ=255)
                                                                                                                               tGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTc                        <  2:4308221/150‑1 (MQ=255)
                                                                                                                                 ccccGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGa                      <  1:996667/150‑1 (MQ=255)
                                                                                                                                  cccGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGAt                     >  2:3649116/1‑150 (MQ=255)
                                                                                                                                  cccGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGAt                     <  2:3192105/150‑1 (MQ=255)
                                                                                                                                   ccGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCaa                                     >  2:4617289/1‑133 (MQ=255)
                                                                                                                                    cgcgAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGa                   >  2:3459415/1‑150 (MQ=255)
                                                                                                                                    cgcgAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGa                   <  1:3302955/150‑1 (MQ=255)
                                                                                                                                      cgagcgaCCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGAcc                 <  1:4084958/150‑1 (MQ=255)
                                                                                                                                      cgagcgaCCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGAcc                 <  2:2304640/150‑1 (MQ=255)
                                                                                                                                         gcgaCCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTc              <  2:288938/150‑1 (MQ=255)
                                                                                                                                           gaCCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGACGGTCATCGACCATCGc            <  1:1528966/150‑1 (MQ=255)
                                                                                                                                              cGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCTGGGCGCAAGGTAGCGACGGTCATCGACCATCGctac         <  2:2542062/150‑4 (MQ=255)
                                                                                                                                                   cTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCAGGTCgg    <  2:266900/150‑1 (MQ=255)
                                                                                                                                                    ttGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCAGGTCGGt   >  1:1190210/1‑150 (MQ=255)
                                                                                                                                                     tGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCAGGTCGGTg  <  2:1163247/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAATGGGTGAAGCATATGTAGTGCGGGCAAAATTGGCCTACGACGCGACCAAGGGTTAATGTGCGCGGCGTTCAGGGTTGATTAGACTGCACCCGAATTTTCTCATGGCGCCCCGGCGGCATGATCGTGCCCCGCGAGCGACCGTGGCTTGCGGGCCAGTTCTGAAGGAGTACGCATGGCTGTCTACAACTACGACGTAGTGGTGCTGGGTTCCGGCCCGGCCGGTGAAGGCGCGGCAATGAACGCCGCCAAAGCAGGGCGCAAGGTAGCGATGGTCGATGACCGTCGCCAGGTCGGTG  >  NC_002947/2459923‑2460221

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: