Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 553,473 C→A 7.4% M156I (ATG→ATT lpxH ← UDP‑2,3‑diacylglucosamine pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913553,4730CA 7.4% 185.9 / 3.0 53M156I (ATG→ATTlpxHUDP‑2,3‑diacylglucosamine pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base A (2/2);  total (26/27)
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

AAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAT  >  NC_000913/553329‑553600
                                                                                                                                                |                                                                                                                               
aaGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCAtt                                                                                                                              <  4:411960/148‑1 (MQ=255)
                       gggTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTgc                                                                                                           <  4:106070/144‑1 (MQ=255)
                       gggTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTgc                                                                                                           >  3:106070/1‑144 (MQ=255)
                          tATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTg                                                                                                            <  4:790627/140‑1 (MQ=255)
                          tATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCGGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTg                                                                                                            >  3:790627/1‑140 (MQ=255)
                                cGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTc                                                                                                                             <  1:1063808/117‑1 (MQ=255)
                                cGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTc                                                                                                                             >  2:1063808/1‑117 (MQ=255)
                                      tCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGtt                                                                                                              >  1:789717/1‑126 (MQ=255)
                                      tCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGtt                                                                                                              <  2:789717/126‑1 (MQ=255)
                                         gCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgcg                                                                                                                          <  1:696946/111‑1 (MQ=255)
                                         gCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgcg                                                                                                                          >  2:696946/1‑111 (MQ=255)
                                             ttGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCAt                                                                                                                               <  2:934082/102‑1 (MQ=255)
                                             ttGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCAt                                                                                                                               >  1:934082/1‑102 (MQ=255)
                                              tGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGt                                                                                                               <  2:516660/117‑1 (MQ=255)
                                              tGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGt                                                                                                               >  1:516660/1‑117 (MQ=255)
                                              tGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATa                                                                                >  3:974630/1‑148 (MQ=255)
                                               tcACCTGACGTTTTACCATCGCACTGACCACCGCGTTTTGGTTAACGTCCCTGATCGCCAGTGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCTCCCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAg                                                                               <  4:702384/147‑1 (MQ=255)
                                                  ccTGATGTTATTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGGTTTACTGCTGTTGGCTTCTTTGCAGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCga                                                                            >  3:702384/1‑147 (MQ=255)
                                                    tGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGAATgcgcgc                                                                                                                         >  3:737430/1‑101 (MQ=255)
                                                    tGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGAATgcgcgc                                                                                                                         <  4:737430/101‑1 (MQ=255)
                                                            ttCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAg                                                                  >  3:469128/1‑148 (MQ=255)
                                                              ccATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAggg                                                                <  4:469128/148‑1 (MQ=255)
                                                                aTCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGcc                                                                    >  1:1086089/1‑142 (MQ=255)
                                                                aTCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGcc                                                                    <  2:1086089/142‑1 (MQ=255)
                                                                aTCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGtt                                                              >  4:346671/1‑148 (MQ=255)
                                                                  cGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtg                                                            <  4:974630/148‑1 (MQ=255)
                                                                  cGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGAATGcgcgcga                                                                                                                       >  4:239246/1‑88 (MQ=255)
                                                                  cGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGAATGcgcgcga                                                                                                                       <  3:239246/89‑2 (MQ=255)
                                                                               gcgTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGgcgc                                               >  2:1036672/1‑148 (MQ=255)
                                                                                   tttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGcc                                                                    >  1:255948/1‑123 (MQ=255)
                                                                                   tttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGcc                                                                    <  2:255948/123‑1 (MQ=255)
                                                                                    ttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGaaaa                                          <  3:346671/148‑1 (MQ=255)
                                                                                        ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCt                                      >  4:264894/1‑148 (MQ=255)
                                                                                                cATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGtt                                                                                                              >  1:492300/1‑68 (MQ=255)
                                                                                                cATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGtt                                                                                                              <  2:492300/68‑1 (MQ=255)
                                                                                                      cGCCAGCGATTTACTTCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCa                                                                   >  2:853084/1‑105 (MQ=255)
                                                                                                      cGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCa                                                                   <  1:853084/105‑1 (MQ=255)
                                                                                                         cAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGt                                                                            >  3:654021/1‑93 (MQ=255)
                                                                                                         cAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGt                                                                            <  4:654021/93‑1 (MQ=255)
                                                                                                          aGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCgg                    <  1:1036672/148‑1 (MQ=255)
                                                                                                                    tgctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtg                                                            >  1:630855/1‑98 (MQ=255)
                                                                                                                    tgctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtg                                                            <  2:630855/98‑1 (MQ=255)
                                                                                                                      ctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCgtgt        <  3:264894/148‑1 (MQ=255)
                                                                                                                          tGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAg                                                                                           >  1:827994/1‑61 (MQ=255)
                                                                                                                          tGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAg                                                                                           <  2:827994/61‑1 (MQ=255)
                                                                                                                               tcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAAt                                                                                 >  1:832333/1‑66 (MQ=255)
                                                                                                                               tcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAAt                                                                                 <  2:832333/66‑1 (MQ=255)
                                                                                                                               tcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAg                                         >  1:936492/1‑106 (MQ=255)
                                                                                                                               tcttTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAg                                         <  2:936492/106‑1 (MQ=255)
                                                                                                                                      tGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAt  >  4:1122735/1‑138 (MQ=255)
                                                                                                                                      tGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAt  <  3:1122735/138‑1 (MQ=255)
                                                                                                                                       gTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGt                                                                            >  3:673842/1‑63 (MQ=255)
                                                                                                                                       gTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGt                                                                            <  4:673842/63‑1 (MQ=255)
                                                                                                                                                cATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCcgc                            >  4:27950/1‑102 (MQ=255)
                                                                                                                                                cATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCcgc                            <  3:27950/102‑1 (MQ=255)
                                                                                                                                                cATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGt                   <  3:624161/111‑1 (MQ=255)
                                                                                                                                                cATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGt                   >  4:624161/1‑111 (MQ=255)
                                                                                                                                                |                                                                                                                               
AAGTTCATGCACCGCCGGGCGATGGGTATGCCCGTGGATCAGCCATTGCACCTGATGTTTTTCCATCGCACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGGACCTTGGCGCGAAAAGCCTGATAACCCGCGTCATCGGTGCACAGCGTGTCGCCAT  >  NC_000913/553329‑553600

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