Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,796,375 T→A 9.8% 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 9.8% 135.6 / 6.2 41I138F (ATT→TTT) ygaUuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base A (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

ATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGT  >  NC_000913/2796254‑2796522
                                                                                                                         |                                                                                                                                                   
aTACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAt                                                                                                                                                     >  2:815988/1‑122 (MQ=255)
aTACCGCGCGCGCCAGTCAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACCTGCCCCGGATAAAt                                                                                                                                                     <  1:815988/122‑1 (MQ=255)
                         gcgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAtttt                                                                                                                                                  <  2:980361/100‑1 (MQ=255)
                         gcgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAtttt                                                                                                                                                  >  1:980361/1‑100 (MQ=255)
                                  tttACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGc                                                                                         >  2:407561/1‑148 (MQ=255)
                                          cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCg                                                                                                                                             >  3:161870/1‑88 (MQ=255)
                                          cGCCAATGCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCg                                                                                                                                             <  4:161870/88‑1 (MQ=255)
                                            ccAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCgg                                                                                                                             <  2:623538/102‑1 (MQ=255)
                                            ccAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCgg                                                                                                                             >  1:623538/1‑102 (MQ=255)
                                            ccAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGa                                                                                                            >  1:297141/1‑119 (MQ=255)
                                            ccAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGa                                                                                                            <  2:297141/119‑1 (MQ=255)
                                            ccAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCg                                                                                                                                             <  1:1042174/86‑1 (MQ=255)
                                            ccAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCg                                                                                                                                             >  2:1042174/1‑86 (MQ=255)
                                              aaTTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                                                                             >  3:440756/1‑148 (MQ=255)
                                                   ccATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGtt                                                                        <  1:407561/148‑1 (MQ=255)
                                                    cATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCg                                                                                                                <  1:2621/107‑1 (MQ=255)
                                                    cATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCg                                                                                                                >  2:2621/1‑107 (MQ=255)
                                                     aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTAt                                                                                                                        >  2:576426/1‑98 (MQ=255)
                                                     aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTAt                                                                                                                        <  1:576426/98‑1 (MQ=255)
                                                     aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAg                                                                               <  4:36268/139‑1 (MQ=255)
                                                     aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAg                                                                               >  3:36268/1‑139 (MQ=255)
                                                     aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTg                                                                      >  3:769991/1‑148 (MQ=255)
                                                          ccAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTAc                                                                                                 >  2:168032/1‑116 (MQ=255)
                                                          ccAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTAc                                                                                                 <  1:168032/116‑1 (MQ=255)
                                                                        ctggcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAgtgt                                                   <  4:440756/148‑1 (MQ=255)
                                                                          ggcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGt                                                                         <  1:207527/124‑1 (MQ=255)
                                                                          ggcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGt                                                                         >  2:207527/1‑124 (MQ=255)
                                                                              aTTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAtttt                                                                                                                                                  <  3:176530/47‑1 (MQ=255)
                                                                              aTTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAtttt                                                                                                                                                  >  4:176530/1‑47 (MQ=255)
                                                                                tAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAg                                                                                                                                   >  1:101509/1‑60 (MQ=255)
                                                                                tAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAg                                                                                                                                   <  2:101509/60‑1 (MQ=255)
                                                                                    tactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAgtgt                                                   <  3:1090436/136‑1 (MQ=255)
                                                                                    tactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAgtgt                                                   >  4:1090436/1‑136 (MQ=255)
                                                                                     actCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTa                                        >  2:743554/1‑146 (MQ=255)
                                                                                     actCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTa                                        <  1:743554/146‑1 (MQ=255)
                                                                                          tCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTAc                                                                                                 <  4:844466/84‑1 (MQ=255)
                                                                                          tCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTAc                                                                                                 >  3:844466/1‑84 (MQ=255)
                                                                                             ggAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATaaa                              <  4:769991/148‑1 (MQ=255)
                                                                                              gAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATaaaa                             >  2:158320/1‑148 (MQ=255)
                                                                                                   cGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActggc                        <  1:158320/148‑1 (MQ=255)
                                                                                                       acacTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAg                                                                                               >  1:1093327/1‑73 (MQ=255)
                                                                                                       acacTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAg                                                                                               <  2:1093327/73‑1 (MQ=255)
                                                                                                       acacTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGa                                                                              <  1:314571/90‑1 (MQ=255)
                                                                                                       acacTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGa                                                                              >  2:314571/1‑90 (MQ=255)
                                                                                                                 cGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGgt          >  3:965824/1‑148 (MQ=255)
                                                                                                                 cGGATAAAATTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActg                          >  3:66771/1‑132 (MQ=255)
                                                                                                                 cGGATAAAATTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActg                          <  4:66771/132‑1 (MQ=255)
                                                                                                                         ttttATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCt             <  3:552671/137‑1 (MQ=255)
                                                                                                                         ttttATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCt             >  4:552671/1‑137 (MQ=255)
                                                                                                                         ttttATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGt  <  4:965824/148‑1 (MQ=255)
                                                                                                                         |                                                                                                                                                   
ATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGT  >  NC_000913/2796254‑2796522

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: