Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,329,419 G→A Q127* (CAA→TAA)  yfaQ ← tandem DUF2300 domain‑containing protein YfaQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,329,4190GA97.9% 153.9 / ‑4.3 48Q127* (CAA→TAA) yfaQtandem DUF2300 domain‑containing protein YfaQ
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (27/20);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

AAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTC  >  NC_000913/2329283‑2329547
                                                                                                                                        |                                                                                                                                
aaGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTAccc                                                                                                                               >  1:273077/1‑140 (MQ=255)
     ttCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCcaca                                                                                                                          >  2:407406/1‑140 (MQ=255)
       ccTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCcacaca                                                                                                                        >  2:157541/1‑140 (MQ=255)
       ccTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCaacaca                                                                                                                        >  1:68338/1‑140 (MQ=255)
            gCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTgc                                                                                                                   <  2:273077/140‑1 (MQ=255)
                tGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCa                                                                                                               >  2:1078640/1‑140 (MQ=255)
                tGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCa                                                                                                               <  1:531338/140‑1 (MQ=255)
                       cgcgGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCa                                                                                                        >  1:536546/1‑140 (MQ=255)
                           ggTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGcc                                                                                                                <  2:517663/128‑1 (MQ=255)
                           ggTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGcc                                                                                                                >  1:517663/1‑128 (MQ=255)
                                cGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCCACACATCTGCGCCAATATGCAACCGCTgcg                                                                                               <  1:864315/140‑1 (MQ=255)
                                cGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTgcg                                                                                               <  1:696498/140‑1 (MQ=255)
                                 gTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTgcgg                                                                                              >  2:282585/1‑140 (MQ=255)
                                     cTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCa                                                                                          <  2:719876/140‑1 (MQ=255)
                                      tAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCaa                                                                                         >  2:1187663/1‑140 (MQ=255)
                                            aTCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATaga                                                                                   <  2:299980/140‑1 (MQ=255)
                                            aTCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATaga                                                                                   <  1:608726/140‑1 (MQ=255)
                                            aTCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATaga                                                                                   <  1:1187663/140‑1 (MQ=255)
                                            aTCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATaga                                                                                   <  2:536546/140‑1 (MQ=255)
                                             tCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTg                                                                                                 >  1:430283/1‑125 (MQ=255)
                                             tCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTg                                                                                                 <  2:430283/125‑1 (MQ=255)
                                             tCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATagag                                                                                  >  2:328647/1‑140 (MQ=255)
                                                  cagcGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCcac                                                                             <  2:570823/140‑1 (MQ=255)
                                                       aTTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGaa                                                                        >  2:1937/1‑140 (MQ=255)
                                                            aCCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTAc                                                                   <  1:1189853/140‑1 (MQ=255)
                                                                  aCCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAgcgc                                                             >  2:332571/1‑140 (MQ=255)
                                                                      aGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGTTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTa                                                         >  1:645101/1‑140 (MQ=255)
                                                                       ggTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTAt                                                        >  1:1183473/1‑140 (MQ=255)
                                                                            cTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATcgcgc                                                   >  2:314847/1‑140 (MQ=255)
                                                                                  tCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATg                                             >  2:456762/1‑140 (MQ=255)
                                                                                    aGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGct                                           >  2:903613/1‑140 (MQ=255)
                                                                                     gCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGctc                                          >  1:945930/1‑140 (MQ=255)
                                                                                     gCGTAGTCAGAGACGCTAACCACGCAGCTGCCTGACGCTGTTGCCAGCATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGctc                                          >  2:509700/1‑140 (MQ=255)
                                                                                        tAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTcc                                       <  2:1183473/140‑1 (MQ=255)
                                                                                         aGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGct                                           <  2:1161114/135‑1 (MQ=255)
                                                                                         aGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGct                                           >  1:1161114/1‑135 (MQ=255)
                                                                                         aGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTccc                                      <  1:314847/140‑1 (MQ=255)
                                                                                                 aCGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGa                              >  2:133614/1‑140 (MQ=255)
                                                                                                  cGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGac                             >  1:288834/1‑140 (MQ=255)
                                                                                                  cGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGac                             <  1:1110518/140‑1 (MQ=255)
                                                                                                     tAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGacaac                          >  1:114608/1‑140 (MQ=255)
                                                                                                      aaCCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGacaaca                         >  2:20228/1‑140 (MQ=255)
                                                                                                       aCCACGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCAACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGAcaacaa                        <  2:68338/140‑1 (MQ=255)
                                                                                                           cGCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATaaa                    >  2:880088/1‑140 (MQ=255)
                                                                                                            gCAGGTGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAAc                   <  1:407406/140‑1 (MQ=255)
                                                                                                                 tGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGac                             >  1:933920/1‑125 (MQ=255)
                                                                                                                 tGCCTGACGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGac                             <  2:933920/125‑1 (MQ=255)
                                                                                                                        cGCTGTTGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGa       >  2:707829/1‑140 (MQ=255)
                                                                                                                             ttGCCAGTATTACCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTc  <  2:936120/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                
AAGACTTCCTGAGCTTTGTTTTGCGCGGGTAACGTCGCTAATGAATCAAGCAGCGATTTAACCGTTACCGAGGTTTCTGGCTTCAGCGTAGTCAGAGACGCTAACCACGCAGGTGCCTGACGCTGTTGCCAGTATTGCCCCCACACATCTGCGCCAATATGCAACCGCTGCGGCGCAAAATAGAGCCCACAGGAACGTACCAGCGCGGTATCGCGCGTTATGCTCTCCCCAGCCTGACAACAATAAACCTCTTCCGGGGAGTTTC  >  NC_000913/2329283‑2329547

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