Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,371,433 (G)5→4 coding (526/1653 nt) arnT → lipid IVA 4‑amino‑4‑deoxy‑L‑arabinosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,371,4290G.88.0% 76.3 / 7.4 25coding (522/1653 nt)arnTlipid IVA 4‑amino‑4‑deoxy‑L‑arabinosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base . (14/8);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CCGTGCTCGATCCGTTTATTGCATTCTGGCTGGTGGCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCGGTTATCAGTTGCGTACT  >  NC_000913/2371302‑2371557
                                                                                                                                |                                                                                                                                
ccGTGCTCGATCCGTTTATTGCATTCTGGCTGGTGGCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATggggtg                                                                                                                             <  1:9367/133‑3 (MQ=255)
ccGTGCTCGATCCGTTTATTGCATTCTGGCTGGTGGCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATggggtg                                                                                                                             <  2:553333/133‑3 (MQ=255)
ccGTGCTCGATCCGTTTATTGCATTCTGGCTGGTGGCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATggggtg                                                                                                                             >  1:48897/1‑131 (MQ=255)
       cGATCCGTTTATTGCATTCTGGCTGGTGGCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCa                                                                                                                     <  2:526050/133‑1 (MQ=255)
              tttATTGCATTCTGGCTGGTGGCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGttt                                                                                                              >  1:393090/1‑133 (MQ=255)
                   tGCATTCTGGCTGGTGGCGGGAATGTGCATGTTCTGGGTGGCAATGCAGGGACAGACGAGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCCTGGGGGGATGACCAAAGGctttattg                                                                                                         >  2:390561/1‑129 (MQ=255)
                   tGCATTCTGGCTGGAGGCGGGAATGTGGAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAAGCGAAAACGCAAGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGATTTCTCg                                                                                                         >  2:622630/1‑133 (MQ=255)
                    gCATTCTGGCTGGTGGCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGc                                                                                                        >  1:586611/1‑133 (MQ=255)
                              tggtggCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTtgccgtg                                                                                              <  2:586611/133‑1 (MQ=255)
                               ggtggCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGGGATAACCCAAGGTTTTCTCGCACTtgccgtgc                                                                                             >  2:241029/1‑133 (MQ=255)
                                gtggCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTtgccgtgcc                                                                                            <  2:414054/133‑1 (MQ=255)
                                                    ctggctggCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGcc                                                                          >  1:507118/1‑133 (MQ=255)
                                                    ctggctggCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCAt                                                                        >  2:391096/1‑133 (MQ=255)
                                                                       caGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACg                                                     <  2:42774/133‑1 (MQ=255)
                                                                        aGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAAc                                                      <  2:507118/133‑1 (MQ=255)
                                                                         gACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCt                                                   >  1:508519/1‑133 (MQ=255)
                                                                              ggAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGaaa                                              >  1:83175/1‑133 (MQ=255)
                                                                              ggAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGaaa                                              >  2:50250/1‑133 (MQ=255)
                                                                                                    tttACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACggctggc                        <  2:451383/133‑1 (MQ=255)
                                                                                                      tACTGCTGGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGgctggctg                      <  1:390561/133‑1 (MQ=255)
                                                                                                          gctgGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCGAAAACGCTGGAAGGATGTCTTTATTGACGGGTGGCTGGCgg                  >  2:167454/1‑133 (MQ=255)
                                                                                                          gctgGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCgg                  <  1:391096/133‑1 (MQ=255)
                                                                                                          gctgGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCgg                  >  1:510637/1‑133 (MQ=255)
                                                                                                          gctgGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCgg                  >  1:108160/1‑133 (MQ=255)
                                                                                                          gctgGGAATCACCTGCGGCATGGGGTGATGACCAAAGGTTTCCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCgg                  >  1:587812/1‑133 (MQ=255)
                                                                                                                     ccTGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCGGTTATCAGTTGc       >  2:68811/1‑133 (MQ=255)
                                                                                                                       tGCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCGGTTATCAGTTGCGt     <  2:557000/133‑1 (MQ=255)
                                                                                                                        gCGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGTCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCGGTTATCAGTTGCGTa    >  1:452198/1‑133 (MQ=255)
                                                                                                                         cGGCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCGGTTATCAGTTGCGTAc   >  2:482346/1‑133 (MQ=255)
                                                                                                                          ggCATGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGTGGGTATCAGTTGCGTACt  >  1:209832/1‑133 (MQ=255)
                                                                                                                                |                                                                                                                                
CCGTGCTCGATCCGTTTATTGCATTCTGGCTGGTGGCGGGAATGTGCAGCTTCTGGCTGGCAATGCAGGCACAGACGTGGAAAGGCAAAAGCGCAGGATTTTTACTGCTGGGAATCACCTGCGGCATGGGGGTGATGACCAAAGGTTTTCTCGCCCTTGCCGTGCCGGTATTAAGCGTGCTGCCATGGGTAGCGACGCAAAAACGCTGGAAAGATCTCTTTATTTACGGCTGGCTGGCGGTTATCAGTTGCGTACT  >  NC_000913/2371302‑2371557

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