Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,054,906 C→T 19.1% intergenic (+46/‑55) quiA → / → oprB‑III quinate dehydrogenase/carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,054,9060CT19.1% 556.9 / 121.3 260intergenic (+46/‑55)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (89/120);  new base T (15/34);  total (104/156)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

Only 100 of 271 total aligned reads displayed.

ACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCC  >  NC_002947/4054757‑4055054
                                                                                                                                                     |                                                                                                                                                    
 cGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCt                                                                                                                                                     <  1:2512109/150‑1 (MQ=255)
  gTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTg                                                                                                                                                    <  2:3872430/150‑1 (MQ=255)
   tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGt                                                                                                                                                   <  2:3933901/150‑1 (MQ=255)
   tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGt                                                                                                                                                   <  2:4865766/150‑1 (MQ=255)
   tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGt                                                                                                                                                   <  1:1773622/150‑1 (MQ=255)
   tctcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGt                                                                                                                                                   <  1:5456478/150‑1 (MQ=255)
     tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAg                                                                                                                                                 <  2:2033121/150‑1 (MQ=255)
     tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                                 <  1:2297607/150‑1 (MQ=255)
     tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAg                                                                                                                                                 >  2:3453630/1‑150 (MQ=255)
       gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                               <  2:5945292/150‑1 (MQ=255)
       gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                               <  2:2115854/150‑1 (MQ=255)
       gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                               <  1:3253507/150‑1 (MQ=255)
       gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                               <  2:269296/150‑1 (MQ=255)
       gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACACGGCCCCTggagga                                                                                                                                               >  2:4271732/1‑150 (MQ=255)
       gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGa                                                                                                                                               <  2:6675436/150‑1 (MQ=255)
       gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGa                                                                                                                                               <  1:6484018/150‑1 (MQ=255)
         cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                             <  2:2919881/150‑1 (MQ=255)
         cGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                             <  2:320388/150‑1 (MQ=255)
          gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCg                                                                                                                                            <  2:2388931/150‑1 (MQ=255)
           aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                           <  2:2297413/150‑1 (MQ=255)
           aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCgg                                                                                                                                           <  1:2568085/150‑1 (MQ=255)
           aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCgg                                                                                                                                           <  2:5465415/150‑1 (MQ=255)
            aGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCggg                                                                                                                                          <  2:2414511/150‑1 (MQ=255)
            aGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGATTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCggg                                                                                                                                          <  2:4604802/150‑1 (MQ=255)
             gACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGc                                                                                                                                         <  1:3066556/150‑1 (MQ=255)
               ccGGCACGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTCCGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTCCAAGGCCCCTGTAGGAGCGGGCtt                                                                                                                                       <  2:110101/150‑1 (MQ=255)
               ccGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCtt                                                                                                                                       <  2:2440132/150‑1 (MQ=255)
                 ggCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGc                                                                                                                                     <  2:2856612/150‑1 (MQ=255)
                     aGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCcgc                                                                                                                                 <  2:4138293/150‑1 (MQ=255)
                           aCGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAACCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                           <  1:503372/150‑1 (MQ=255)
                                     cACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCCAGCCCCCGCCGACAAGGCCCCGGTAGGGGCGGGCTTGCCCGCGAATGCCCCACCAAGa                                                                                                                 >  1:3215989/1‑150 (MQ=255)
                                     cACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGa                                                                                                                 <  2:6282880/150‑1 (MQ=255)
                                      aCTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGGAATCCCCCTCCGACCAAGCCCCGGTAGGAGCCGGCGGGGCCCCGAATGCCCCCCCAAGAc                                                                                                                >  1:5169422/1‑150 (MQ=255)
                                        tgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACaa                                                                                                              >  2:4514996/1‑150 (MQ=255)
                                        tgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCGGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACaa                                                                                                              >  2:4913411/1‑150 (MQ=255)
                                        tgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCACCAACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACaa                                                                                                              >  1:3911329/1‑150 (MQ=255)
                                          ctgGCGGGCGCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCAACAAGGCCCCGGTAGGAGCGGGCGTGCCCGCGAAAGCCCCGCAAAGACAAcc                                                                                                            >  2:1146678/1‑150 (MQ=255)
                                              cggcgCCCGCCAGTCGACTGACCGTGGCGACTATGTTATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCgag                                                                                                        <  1:123244/150‑1 (MQ=255)
                                                   cccgccAGTCGACTGACCGTGGCGACTATGTTTTTTCTTACGCTCTGCCTTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGATCCTACAAGGCCCGTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCt                                                                                                   <  1:3636229/150‑1 (MQ=39)
                                                   cccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCt                                                                                                   <  2:1710802/150‑1 (MQ=255)
                                                   cccgccAGTCGACTGACCGTGGCGACTAGGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCGGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCt                                                                                                   <  2:5872798/150‑1 (MQ=255)
                                                        cAGTCGACTGACCGTGGCGACTATGTGATTGCTGACGCTCTGCCGTACAACCGCGTCGCCTGCTTCGCGGGCAACCACGCACCAACACGGCCCCAGGAGGAGCGGGCCTGCCCGCGAAACCCCCTCCAAGACCACCCCACAATCTGTaaa                                                                                              >  2:2245485/1‑149 (MQ=255)
                                                        cAGGCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAAt                                                                                              <  2:5564196/150‑1 (MQ=255)
                                                          gTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCGGGAGGAGCGGGCTTGCCCGCGAATGCCCCGCAAAGACAACCCGAGATTCTGTAATGc                                                                                            >  2:5001768/1‑150 (MQ=255)
                                                                  aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCa                                                                                    >  2:6204589/1‑150 (MQ=255)
                                                                  aCCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCTCCTACAAGGCCCCGGTAGGAGCGGGCTTGCCCGCGAATGCCCCGCAAAGACAACCCCAGATTCTGTAATGCCATCCGCa                                                                                    >  2:2469403/1‑150 (MQ=255)
                                                                      tGGCGACTATGTGATTTCTTACGCTCTGCCGTACAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTc                                                                                <  2:5205293/150‑1 (MQ=255)
                                                                      tGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTc                                                                                <  1:4412575/150‑1 (MQ=255)
                                                                        gCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCACGCCCGCTCCTACAAGGCCCCTGTTGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGc                                                                              <  1:6612217/150‑1 (MQ=255)
                                                                          gACTATGTGATTTCTTACGCTCTGCCTTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCAt                                                                            <  1:5196496/150‑1 (MQ=255)
                                                                          gACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCAt                                                                            <  2:3907589/150‑1 (MQ=255)
                                                                            cTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTCCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCa                                                                          <  1:409971/150‑1 (MQ=255)
                                                                             tctgtgATTTCTTCCGCTCTGCCGTAAACCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCAc                                                                         <  1:4913399/148‑1 (MQ=255)
                                                                             tATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCAACAAGGCCCCAGTAGGAGCGGGCGTGCCCGCGAATGCCCCACAAATACAACCCCAGCTTCTGTAATTCCATCCGCAAATCGCAacac                                                                         >  1:1113139/1‑150 (MQ=255)
                                                                                 tgATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCGCCTACAAGGCCCCGGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCAccccc                                                                     >  1:3895053/1‑150 (MQ=255)
                                                                                  gATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCGCCCACGCCCCCGGAGGCGCGGGCGTGCCCGCCAATGCCCCGCCAAGACAACCCCCGAATCAGTAATGCCCTCCGCCCTGCGCCGCAcccccc                                                                    >  2:5151710/1‑149 (MQ=255)
                                                                                        ttccGCTCTGCCGTAAACCCGCGTCGCCTGCGTCGCGGGCCAGCCCGCTCCTCCACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCctgc                                                              <  1:2936380/147‑1 (MQ=39)
                                                                                        ttACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCCCTACTACAAGGCCCCCGTAGGCGCCGGCTTGCCCGCGAATGCCCCACAAATACAAACCCACAATATGTAAATCCATCAGCAATTCGCATAACCCCCAcccacc                                                              >  2:3604452/1‑147 (MQ=255)
                                                                                              tctGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCACCTACAAGGCCCCTGTAGGAGCGGGCGTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCCCCCCCACCCTGCGGctgg                                                        >  2:4950414/1‑150 (MQ=255)
                                                                                               ctGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggc                                                       <  2:3683838/150‑1 (MQ=255)
                                                                                                  ccGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCACAAAAACAACCCGAGATTCTGTAAAGCCAACCGCAAATGGCCTCACCCCCACCCGGCAGctggccca                                                    >  2:4408568/1‑149 (MQ=37)
                                                                                                  ccGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCGACAAGGCCCCCGTTGGAGCCGGCTTGCCCGCGAATTCCCCACAAATACAACCCCAGATTCGGTAAAGCCCAACGCAAATCGCATCACCCCCACCCGGCTGctggcccg                                                    >  2:1983953/1‑149 (MQ=255)
                                                                                                   cGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctg                                                   <  1:6269095/150‑1 (MQ=255)
                                                                                                   cGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctg                                                   <  1:2114717/150‑1 (MQ=255)
                                                                                                    ggaaaaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCGTGCCCGCGGATGCCCCTCAAAGAAAACCCGAGATTCTGTAATGCCATCCGCCAGTCGCATAACCCCCAACCGGCTGctggccctgg                                                  >  1:5559922/3‑150 (MQ=38)
                                                                                                    gTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctgg                                                  <  2:2519653/150‑1 (MQ=255)
                                                                                                      aaaaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctggcc                                                <  2:1603195/150‑1 (MQ=255)
                                                                                                         aCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCGTGCCCGCGGATGCCCCGCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCGGCTGCTGGCCCGGGCCAGc                                             >  2:4776273/1‑150 (MQ=38)
                                                                                                           cgcgTCGCCTGCTTCGCGGGCAAGCCCGCACCTACACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcac                                           <  1:3562749/150‑1 (MQ=255)
                                                                                                             cgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACACGGCCCCTGGAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcacca                                         <  1:1692687/150‑1 (MQ=255)
                                                                                                                  ccTGCTTCGCGGGCAAGCCCGCTCCTACAACGCCCCGGGAGGAGCGGGCGTGCCCGCGGAAGCCCCGCCAAGACAACCCCCGATTCTGTAATGCCATCCGCAATTAGCATCACCCCCACCCAGCAGCTGGCCCGGGCCCGCACCCCTGcc                                    >  2:4636393/1‑150 (MQ=255)
                                                                                                                            gggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcg                          >  1:4832344/1‑150 (MQ=255)
                                                                                                                             ggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGCCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                         <  2:677038/150‑1 (MQ=255)
                                                                                                                              gCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgac                        >  2:4028022/1‑150 (MQ=255)
                                                                                                                                  gCCCGCTCCTACAAGGCCCCTGTAGGCGGGGGCTTGCCCGCGAATGCCCCTCAAACACAACCCGAGATTCTGTAATGCCATCCGCAATTCGAATCACCCCCACCCTGCTGCGGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtc                    <  1:4617067/150‑1 (MQ=255)
                                                                                                                                   cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct                   <  2:3375491/150‑1 (MQ=255)
                                                                                                                                    ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc                  >  1:3784664/1‑150 (MQ=255)
                                                                                                                                     cGCTCCTACCAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGAATCACCCCCACCCTACTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                 <  1:1485852/150‑1 (MQ=255)
                                                                                                                                     cGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                 >  2:5638227/1‑150 (MQ=255)
                                                                                                                                      gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                >  1:1894024/1‑150 (MQ=255)
                                                                                                                                      gCTCCTACAGGGGCCTGGTAGGAGCGGGCTTGCCCGCGAAGC‑‑‑‑‑‑‑AGGCGACGCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca        >  1:4772282/1‑150 (MQ=21)
                                                                                                                                      gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCACTTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                >  1:492863/1‑150 (MQ=255)
                                                                                                                                      gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                >  2:769213/1‑150 (MQ=255)
                                                                                                                                      gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                >  2:3845582/1‑150 (MQ=255)
                                                                                                                                        tCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa              >  1:4332664/1‑150 (MQ=255)
                                                                                                                                         ccTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc             <  2:4580625/150‑1 (MQ=255)
                                                                                                                                          cTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAcc            >  1:6019869/1‑150 (MQ=255)
                                                                                                                                          cTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAcc            >  2:6560492/1‑150 (MQ=255)
                                                                                                                                           tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc           >  1:6344966/1‑150 (MQ=255)
                                                                                                                                            aCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg          >  1:1155751/1‑150 (MQ=255)
                                                                                                                                            aCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg          >  1:6024142/1‑150 (MQ=255)
                                                                                                                                              atgcccccTGTTGGCGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca        <  1:806054/146‑1 (MQ=255)
                                                                                                                                              aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca        >  2:5522993/1‑150 (MQ=255)
                                                                                                                                               aGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAgcagcagg                                       >  2:3166116/1‑116 (MQ=255)
                                                                                                                                               aGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAgcagcagg                                       <  1:3166116/118‑3 (MQ=255)
                                                                                                                                                 gCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcagca     >  2:5355641/1‑150 (MQ=255)
                                                                                                                                                 gCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcagca     >  1:5670850/1‑150 (MQ=255)
                                                                                                                                                   cccTGTAGGCGCGGGCTTGCCCGCGAATGCCCCTCAAAGACACCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc   <  1:4271732/150‑1 (MQ=255)
                                                                                                                                                   cccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAACGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc   <  2:972209/150‑1 (MQ=255)
                                                                                                                                                   cccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc   >  1:1419797/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCC  >  NC_002947/4054757‑4055054

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