Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,649,718 A→G E37G (GAA→GGA)  dicB → Qin prophage; cell division inhibition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,649,7180AG95.2% 135.5 / ‑1.8 42E37G (GAA→GGA) dicBQin prophage; cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (19/21);  total (20/22)
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.79e-01

CTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTG  >  NC_000913/1649575‑1649857
                                                                                                                                               |                                                                                                                                           
cTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACggg                                                                                                                                        <  8:85038/149‑1 (MQ=255)
             ttCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATtaaata                                                                                                                            >  5:107276/1‑148 (MQ=255)
              tCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATtaaataa                                                                                                                           >  1:31588/1‑148 (MQ=255)
               cATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTGACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATtaaataaa                                                                                                                          <  2:170750/148‑1 (MQ=255)
                      gAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCa                                                                                                                   <  1:122692/148‑1 (MQ=255)
                       aGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCAtt                                                                                                                 >  6:133346/1‑149 (MQ=255)
                            gcgATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGttt                                                                                                             >  3:9353/1‑148 (MQ=255)
                              gATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACCAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCaa                                                                                                          <  1:82864/149‑1 (MQ=255)
                               aTAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAAt                                                                                                         <  7:77555/149‑1 (MQ=255)
                                   tGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGt                                                                                                               >  5:28612/1‑139 (MQ=255)
                                   tGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGt                                                                                                               <  6:28612/139‑1 (MQ=255)
                                           ttattaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGttt                                                                                             <  6:107276/149‑1 (MQ=255)
                                                  cAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAAt                                                                                                         >  1:107299/1‑130 (MQ=255)
                                                  cAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAAt                                                                                                         <  2:107299/130‑1 (MQ=255)
                                                       gTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGc                                                                                  >  7:121266/1‑148 (MQ=255)
                                                           aTACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCaaaa                                                                             <  2:31588/149‑1 (MQ=255)
                                                               gTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCAtt                                                                                                                 <  6:39058/109‑1 (MQ=255)
                                                               gTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCAtt                                                                                                                 >  5:39058/1‑109 (MQ=255)
                                                                  tGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGAtgtg                                                                      >  1:164645/1‑149 (MQ=255)
                                                                     aGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGca                                                                    >  1:132663/1‑148 (MQ=255)
                                                                       gttgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACaa                                                                 <  6:125171/149‑1 (MQ=255)
                                                                           ttttGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGaa                                                             >  5:92787/1‑149 (MQ=255)
                                                                            tttGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGAt                                                                                    <  2:159790/125‑1 (MQ=255)
                                                                            tttGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGAt                                                                                    >  1:159790/1‑125 (MQ=255)
                                                                                      gtgtCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCAttt                                                   <  6:92787/148‑1 (MQ=255)
                                                                                           aCTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCtt                                             <  2:164645/149‑1 (MQ=255)
                                                                                                 aGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGCCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCt                                        <  5:133346/148‑1 (MQ=255)
                                                                                                      cccAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGt                                  >  6:33545/1‑149 (MQ=255)
                                                                                                       ccAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTtgt                                                 >  7:83843/1‑133 (MQ=255)
                                                                                                       ccAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTtgt                                                 <  8:83843/133‑1 (MQ=255)
                                                                                                             tttACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCGCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCt                            <  8:121266/148‑1 (MQ=255)
                                                                                                                   gAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCa                                                         <  4:1181/113‑1 (MQ=255)
                                                                                                                   gAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCa                                                         >  3:1181/1‑113 (MQ=255)
                                                                                                                    aaGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCg                    <  5:33545/149‑1 (MQ=255)
                                                                                                                            aTTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTcc                     <  2:11384/140‑1 (MQ=255)
                                                                                                                            aTTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTcc                     >  1:11384/1‑140 (MQ=255)
                                                                                                                                  aaGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGgaga                >  5:148945/1‑139 (MQ=255)
                                                                                                                                  aaGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGgaga                <  6:148945/139‑1 (MQ=255)
                                                                                                                                  aaGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGa       >  4:50473/1‑148 (MQ=255)
                                                                                                                                   agagAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGata     >  7:40339/1‑149 (MQ=255)
                                                                                                                                      gaAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTg  <  4:9353/149‑1 (MQ=255)
                                                                                                                                      gaAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTg  >  3:84423/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                           
CTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTG  >  NC_000913/1649575‑1649857

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