Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,191,701 T→C 14.1% L375L (TTA→CTA) ‡ icd → e14 prophage; isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,191,7010TC14.1% 118.9 / 10.4 52L375L (TTA→CTA) icde14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base T (20/23);  new base C (4/3);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

ACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTT  >  CP009273/1191558‑1191849
                                                                                                                                               |                                                                                                                                                    
aCATCGGTTGCGAATGCGCTCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAtt                                                                                                                                                 <  1:502079/149‑1 (MQ=255)
     ggttaagAATGCGCCCTGTTTGAAGCAACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTAGTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaa                                                                                                                                            <  2:961063/143‑1 (MQ=255)
     ggcgACGAAGGCGCCCTGTTCGACGCCACCCACGGTACTGCGCCGAAATATGCCGCTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaa                                                                                                                                            <  2:186357/146‑1 (MQ=255)
      gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGGCGAAAAATGCAGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCCCTCCGCTGAGAGGATCGTTCGCCACATGGGTTGGACCGAAGCGGCTGACGTAATTGTTaaa                                                                                                                                           >  2:518984/1‑149 (MQ=255)
      gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGGTGATGCTGCGCCACATGGGTTGGACCGTAGCGGCTGGGTTAGTTGTTgaa                                                                                                                                           >  1:1183387/1‑149 (MQ=255)
      gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                                                                                                           >  2:280452/1‑149 (MQ=255)
      gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                                                                                                           >  1:823063/1‑149 (MQ=255)
      gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                                                                                                           <  1:1026520/149‑1 (MQ=255)
       cgACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                                                                                                          <  2:55891/148‑1 (MQ=38)
          cGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTa                                                                                                                                       <  2:69167/149‑1 (MQ=255)
           gAATGCGCCCTGTTTGAAGCCACCCACGGTACGGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTAt                                                                                                                                      <  2:1106414/149‑1 (MQ=255)
             aTGCGCCCTGTTTGAAGCCACCCACGGTACTGGGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                                                                                                                                    <  2:557926/149‑1 (MQ=255)
                   ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgcg                                                                                                                              <  1:682064/149‑1 (MQ=255)
                          gAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAAGCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGCACCGAAGCGGCTGTCTTAATTGTTAAAGGTATGGAAGGCGCAATCaa                                                                                                                       <  1:797966/149‑1 (MQ=255)
                           aaGCCACCCACGGTACTGCGCCGACATATGCCGGTCAGGACACAGTAAATCCTGGCTCTATTATTCTCGCCGCTGAGATGATGATGCGGCACATGGGATGGACCGACGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCaag                                                                                                                      >  2:587820/1‑148 (MQ=38)
                                acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaa                                                                                                                 <  2:1074076/149‑1 (MQ=255)
                                         aCTGCGCCGAAATATGCCGCTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGTCGCAATCAACGCGAAAACCGTAAc                                                                                                        <  1:918715/149‑1 (MQ=255)
                                           ggcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCAGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACATAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAcccc                                                                                                      >  2:18761/2‑148 (MQ=255)
                                               ccGAAATATGCCGGTCAGGACAAAGTTAATCCTGGCTCTATTATACTCTCCTCTGAGTGGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGa                                                                                                  <  1:711401/149‑1 (MQ=37)
                                               ccGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGCTCGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTACCTTATGa                                                                                                  >  2:751610/1‑149 (MQ=18)
                                                       tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGCATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGc                                                                                          >  1:1019743/1‑149 (MQ=255)
                                                          cGGTCAGGACAAAGTAAATCCTGGCTCTATTATTTTCCCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTc                                                                                       <  1:29543/149‑1 (MQ=255)
                                                           ggTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGATGCGTCAAATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt                                                                                      >  1:852845/1‑149 (MQ=17)
                                                             tcaggcatacgTAAGTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTga                                                                                    >  2:1134859/11‑149 (MQ=255)
                                                                 gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                                                <  1:1134163/149‑1 (MQ=255)
                                                                    aaactcaaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGatggatgg                                                                             <  2:221490/143‑1 (MQ=255)
                                                                    aaaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTTCGCCACATGGGTTGGTCCGAAGCGGCAGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAGCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGatggatgg                                                                             <  2:481477/149‑1 (MQ=255)
                                                                    aaaGGAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAATCCGTAACCTATGACTTCGAGCGTCTGatggatgg                                                                             <  1:688157/149‑1 (MQ=255)
                                                                      aGTAAGTCCTGGCTCTGTTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgcg                                                                           <  2:583801/149‑1 (MQ=11)
                                                                      aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgcg                                                                           <  2:707350/149‑1 (MQ=11)
                                                                       gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgcgc                                                                          >  1:1093492/1‑149 (MQ=11)
                                                                          aaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTaa                                                                       >  1:1179295/1‑149 (MQ=255)
                                                                           aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTaaa                                                                      >  1:1067930/1‑149 (MQ=255)
                                                                                 gCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTAGAACGTCTGATGGAAGGCGCTAAGctgctg                                                                >  1:368566/1‑149 (MQ=11)
                                                                                     tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGGTGGCGCTAAACTGCTGGAAt                                                            >  2:178273/1‑149 (MQ=255)
                                                                                     tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAAt                                                            <  1:110620/149‑1 (MQ=255)
                                                                                      attattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTCAACGGCTGAAATg                                                           >  2:528477/1‑149 (MQ=255)
                                                                                         attCTCTCCGCTGAGATGATGCCGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTc                                                        <  1:953126/149‑1 (MQ=255)
                                                                                              ctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGt                                                   >  2:935628/1‑149 (MQ=255)
                                                                                                          gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTACACTGATAAAATGATAAGAGTTTGGTGACGCGa                                       >  1:224367/1‑149 (MQ=18)
                                                                                                          gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGATTGGTGAggcgc                                       >  1:278578/1‑148 (MQ=255)
                                                                                                           atgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                                      <  1:351126/149‑1 (MQ=255)
                                                                                                             gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatca                                    >  2:835479/1‑149 (MQ=255)
                                                                                                               tgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                  >  1:1008762/1‑149 (MQ=255)
                                                                                                                  gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaa                               >  2:1150048/1‑149 (MQ=255)
                                                                                                                           gTTGGACCGAAGCGGCTGACTTAATTGTTAACGGTATGGAAGGCGCAACCACCGTGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                      <  2:1147048/149‑1 (MQ=255)
                                                                                                                            ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAAAGTTTGGTGACGCGATCATCGAAAACATGTAAt                     >  1:1059390/1‑149 (MQ=255)
                                                                                                                             tGGACCGAAGCGGCTGACTTAAGTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATg                    <  1:645801/149‑1 (MQ=255)
                                                                                                                                 ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGt                >  1:941991/1‑149 (MQ=255)
                                                                                                                                  ggtaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTa               <  2:924384/146‑1 (MQ=255)
                                                                                                                                           tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGGCGCGATCATCGAAAACATATAATGCCGTAGTTTGTTaa      >  2:473723/1‑149 (MQ=255)
                                                                                                                                           tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTaa      >  2:509939/1‑149 (MQ=255)
                                                                                                                                               ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAACCCGTAACCTCTGACTTCAAGCGTCTGATGGATGGCGCTAAGCAGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAttt  >  1:247936/1‑149 (MQ=18)
                                                                                                                                               ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGGGCGTCTGATGGATGGCGCTAAACTGCGGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGGAATGCCGTAGTTGGTTAAAttt  >  1:290196/1‑149 (MQ=255)
                                                                                                                                               ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAttt  >  2:548050/1‑149 (MQ=255)
                                                                                                                                               ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAttt  >  1:458122/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTT  >  CP009273/1191558‑1191849

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