Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,454,626 A→T L48L (CTA→CTT USA300HOU_2315 → hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,454,6260AT86.4% 54.5 / ‑0.4 22L48L (CTA→CTTUSA300HOU_2315hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base T (15/4);  total (18/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAACAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTC  >  CP000730/2454487‑2454762
                                                                                                                                           |                                                                                                                                        
agaaCAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCTAATACTTTCTTTTTGTAGTTACCGATCTAg                                                                                                                                         >  1:3998/1‑141 (MQ=255)
        aaaaTTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTTAAGTTACAGATCTTGCACCATTa                                                                                                                                 >  2:472222/1‑141 (MQ=255)
         aaaTTTCTGGTTTTCCATGGGCAATGACGAGTATCGTCTTCTTTGTCATTACAATGGCGGAATCTAGTATGCTCAGAGATTTTCAGGCAGTAGTTGGTGTCCAACACTTTATATTGGAAGTTACAGGTCTTGccacatgaa                                                                                                                                >  2:113498/1‑133 (MQ=255)
           aTTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAAtt                                                                                                                              <  2:375330/141‑1 (MQ=255)
              tCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATtgct                                                                                                                           >  2:333290/1‑141 (MQ=255)
                 ggTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCa                                                                                                                        >  2:41720/1‑141 (MQ=255)
                     ttCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCa                                                                                                                    >  1:8093/1‑141 (MQ=255)
                         aTGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCAtttt                                                                                                                >  2:465763/1‑140 (MQ=255)
                           gggCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGta                                                                                                              >  2:304969/1‑141 (MQ=255)
                                aTGTCGATTGTCTTCTTCTTTGTGATTACAAGTTCTTTATCCATTATGCTCAGAGACTTCCAGTAGATAATTGTTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTc                                                                                                         <  2:487057/141‑1 (MQ=255)
                                                 ttttGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCGTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATaaaaa                                                                                        <  2:152843/140‑1 (MQ=255)
                                                              aTGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCa                                                                           >  1:9690/1‑141 (MQ=255)
                                                                        ccATTATGCTCAGAGATTTCCAGTCCATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGGCCAACTTGCAGGTTTACAAt                                                                 >  2:114409/1‑141 (MQ=255)
                                                                                          tCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAg                                               >  2:364991/1‑141 (MQ=255)
                                                                                                             cAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCt                            >  2:279441/1‑141 (MQ=255)
                                                                                                             cAAACACTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCt                            >  2:461392/1‑141 (MQ=255)
                                                                                                                   cTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCACTCATTTGTATACTCGTTTTCAAATATAAAAAAGTCCAACTTGCAGGTTTAAAAGTCTCAATCAGCCCGAAAGTAATTGAACGTCTATTGCTAGCttt                      >  2:364192/1‑141 (MQ=255)
                                                                                                                   cTTTATATTTGAAGTTACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCttt                      >  2:190916/1‑141 (MQ=255)
                                                                                                                              aaGTTACAGATCTTGCACCATTACTTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCtt           >  2:447734/1‑141 (MQ=255)
                                                                                                                                  tACAGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTCTTGCTAGCTTTAATTTTACCTTTAat       >  2:70085/1‑141 (MQ=255)
                                                                                                                                     aGATCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAattat    >  2:74899/1‑141 (MQ=255)
                                                                                                                                       aTCTTGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTc  <  2:386115/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
AGAACAATAAAATTTCTGGTTTTCAATGGGCAATGACGATTTTCGTCTTCTTTGTCATTACAATGGCGTTATCCATTATGCTCAGAGATTTCCAGTCTATAATTGGTGTCAAACACTTTATATTTGAAGTTACAGATCTAGCACCATTAATTGCTGCAATCATTTGTATACTCGTTTTCAAATATAAAAAGGTCCAACTTGCAGGTTTAAAATTCTCAATCAGCCTGAAAGTAATTGAACGTCTATTGCTAGCTTTAATTTTACCTTTAATTATTC  >  CP000730/2454487‑2454762

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: