Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 1,195,717 A→C K199Q (AAG→CAG)  icd → isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,195,7170AC97.3% 103.1 / NA 37K199Q (AAG→CAG) icdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (15/21);  minor base G (0/1);  total (15/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 = 1.00e+00

TTATCTTCCGTGAAAACTCGGAAGACATTTATGCGGGTATCGAATGGAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTAAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAACATCATGAAGTTCACCGAAGGAGCGTTTAAAGACTGGGG  >  NC_000913_3_pae_pgi/1195571‑1195856
                                                                                                                                                  |                                                                                                                                           
ttATCTACCGTGAAAACTCGGAAGCCATTTATGCGGGTATCGAATGAAAAGCAGACTCTGCCGACCCCGAGAAAGTCATTAAATTCCTGCCTGAAGAGATGGCGGTGAATAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAg                                                                                                                                           <  2:319085/149‑1 (MQ=255)
         ttGAAAACTCGGAAGACATCCAAGCCGGTATCGAAATGAAAGCAGACCTTGCCGACGCCGAGAAAGTGATTAAATTCCTCCGTGAATAGATGGTGGTGAATAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCg                                                                                                                                  <  1:63919/148‑1 (MQ=255)
         gTGAAAACTCGGAAGACATTTATGGGGGTATCGAATGGAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTTGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCg                                                                                                                                  <  1:488709/149‑1 (MQ=255)
               aCTCGGAAGACATTGATGCGGGTATTGAATGGAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTTCTGCGTGAAGAGATGGGGGTTAATAAAATTCGCTTCCCGGAACATTGTTGTATCGGTATTCAGCCGTGTTCGgaagaa                                                                                                                            <  1:369326/149‑1 (MQ=255)
                 tCGGAAGACATTTATGCGGGTATCGAATGGAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAGATGCGGGTGACGCAAATTCGCTTCCCGGAACATTGTGGTGTCGGTGTGCGGGCGTGTTCGGAAGAAgg                                                                                                                          <  2:2810/149‑1 (MQ=255)
                          aTTTATGCGGGTATCGAATGGAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGAAAGAAGGCCCCAAACg                                                                                                                 >  1:239857/1‑149 (MQ=255)
                           tttATGCGGGTATCGAATGGAAAGCAGACTCTGCAGACGCCGAGAACGTGATTAAATTACTGCGTGACGAGATGGGGGTGAAGAGAATTCGATTCCCGGAACATTGTGGTATCGGTATTCAGCAGTGTTCGGAAGAAGGCACGAAACGt                                                                                                                >  1:722530/1‑149 (MQ=255)
                             taagCGGGTAACGAATGGAAAGCAAATTCTGCCGACGCCGAGAAAGTTATTAAATTCCTGAGTGAAGAGATGGGTGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGTACCAAACGTCt                                                                                                              <  2:777348/146‑1 (MQ=255)
                             tATGCGGGTATCGAATGGAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGTAAATTCGCTTCCCGGAATATTGTGGTATCGGTATTCAGCCGTGCTCGGAAGAAGGCACCAAACGTCt                                                                                                              <  1:393740/149‑1 (MQ=255)
                                  gggTATCGAATGGAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGGATTCAGCCGTGTTCGGAAGAAGGCACCAAACGCCTGGTTc                                                                                                         >  2:248505/1‑149 (MQ=255)
                                   ggTATCGAATGGAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGACACCAAACGTCTGGTTCg                                                                                                        <  2:647995/149‑1 (MQ=255)
                                             ggAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTTTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATc                                                                                              <  1:87965/149‑1 (MQ=255)
                                                aaGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAAATGGGGGTAAAGAAAATTCGCTTCCCTAAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCGAACGTCTGGTTCGTGCAGCGATCGaa                                                                                           >  1:698541/1‑149 (MQ=255)
                                                    aGACTCTGCCGACACCGAGAAAGTGATTAAATTCCTGCGAAACCATGGGGGGGTCTAGACAATCCGCCTCCCGGAACATTGTGGTCTCGGTATTCAGCCGTGTTCGAAGGAAGGCCCCAACCGTCTGGTGAGTGACGCGATCGaaagcg                                                                                       >  2:361910/1‑145 (MQ=255)
                                                       ctttgacgacgacgaGAAAGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAATAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCaa                                                                                    <  1:361910/143‑1 (MQ=255)
                                                                   cGAGAAAGGGATAAAATTCCTGCGTGAAGAGAGGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATc                                                                        <  2:450269/149‑1 (MQ=255)
                                                                    gagaAAGTGATTAAATTACTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCg                                                                       >  2:559492/1‑149 (MQ=255)
                                                                         aGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCgtgact                                                                  >  2:121325/1‑149 (MQ=255)
                                                                            gATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCgtgactctg                                                               <  1:391146/149‑1 (MQ=255)
                                                                               tAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCgtgactctgtga                                                            <  1:121325/149‑1 (MQ=255)
                                                                                    tCCCGCGTGAAGAGATGGGGGTTAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTTTTCAGCCGTGTTCGTAAGAAGGCACCTTTTGTCTGGTTTGTGCAGCGATCGAATACGCAATTGCTAACGATCGtgactctgtgactctg                                                       <  1:78647/149‑1 (MQ=255)
                                                                                       tGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTg                                                    >  2:510698/1‑149 (MQ=255)
                                                                                         cGTGAAGAGGTGGGGGTGAAGAAAATTCGCTCCCCGAAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGTCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGca                                                  <  1:115796/149‑1 (MQ=255)
                                                                                               gagaTGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAgg                                            >  1:719344/1‑149 (MQ=255)
                                                                                                    gggggTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAAca                                       >  2:49369/1‑149 (MQ=255)
                                                                                                     ggggTGAAGAAAATTCGTTTTCCTGAACATTGTTGTATCAGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAAcat                                      <  2:94644/149‑1 (MQ=255)
                                                                                                       ggTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGTTAACGATAGTGACTCTGTGACTCGGGTGCACAAAGGCAAcatca                                    <  1:196379/149‑1 (MQ=255)
                                                                                                        gTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAAcatcat                                   >  1:428342/1‑149 (MQ=255)
                                                                                                        gTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAAcatcat                                   >  2:663506/1‑149 (MQ=255)
                                                                                                           aagaaAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGCTGCACAAAGGCACCATCATGaa                                >  2:695759/1‑149 (MQ=255)
                                                                                                            agaaAATTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACACAGGCAACATCATGAAg                               <  1:695759/149‑1 (MQ=255)
                                                                                                                aaTTCGCTTCCCGGAACATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCGATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAACATCATGAAGTTCa                           >  2:430847/1‑149 (MQ=255)
                                                                                                                               aCATTGTGGTATCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAACATCATGAAGTTCACCGAAGGAGCGTTTa            <  2:719344/149‑1 (MQ=255)
                                                                                                                               aCATTGGTGTATCGGTATTCAGTCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGATTCTGTGACTCTGGTGCACAAAGGCAACATCATGAAGTTCACCGAAGGAGCGTTTa            <  1:125390/149‑1 (MQ=255)
                                                                                                                                        tagcGGTATTGAGCCGGGTGCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAACATCATGAAGTTCACCGAAGGTGCGTTTAAAGACTggg   <  2:662250/146‑1 (MQ=255)
                                                                                                                                         aTCGGTATTCAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGAGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAACATCATGAAGTTCACAGAAGGAGCGTTTACAGACTtggg  >  2:84975/1‑149 (MQ=255)
                                                                                                                                         aTCGGTATTCAGCCGTGTTAGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAACATCATGAAGTTCACCGAAGGAGCGTTTAAAGACTgggc  <  1:760562/149‑2 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
TTATCTTCCGTGAAAACTCGGAAGACATTTATGCGGGTATCGAATGGAAAGCAGACTCTGCCGACGCCGAGAAAGTGATTAAATTCCTGCGTGAAGAGATGGGGGTGAAGAAAATTCGCTTCCCGGAACATTGTGGTATCGGTATTAAGCCGTGTTCGGAAGAAGGCACCAAACGTCTGGTTCGTGCAGCGATCGAATACGCAATTGCTAACGATCGTGACTCTGTGACTCTGGTGCACAAAGGCAACATCATGAAGTTCACCGAAGGAGCGTTTAAAGACTGGGG  >  NC_000913_3_pae_pgi/1195571‑1195856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑