Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 361,682 Δ1 bp coding (179/612 nt) lacA ← galactoside O‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,6790T.95.2% 74.1 / ‑1.8 21coding (182/612 nt)lacAgalactoside O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (7/13);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCA  >  NC_000913/361547‑361806
                                                                                                                                    |                                                                                                                               
gCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGt                                                                                                                                 <  2:465307/133‑1 (MQ=255)
       gTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCCAGCGTTTCcccc                                                                                                                          >  1:137424/1‑133 (MQ=255)
       gTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTcccct                                                                                                                          >  2:566776/1‑132 (MQ=255)
            ttATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGt                                                                                                                    >  2:566595/1‑133 (MQ=255)
                cACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCa                                                                                                                >  1:369205/1‑133 (MQ=255)
                cACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTAACAAGGAGTATCCCCTGCCGCGGCa                                                                                                                >  1:301777/1‑133 (MQ=255)
                    gATTGTTACCGTGTACTCTTCGCCAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACa                                                                                                            <  1:260735/133‑1 (MQ=255)
                       tGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTAACCAGGCGGTTTCCCTACCGTGGGAAAAAttt                                                                                                         >  2:120406/1‑133 (MQ=255)
                         ttACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCCGAGGTTTTCCCTACCGGGGCAAACAttttt                                                                                                       >  2:346401/1‑131 (MQ=255)
                                      aTCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTATCCAGGCGGTTCCCCTACCGCGGCAAAAATTTCTTTAATTAGGCttt                                                                                          >  2:535998/1‑133 (MQ=255)
                                       tCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACCTTTCTTTAATCCGGCTTtc                                                                                         <  2:369205/133‑1 (MQ=255)
                                              tGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCttttt                                                                                  >  2:495670/1‑133 (MQ=255)
                                                      ttGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTc                                                                          <  1:132569/133‑1 (MQ=255)
                                                        gtaaTTTGCATAAAACTTGCGGCCTAAAAGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTg                                                                        <  2:301777/131‑1 (MQ=255)
                                                        gAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTg                                                                        <  1:566595/133‑1 (MQ=255)
                                                        gAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTg                                                                        <  1:120406/133‑1 (MQ=255)
                                                              tGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTgatggat                                                                  <  2:5038/133‑1 (MQ=255)
                                                              tGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTgatggat                                                                  <  1:346401/133‑1 (MQ=255)
                                                                 aTAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCg                                                               <  1:566776/133‑1 (MQ=255)
                                                                       aTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGAt                                                         <  1:495670/133‑1 (MQ=255)
                                                                               agatatGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCa                                                 <  2:147091/131‑1 (MQ=255)
                                                                                        tGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAAc                                        >  1:84298/1‑133 (MQ=255)
                                                                                            ggAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGttt                                    >  1:569931/1‑133 (MQ=255)
                                                                                            ggAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGGGGCAAACATTTCTTTTATCAGGCTTTTTCTTTTTTTAACTTCTGATGGATGCGAGGGGTTAAACTCATACATTAACGttt                                    >  1:413100/1‑133 (MQ=255)
                                                                                                   tAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTTAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACg                             >  2:513791/1‑133 (MQ=255)
                                                                                                           aTAGACAGGCGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTTTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATAAACGTTTTTCCACGAAGTCttt                     >  1:244520/1‑133 (MQ=255)
                                                                                                                    cGGTTCTACCCAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTa            <  2:84298/133‑1 (MQ=255)
                                                                                                                              cAGGCGTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGca  <  2:413100/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                                               
GCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCA  >  NC_000913/361547‑361806

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