Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 708,774 G→A A147A (GCG→GCA chiP → chitooligosaccharide outer membrane channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913708,7740GA95.5% 58.0 / ‑4.1 23A147A (GCG→GCAchiPchitooligosaccharide outer membrane channel
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (13/8);  total (14/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 = 6.49e-01

ACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCT  >  NC_000913/708539‑709003
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
aCTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACAtt                                                                                                                                                                                                                        >  1:13241/1‑251 (MQ=255)
            aTGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGt                                                                                                                                                                                                            <  2:137238/251‑1 (MQ=255)
                  aCCTCCATTTTCAGTCCGGTTATGCTGGTGATATTTTCCGGCCTTATATTGCTGCGTTTACGGCGAATTAAATTGCAGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCCTTTTCAAAAAGTAATAAAGACTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACg                                                                                                                                                                                                      <  2:114639/251‑1 (MQ=255)
                                                  atGTTCGGCCTTGATATTGCTGCGTTGACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAgg                                                                                                                                                                      <  1:11551/251‑1 (MQ=255)
                                                          gCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACATCTCCCCCCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACCAAAGCGGTATCAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGTATACATTCAGCCAACTGGTCAAACGCTTTTATCCCCGCACTGGAGCTTTATGCCAGAGACTTAt                                                                                                                                                               >  2:173435/1‑251 (MQ=255)
                                                                              aCGGCGATTGAAAGGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTAGCAGGGGGCGGAAGCCGGGGc                                                                                                                                           <  1:381554/251‑1 (MQ=255)
                                                                               cGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCg                                                                                                                                          <  1:320903/251‑1 (MQ=255)
                                                                                    aTTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGCCTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGAGGCGAAttt                                                                                                                                     >  1:36367/1‑251 (MQ=255)
                                                                                                         gACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGc                                                                                                                >  1:10826/1‑251 (MQ=255)
                                                                                                           cAGCTCCCACCCGAACGAAATCGCGTTTTAAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCAAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCgg                                                                                                                                              >  1:305886/1‑219 (MQ=255)
                                                                                                                                             aGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACa                                                                             >  1:65364/1‑250 (MQ=255)
                                                                                                                                                     atccTAAGACGAAGACTGGTCCGGCGCCAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAACTTTAAATACGGTCCGGTTTGGGCAAGGGCAGGGTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGAAGCTTTATGCCAGGGACTTATCAGGGAGCGGACGCCGGGGCGAATTTTGATTACGGCGATGATGGTGCGTTGCGTTTCTCCGACATGTGGGCCACCGAATACAAAgagacgt                                                                    >  2:276422/3‑245 (MQ=255)
                                                                                                                                                        cTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGc                                                                 >  1:203994/1‑251 (MQ=255)
                                                                                                                                                            gACGAAGACTGGTCCGTGCACAAAAGCGGTATAATTCTGTATAAAGCCGCGGCCAAATTTAAACACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTTCTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATGACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGGGGACCAACGAATACAAAGCGCCGTGGCATCt                                                             <  2:250452/251‑1 (MQ=255)
                                                                                                                                                            gACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCt                                                             >  1:259396/1‑251 (MQ=255)
                                                                                                                                                                          cGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCg                                                                                                                                          >  1:303983/1‑160 (MQ=255)
                                                                                                                                                                          cGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCg                                                                                                                                          <  2:303983/160‑1 (MQ=255)
                                                                                                                                                                                          tataAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGAt                                 >  1:284170/1‑250 (MQ=255)
                                                                                                                                                                                                           cTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCTACTGGTCAAACGCTGTTAGCGCCGCACTCGATCTTTATGCCAGGTACTTACCAGGGGGCGGAAGCCGGGGCGCATTTTGATTAc                                                                                                                              <  2:166715/140‑1 (MQ=255)
                                                                                                                                                                                                           cTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTAc                                                                                                                              >  1:166715/1‑140 (MQ=255)
                                                                                                                                                                                                                    aaTTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGGGGCATCTGGAAATGGATGAGTTTTATCAGAACGATCAAACCACCAAAGTTGATTATCTGCAc      >  1:121372/1‑251 (MQ=255)
                                                                                                                                                                                                                     aTTTAAATTCGGTCCGGTTTGGGCCAGTGCTGGTTTCTTTCATCCAACTTGTCAAACGTGGTTAGCGCCGCACTGGAGCGTTATGCCAGGTACTTATCAGGGTGCGGAAGCCGGGTCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATTTGGTCCAACGAATACAAAGCGCCGTGGCATCTTGAAAGGGATGAGTTTTATCATAACGATAAAGCCACCAAAGTTGATTATTTGCACt     <  2:273064/251‑1 (MQ=255)
                                                                                                                                                                                                                        tAAATACGGTCCGGTTGGGGCAAGGGCAGGATACATTCAGACAACTGGTCAAACGCTGTTAGCGACGCACTGGAGCTTTATGCCTGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGGTGATGCGTTGGGTTTCCCCTACATGTGGTCCAACGTATAAAAAGTGCCGGGGAATCTGGAAAGGGATGCGGTTTATCAGCACGATAAAACCACAAACGTTGATTAGTTGCACTCCt  >  2:342671/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     
ACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCT  >  NC_000913/708539‑709003

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