Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,901,848 (C)5→6 coding (1257/1278 nt) yqcE → putative transport protein YqcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,901,8431.C95.5% 147.2 / ‑0.1 44coding (1252/1278 nt)yqcEputative transport protein YqcE
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (19/23);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

CCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAGTGGTGATGTACGATCTGGA  >  NC_000913/2901704‑2901966
                                                                                                                                            |                                                                                                                           
cccAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAgcg                                                                                                                              >  1:75800/1‑140 (MQ=255)
     ttATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTAGCTGTTCAAAAGAAGTCGTACGGCTGATAGCGccccc                                                                                                                        >  1:267496/1‑140 (MQ=255)
         cGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGccccccgca                                                                                                                    <  1:538619/140‑5 (MQ=255)
               cTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTCCCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCt                                                                                                               >  2:808674/1‑140 (MQ=255)
               cTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCt                                                                                                               <  1:118316/140‑1 (MQ=255)
               cTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCt                                                                                                               >  2:118316/1‑140 (MQ=255)
               cTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCt                                                                                                               <  1:740995/140‑1 (MQ=255)
               cTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGccccccgca                                                                                                                    >  2:40335/1‑130 (MQ=255)
               cTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGccccccgca                                                                                                                    <  1:40335/134‑5 (MQ=255)
                   cAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTagca                                                                                                           <  2:812025/140‑1 (MQ=255)
                     aGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTagcagc                                                                                                         <  1:590004/140‑1 (MQ=255)
                          cccTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCGCAATGGCTAGCAGCAAGTaa                                                                                                   >  1:810922/1‑140 (MQ=255)
                            cTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAg                                                                                                  >  2:135549/1‑140 (MQ=255)
                             tGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCGCAATGGCTAGCAGCAAGTAAGcc                                                                                                <  2:810922/140‑1 (MQ=255)
                                   cTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGaa                                                                                           <  1:54145/140‑1 (MQ=255)
                                        gAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTa                                                                                      <  1:135549/140‑1 (MQ=255)
                                              tACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTCCCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGgag                                                                                <  1:808674/140‑1 (MQ=255)
                                                            gCTGATGGGCATGGCGGCTCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGa                                                                  <  1:213232/140‑1 (MQ=255)
                                                             cTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGaa                                                                 >  1:34311/1‑140 (MQ=255)
                                                              tGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGaaa                                                                >  1:442372/1‑140 (MQ=255)
                                                                aTGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTCGCAGCAAGTAAGGCAGACAAATTACCGGAGTTAAATCATTTCGaaaaa                                                              >  2:402855/1‑140 (MQ=255)
                                                                aTGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGaaaaa                                                              >  2:713333/1‑140 (MQ=255)
                                                                  gggCTTGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAAt                                                            >  1:116427/1‑140 (MQ=255)
                                                                  gggCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAAt                                                            >  2:427375/1‑140 (MQ=255)
                                                                   ggCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATa                                                           <  2:580271/140‑1 (MQ=255)
                                                                    gCATGGCGGCGCTTGACATGGTCATCGTCGTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATAc                                                          <  1:402855/140‑1 (MQ=255)
                                                                        ggcggcGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATAcatca                                                      <  2:31851/140‑1 (MQ=255)
                                                                        ggcggcGCTTGGCATGGTGATTCTCTTTACCTTTTCGCTATTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATAcatca                                                      <  2:267496/140‑1 (MQ=255)
                                                                                  ggCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGa                                            <  2:741447/140‑1 (MQ=255)
                                                                                     aTGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGAAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATgg                                         <  1:580420/140‑1 (MQ=255)
                                                                                            ttGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGt                                  >  2:811281/1‑140 (MQ=255)
                                                                                             tGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTc                                 >  1:606102/1‑140 (MQ=255)
                                                                                                 tttACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCagag                             <  1:427375/140‑1 (MQ=255)
                                                                                                   tACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGca                           >  1:106588/1‑140 (MQ=255)
                                                                                                      cTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACaa                        <  1:638840/140‑1 (MQ=255)
                                                                                                         tttGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTGAGAGCACAAAAg                     >  2:544360/1‑140 (MQ=255)
                                                                                                         tttGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAg                     <  2:606102/140‑1 (MQ=255)
                                                                                                         tttGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAg                     >  1:560702/1‑140 (MQ=255)
                                                                                                         tttGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATAAATCATAGGGATGGAGGGCGGAAGTCAGAGCACAAAAg                     >  1:232546/1‑140 (MQ=255)
                                                                                                         tttGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATAAATCATAGGGATCGAAGGCGGAAGTCAGAGCACAAAAg                     >  1:416863/1‑140 (MQ=255)
                                                                                                             cTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCa                                >  2:52547/1‑125 (MQ=255)
                                                                                                             cTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCa                                <  1:52547/125‑1 (MQ=255)
                                                                                                             cTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAgtggt                 <  1:713333/140‑1 (MQ=255)
                                                                                                                     aaaaaTTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAGTGGTGATGTACg         <  1:811281/140‑1 (MQ=255)
                                                                                                                     aaaaaTTCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAGTGGTGATGTACg         <  2:560702/140‑1 (MQ=255)
                                                                                                                          ttCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAGTGGTGATGTACGATCTg    >  1:115389/1‑140 (MQ=255)
                                                                                                                          ttCGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGAGGGCAGAAGTCAGAGCAGAAAAATGGTGATGTACGATCTg    >  2:308811/1‑140 (MQ=255)
                                                                                                                            cGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAGTGGTGATGTACGATCTGGa  <  2:106588/140‑1 (MQ=255)
                                                                                                                            cGTACTGCTGATAGCGCCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAGTGGTGATGTACGATCTGGa  <  2:75800/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                           
CCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCATGGCGGCGCTTGGCATGGTGATTGTCTTTACCTTTTTGCTGTTCCAAAAAATTCGTACTGCTGATAGCGCCCCCGCAATGGCTAGCAGCAAGTAAGCCAGGAAAATTACAGGAGAAAAATCATGTCGAAAAAATACATCATAGGGATTGATGGCGGAAGTCAGAGCACAAAAGTGGTGATGTACGATCTGGA  >  NC_000913/2901704‑2901966

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