Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 986,740 A→G G81G (GGT→GGC)  ompF ← outer membrane porin F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,7400AG97.0% 106.6 / NA 33G81G (GGT→GGC) ompFouter membrane porin F
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/15);  minor base T (0/1);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCT  >  NC_000913/986612‑986869
                                                                                                                                |                                                                                                                                 
tAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCAt                                                                                                                               <  1:409976/133‑1 (MQ=255)
 aGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATa                                                                                                                              >  2:112546/1‑133 (MQ=255)
          cGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAg                                                                                                                     >  2:361234/1‑133 (MQ=255)
            tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGAt                                                                                                                   >  1:301381/1‑133 (MQ=255)
            tAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGGTATATTCCCACCGGACATAACAGGTCAGAt                                                                                                                   >  1:173433/1‑133 (MQ=255)
                cGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGa                                                                                                               <  2:105604/133‑1 (MQ=255)
                        ccAACGTCAGCGTATTTAAGGCCCGCGAATACCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGAttt                                                                                                       >  1:425118/1‑133 (MQ=255)
                                     aTTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGCAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCttt                                                                                          >  1:408693/1‑133 (MQ=255)
                                               cgcgAATGTCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                                                                                <  1:189402/133‑1 (MQ=255)
                                               cgcgAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                                                                                >  2:427360/1‑133 (MQ=255)
                                               cgcgAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGa                                                                                <  1:427360/133‑1 (MQ=255)
                                                 cgAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTTAGTTTCCCGTTTAAGACCAAGACg                                                                              >  2:223247/1‑133 (MQ=255)
                                                        cAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAg                                                                       <  1:361234/133‑1 (MQ=255)
                                                              cgTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg                                                                 <  2:375889/133‑1 (MQ=255)
                                                                ttttGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTc                                                               >  2:452871/1‑133 (MQ=255)
                                                                ttttGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTc                                                               >  2:330159/1‑133 (MQ=255)
                                                                     ttACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGTCGGACAGAGGTCAAGTCGCCAt                                                          >  1:251717/1‑133 (MQ=255)
                                                                        ccAGTTTGAAAGTCAGCGCCTTCAGACTCGTTTCCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGc                                                       <  2:246248/133‑1 (MQ=255)
                                                                               gAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTa                                                >  1:481591/1‑133 (MQ=255)
                                                                                      aGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGttt                                         >  2:328524/1‑133 (MQ=255)
                                                                                       gcgcCTTCAGTGTTGCTACCCTGGCAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                        <  2:173433/133‑1 (MQ=255)
                                                                                       gcgcCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                        <  1:328524/133‑1 (MQ=255)
                                                                                       gcgcCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                        <  2:481591/133‑1 (MQ=255)
                                                                                              cAGAGTTGTTACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGtt                                 >  1:57035/1‑133 (MQ=255)
                                                                                                      ttACCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTgg                         <  2:61166/133‑1 (MQ=255)
                                                                                                      ttACCCTGGAAGTTAAATTCCCACGGGCTATATCCGGTCAGTTCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTgg                         <  2:188503/133‑1 (MQ=255)
                                                                                                        aCCCTGGAAGTTATATTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCCCCGTAACTGTTGTCCCCGTTAGCCTTGGaa                       >  1:216672/1‑133 (MQ=255)
                                                                                                                   tataTTCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAg            <  1:112546/133‑1 (MQ=255)
                                                                                                                      aatCCCACTCTCCATGACCGGTCAAATCGGAATTGATTTGATTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGAcc         <  2:251717/131‑1 (MQ=255)
                                                                                                                        tCCCACTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCaa       <  1:290255/133‑1 (MQ=255)
                                                                                                                            aCTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGc   >  2:374054/1‑133 (MQ=255)
                                                                                                                            aCTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGc   >  2:476863/1‑133 (MQ=255)
                                                                                                                             cTGGCCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAATACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCt  <  2:301381/133‑1 (MQ=255)
                                                                                                                                |                                                                                                                                 
TAGTTACGGCCGTAATCGAAAGAACCAACGTCAGCGTATTTAAGACCCGCGAATGCCAGACGCGTTTTGTTACCAGTTTGAGCGTCAGCGCCTTCAGAGTTGTTACCCTGGAAGTTATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCT  >  NC_000913/986612‑986869

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