Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,064,161 G→A 12.2% L42L (CTG→CTA) ‡ yccE → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,064,1610GA12.2% 131.8 / 10.0 49L42L (CTG→CTA) ‡yccEuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (25/18);  new base A (2/4);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

TTAAATTTTATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAATA  >  NC_000913/1064014‑1064299
                                                                                                                                                   |                                                                                                                                          
ttaAATTTTATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTgg                                                                                                                                           >  4:473522/1‑149 (MQ=255)
ttaAATTTTATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCgg                                                                                                                                           >  3:67324/1‑149 (MQ=255)
     ttttATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtgtgt                                                                                                                                      >  2:568537/1‑149 (MQ=255)
       ttATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttt                                                                                                                                    <  2:594956/149‑1 (MQ=255)
       ttATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttt                                                                                                                                    <  3:137643/149‑1 (MQ=255)
             aCGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCg                                                                                                                              <  3:446601/149‑1 (MQ=255)
                 tatataTGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTta                                                                                                                          >  2:411149/1‑149 (MQ=255)
                      atGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACaa                                                                                                                     <  1:423055/149‑1 (MQ=255)
                       tGGGTAGTAATCTACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAAt                                                                                                                    >  1:514013/1‑149 (MQ=255)
                        gggTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATg                                                                                                                   >  1:687565/1‑149 (MQ=255)
                          gtagtaATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGaa                                                                                                                 <  1:568537/149‑1 (MQ=255)
                           tagtaATATACATGGAATTAGTTGCAATGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAAc                                                                                                                <  4:67324/149‑1 (MQ=255)
                            agtaATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACt                                                                                                               >  2:182988/1‑149 (MQ=255)
                              taATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCa                                                                                                             <  3:473522/149‑1 (MQ=255)
                                atatACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGa                                                                                                           <  1:411149/149‑1 (MQ=255)
                                      aTGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaa                                                                                                     >  2:344827/1‑149 (MQ=255)
                                       tGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaag                                                                                                    >  4:321155/1‑149 (MQ=255)
                                       tGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaag                                                                                                    >  4:683185/1‑149 (MQ=255)
                                        ggAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagt                                                                                                   >  3:546741/1‑149 (MQ=255)
                                         gAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagta                                                                                                  >  2:276316/1‑149 (MQ=255)
                                          aaTTAGTTGCACTGGAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATAAACTCAGACgtaaagtaa                                                                                                 >  2:369548/1‑149 (MQ=255)
                                          aaTTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTtat                                                                                                                         <  3:389417/125‑1 (MQ=255)
                                          aaTTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTtat                                                                                                                         >  4:389417/1‑125 (MQ=255)
                                            ttAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAAt                                                                                               >  4:7412/1‑149 (MQ=255)
                                                 tGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAaga                                                                                          <  1:344827/149‑1 (MQ=255)
                                                       gCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATAc                                                                                    <  2:514013/149‑1 (MQ=255)
                                                          aataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTcc                                                                                 >  1:651637/1‑149 (MQ=255)
                                                            taatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                                                             <  1:146292/103‑1 (MQ=255)
                                                            taatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                                                             >  2:146292/1‑103 (MQ=255)
                                                               ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGc                                                                                                                               <  1:156935/98‑1 (MQ=255)
                                                               ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGc                                                                                                                               >  2:156935/1‑98 (MQ=255)
                                                                tattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCt                                                                           <  1:276316/149‑1 (MQ=255)
                                                                     gAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTcc                                                                                 >  1:106243/1‑138 (MQ=255)
                                                                     gAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTcc                                                                                 <  2:106243/138‑1 (MQ=255)
                                                                       aaCAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                                                                    >  4:30285/1‑149 (MQ=255)
                                                                                            aaaaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAAttt                                               >  3:511209/1‑149 (MQ=255)
                                                                                                                 aaaCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataat                          >  4:450191/1‑149 (MQ=255)
                                                                                                                   aCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataatga                        >  3:57468/1‑149 (MQ=255)
                                                                                                                   aCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataatga                        >  4:1261/1‑149 (MQ=255)
                                                                                                                   aCATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGaatgataatga                        >  3:28071/1‑149 (MQ=255)
                                                                                                                     atatTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATc                      <  4:546741/149‑1 (MQ=255)
                                                                                                                       atTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCATACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTa                    >  4:632568/1‑149 (MQ=255)
                                                                                                                           aaTCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCt                                                                           >  1:494845/1‑90 (MQ=255)
                                                                                                                           aaTCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCt                                                                           <  2:494845/90‑1 (MQ=255)
                                                                                                                                 gCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCAttattta          <  3:683185/149‑1 (MQ=255)
                                                                                                                                    ttttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATga       <  4:28071/149‑1 (MQ=255)
                                                                                                                                    ttttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATga       <  2:651637/149‑1 (MQ=255)
                                                                                                                                          aaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAata  >  3:232801/1‑148 (MQ=255)
                                                                                                                                          aaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAata  <  4:232801/148‑1 (MQ=255)
                                                                                                                                          aaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAata  >  1:51468/1‑148 (MQ=255)
                                                                                                                                          aaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAata  <  2:51468/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
TTAAATTTTATAGACGCTATATATGGGTAGTAATATACATGGAATTAGTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGGAAATGCAGAATGATAATGATCTATCATTATTTATGAGAATA  >  NC_000913/1064014‑1064299

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