Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,344,042 T→C G33G (GGT→GGC)  lptB → lipopolysaccharide transport system ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,344,0420TC97.6% 134.5 / NA 41G33G (GGT→GGC) lptBlipopolysaccharide transport system ATP binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (18/22);  total (18/23)
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

CGGCACAGAAGAAGGGTAATTAATTCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGTCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAACATCATTATTGATGATGACGATATCAGTCTGCTGCCTCTG  >  NC_000913/3343915‑3344159
                                                                                                                               |                                                                                                                     
cGGCACAGAAGAAGGGTAATTAATTCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCctgct                                                                                                                  <  2:368588/133‑1 (MQ=255)
 ggCACAGAAGAAGGGTAATTAATTCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTGCGGGGATATTGTCAGCctgctg                                                                                                                 >  1:200717/1‑133 (MQ=255)
 ggCACAGAAGAAGGGTAATTAATTCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCctgctg                                                                                                                 >  1:387394/1‑133 (MQ=255)
        aagaagAGTAATTAATTCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTTTCGGCCTGCTGGGGCCaa                                                                                                          >  1:223808/1‑133 (MQ=255)
         agaagGGTAATTAATTCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCaaa                                                                                                         <  1:23917/133‑1 (MQ=255)
           aagGGTAATTAATTCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACg                                                                                                       >  1:347034/1‑133 (MQ=255)
             gggTAATTAATTCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGt                                                                                                     <  2:101164/133‑1 (MQ=255)
                  attaattCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCgg                                                                                                >  1:57258/1‑133 (MQ=255)
                  attaattCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCgg                                                                                                >  2:388987/1‑133 (MQ=255)
                   ttaattCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGt                                                                                               <  2:341044/133‑1 (MQ=255)
                        tCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGTGGCAAAACGGTGCGGGTAAGac                                                                                          >  1:177601/1‑133 (MQ=255)
                        tCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGac                                                                                          >  1:430064/1‑133 (MQ=255)
                        tCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGCCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGac                                                                                          >  2:127550/1‑133 (MQ=255)
                         cGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGacc                                                                                         >  2:99667/1‑133 (MQ=255)
                         cGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGacc                                                                                         <  1:388987/133‑1 (MQ=255)
                          gTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGacca                                                                                        <  1:78409/133‑1 (MQ=255)
                              tGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACcaccac                                                                                    <  1:312448/133‑1 (MQ=255)
                                      ttAACTGCAATGAACCTTGCAAAACCCTCTAATGGCCGTCGAGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACa                                                                            <  2:135257/133‑1 (MQ=255)
                                      ttAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACa                                                                            <  1:127550/133‑1 (MQ=255)
                                      ttAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACa                                                                            <  1:178445/133‑1 (MQ=255)
                                             cAAAGAACCTTGCAACAGCCTATACAGGCCGTAGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGt                                                                     <  2:177601/133‑1 (MQ=255)
                                               aaGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAg                                                                   >  1:337909/1‑133 (MQ=255)
                                                     cTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGGCCACCACATTCGACTTGGTTGAAGGCATtg                                                             >  1:178139/1‑133 (MQ=255)
                                                        gCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGc                                                          >  1:403321/1‑133 (MQ=255)
                                                            aaGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCcgcg                                                      >  2:91492/1‑133 (MQ=255)
                                                                 tataAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATg                                                 >  2:438955/1‑133 (MQ=255)
                                                                      aGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGc                                            <  1:99667/133‑1 (MQ=255)
                                                                      aGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGc                                            <  1:315610/133‑1 (MQ=255)
                                                                       ggCCGTCGCGTGGTATAAGACGTCTGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCa                                           <  2:200717/133‑1 (MQ=255)
                                                                         ccGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAAc                                         <  2:420173/133‑1 (MQ=255)
                                                                          cGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAAca                                        >  2:358060/1‑133 (MQ=255)
                                                                          cGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTCTAGGCATCGTGACGAGCGATGGGGGCAAca                                        >  1:165678/1‑133 (MQ=255)
                                                                           gTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAAcat                                       <  2:403321/133‑1 (MQ=255)
                                                                                      gaagaCGTCATCCTGACCGTCAACTCCGAGGAAATTGTCGGTCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAACATCATTATtgatg                            <  2:178139/133‑1 (MQ=255)
                                                                                      gaagaCGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAACATCATTATtgatg                            <  2:427999/133‑1 (MQ=255)
                                                                                          aCGTCAGCCTGACCGTCATCTCGGGGGAAAATGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAACATCATTATtgatgatga                        <  1:180177/133‑1 (MQ=255)
                                                                                          aCGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCTCGCGATGCGGGCAACATCATTATtgatgatga                        <  2:430064/133‑1 (MQ=255)
                                                                                                 ccTGACCGTCAACTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAACATCATTATTGATGATGACGATATc                 <  2:347034/133‑1 (MQ=255)
                                                                                                           aaCTCCGGGGCAATTATCGGCCTGATGGCGCCCAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAACATCATTATTGATGATGACGATATCAGTCtgctgc       <  2:154325/133‑1 (MQ=255)
                                                                                                           aaCTCCGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAACATCATTATTGATGATGACGATATCAGTCtgctgc       <  1:480615/133‑1 (MQ=255)
                                                                                                                cGGGGAAATTGTCGGCCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAACATCATTGTTGATGATGACGATATGAGTCGGCTGCCTCtg  >  2:144026/1‑133 (MQ=255)
                                                                                                                               |                                                                                                                     
CGGCACAGAAGAAGGGTAATTAATTCGTTATGGCAACATTAACTGCAAAGAACCTTGCAAAAGCCTATAAAGGCCGTCGCGTGGTAGAAGACGTCAGCCTGACCGTCAACTCCGGGGAAATTGTCGGTCTGCTGGGGCCAAACGGTGCCGGTAAGACCACCACTTTCTACATGGTTGTAGGCATTGTGCCGCGCGATGCGGGCAACATCATTATTGATGATGACGATATCAGTCTGCTGCCTCTG  >  NC_000913/3343915‑3344159

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