Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,796,375 T→A 8.2% I138F (ATT→TTT)  ygaU ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,796,3750TA 8.2% 161.6 / 5.6 49I138F (ATT→TTT) ygaUuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (21/24);  new base A (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

TGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTT  >  NC_000913/2796230‑2796523
                                                                                                                                                 |                                                                                                                                                    
tGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAtttt                                                                                                                                                   >  1:393615/1‑149 (MQ=255)
  aCTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTAt                                                                                                                                                 <  1:452615/149‑1 (MQ=255)
    tGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATcc                                                                                                                                               <  4:442460/149‑1 (MQ=255)
      cggcggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCgg                                                                                                                                             <  4:231585/149‑1 (MQ=255)
      cggcggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCgg                                                                                                                                             <  4:403671/149‑1 (MQ=255)
             aCTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCttttt                                                                                                                                      <  2:393615/149‑1 (MQ=255)
                ggCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGc                                                                                                                                   >  2:94669/1‑149 (MQ=255)
                         tACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTa                                                                                                                          >  3:52783/1‑149 (MQ=255)
                         tACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTa                                                                                                                          <  1:472306/149‑1 (MQ=255)
                              cgcgcgcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGc                                                                                                                     <  2:282272/149‑1 (MQ=255)
                                 gcgcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTc                                                                                                                  <  4:184963/149‑1 (MQ=255)
                                           gTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCTACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAttttt                                                                                                        <  1:424657/149‑2 (MQ=255)
                                            ttAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTAt                                                                                                       <  3:284226/149‑1 (MQ=255)
                                            ttAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTAt                                                                                                       <  4:327925/149‑1 (MQ=255)
                                              aaTGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGc                                                                                                                     <  1:455498/133‑1 (MQ=255)
                                              aaTGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGc                                                                                                                     >  2:455498/1‑133 (MQ=255)
                                                tGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTAt                                                                                                                         <  4:141726/127‑1 (MQ=255)
                                                tGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTAt                                                                                                                         >  3:141726/1‑127 (MQ=255)
                                                 gcgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTAtt                                                                                                                        <  3:458532/127‑1 (MQ=255)
                                                 gcgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTAtt                                                                                                                        >  4:458532/1‑127 (MQ=255)
                                                  cgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACa                                                                                                 <  3:380556/149‑1 (MQ=255)
                                                  cgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACa                                                                                                 <  3:443508/149‑1 (MQ=255)
                                                     cAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGAt                                                                                              >  1:423898/1‑149 (MQ=255)
                                                       gCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTa                                                                                            <  2:25182/149‑1 (MQ=255)
                                                         aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGc                                                                                          >  2:44736/1‑149 (MQ=255)
                                                         aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGc                                                                                          <  4:456603/149‑1 (MQ=255)
                                                                     cAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                                                                              <  1:94669/149‑1 (MQ=255)
                                                                           ccATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGttt                                                                        <  1:44736/149‑1 (MQ=255)
                                                                                 tCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAAt                                                                  >  4:25643/1‑149 (MQ=255)
                                                                                   cagcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATgg                                                                >  3:245143/1‑149 (MQ=255)
                                                                                    agcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATAGGTTTATTCGCTTCGAAGATTTTATTGTATAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGc                                                               >  4:423657/1‑149 (MQ=255)
                                                                                       agGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACt                                                            >  4:283290/1‑149 (MQ=255)
                                                                                           gATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCaga                                                        <  4:52783/149‑1 (MQ=255)
                                                                                                ctggcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTc                                                   <  2:423898/149‑1 (MQ=255)
                                                                                                     aaTTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAg                                              >  3:66832/1‑149 (MQ=255)
                                                                                                           ctactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGa                                                                               >  4:234575/1‑110 (MQ=255)
                                                                                                           ctactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGa                                                                               <  3:234575/110‑1 (MQ=255)
                                                                                                           ctactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTgg                                                                      >  1:298496/1‑119 (MQ=255)
                                                                                                           ctactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTgg                                                                      <  2:298496/119‑1 (MQ=255)
                                                                                                               tctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGt                                    >  2:384932/1‑149 (MQ=255)
                                                                                                                 ttCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGtat                                  >  3:429782/1‑149 (MQ=255)
                                                                                                                            gCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActggctg                       >  3:223777/1‑149 (MQ=255)
                                                                                                                              aaCACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACtggctggc                     >  4:172318/1‑149 (MQ=255)
                                                                                                                              aaCACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACtggctggc                     >  3:225219/1‑149 (MQ=255)
                                                                                                                              aaCACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACtggctggc                     >  1:473750/1‑149 (MQ=255)
                                                                                                                                cacTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGaaa                                                                   <  3:68758/101‑1 (MQ=255)
                                                                                                                                cacTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGaaa                                                                   >  4:68758/1‑101 (MQ=255)
                                                                                                                                        ccGGATAAAATTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTa                                                                                                   <  4:203041/61‑1 (MQ=255)
                                                                                                                                        ccGGATAAAATTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTa                                                                                                   >  3:203041/1‑61 (MQ=255)
                                                                                                                                        ccGGATAAAATTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACCTAACGGTATAAAACTGGCTGGCAGTGGCTGgt           >  2:436798/1‑149 (MQ=255)
                                                                                                                                                 ttttATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGtt  >  1:452862/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTT  >  NC_000913/2796230‑2796523

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