Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,601,771 A→G E94E (GAA→GAG)  lsrA → Autoinducer‑2 ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,601,7710AG96.8% 100.0 / NA 31E94E (GAA→GAG) lsrAAutoinducer‑2 ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (19/11);  total (20/11)
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

GTAAATCGACGTTAATGAAGATTATTGCCGGTATTACCCCTGCTGATAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCT  >  NC_000913/1601635‑1601899
                                                                                                                                        |                                                                                                                                
gTAAATCGACGTTAATGAAGATCATTGCCGGGATTCCCCTTGCTGATATCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCg                                                                                                                               <  2:398101/140‑1 (MQ=255)
          gTTAATGAAGATTATTGCCGGTATTACCCCTGCTGATAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTcc                                                                                                                     >  2:674065/1‑140 (MQ=255)
           ttAATGAAGATTATTGCCGGTATTACCCCTGCTGATAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGAGTATTTATCTCTTTACCCAGGAGCCGCTGCTTTTccc                                                                                                                    >  1:247013/1‑140 (MQ=255)
               tGAAGATTATTGCCGGTATTACCCCTGCTGATAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGc                                                                                                                >  2:175087/1‑140 (MQ=255)
                   gATTATTGCCGGTATTACCCCTGCTGATAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGt                                                                                                            <  1:427985/140‑1 (MQ=255)
                              acagtACCCCTGCTGGTAGCGGTACGCTGGAGATTGATGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGaa                                                                                                 <  1:326136/136‑1 (MQ=255)
                                         gctgATAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATc                                                                                           >  1:210798/1‑135 (MQ=255)
                                         gctgATAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATc                                                                                           <  2:210798/135‑1 (MQ=255)
                                         gctgATAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGtt                                                                                      >  2:832358/1‑140 (MQ=255)
                                             aTAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAACCGCTGCTCTTCTCAAGCCTGTCGACAAAATAAAACATCCTGttctgt                                                                                  >  2:326136/1‑136 (MQ=255)
                                                         tGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCaaaaaaa                                                                      <  1:832358/140‑1 (MQ=255)
                                                         tGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCaaaaaaa                                                                      >  2:211833/1‑140 (MQ=255)
                                                         tGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTATCGGTAAAAGAAAACACCCTGCTTGGGCTGGCaaaaaaa                                                                      >  2:314610/1‑140 (MQ=255)
                                                                tGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGc                                                                  >  1:202391/1‑137 (MQ=255)
                                                                tGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGc                                                                  <  2:202391/137‑1 (MQ=255)
                                                                  aGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTATCGATAAAGGAACACAACCTGTTTGAGCTGGCAAAAAAACAGCTCTcc                                                             >  2:251187/1‑140 (MQ=255)
                                                                     gCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTATTTGGGCTGGC‑AAAAAACAGCTCTCCATGc                                                         >  2:647191/1‑140 (MQ=255)
                                                                               cGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATga                                                >  1:22464/1‑140 (MQ=255)
                                                                                 tCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATgaag                                              >  1:702411/1‑140 (MQ=255)
                                                                                  cAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGC‑AAAAAACAGCTCTCCATGCAGAAAATgaagaa                                            <  1:647191/140‑1 (MQ=255)
                                                                                              aGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTggcgg                                 >  1:29207/1‑140 (MQ=255)
                                                                                                ttCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGgcggcg                               >  1:126383/1‑140 (MQ=255)
                                                                                                         atcaGCTGGGTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGc                      >  2:674150/1‑140 (MQ=255)
                                                                                                                  gTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGAt             <  1:674150/140‑1 (MQ=255)
                                                                                                                  gTATTTATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGAt             <  1:175087/140‑1 (MQ=255)
                                                                                                                       tATCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCa        <  2:702411/140‑1 (MQ=255)
                                                                                                                        aTCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGttt                >  1:424420/1‑131 (MQ=255)
                                                                                                                        aTCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGttt                <  2:424420/131‑1 (MQ=255)
                                                                                                                        aTCTCGTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCAt       >  2:821105/1‑140 (MQ=255)
                                                                                                                             gTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCt  >  1:522316/1‑140 (MQ=255)
                                                                                                                             gTTCCCCAGGAGCCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCt  >  1:719200/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                
GTAAATCGACGTTAATGAAGATTATTGCCGGTATTACCCCTGCTGATAGCGGTACGCTGGAGATTGAGGGCAACAACTACGTCAGATTAACGCCAGTTCATGCTCATCAGCTGGGTATTTATCTCGTTCCCCAGGAACCGCTGCTTTTCCCAAGCCTGTCGATAAAAGAAAACATCCTGTTTGGGCTGGCAAAAAAACAGCTCTCCATGCAGAAAATGAAGAACTTGCTGGCGGCGCTGGGCTGCCAGTTTGATCTGCATAGTCT  >  NC_000913/1601635‑1601899

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