Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 273,525 T→C 17.5% E330G (GAG→GGG)  afuC ← CP4‑6 prophage; putative ferric transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273273,5250TC17.5% 38.0 / 3.9 23E330G (GAG→GGG) afuCCP4‑6 prophage; putative ferric transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (10/9);  new base C (2/2);  total (12/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 = 5.17e-01

CGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTG  >  CP009273/273400‑273668
                                                                                                                             |                                                                                                                                               
cgcagcgcaTCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGAcgc                                                                                                                          >  1:86257/1‑149 (MQ=255)
         tCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTg                                                                                                                 >  2:101958/1‑149 (MQ=255)
          cAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAACATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGa                                                                                                                <  1:113233/149‑1 (MQ=255)
            ggCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGAcc                                                                                                              >  2:145597/1‑149 (MQ=255)
                ttGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCa                                                                                                          <  2:223397/149‑1 (MQ=255)
                ttGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCa                                                                                                          <  2:82040/149‑1 (MQ=255)
                       aTGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCAGACTTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAatat                                                                                                   <  1:81675/149‑1 (MQ=255)
                                   gcgcCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTg                                                                                       <  2:61904/149‑1 (MQ=255)
                                    cgcCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCCCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGc                                                                                      >  1:152689/1‑149 (MQ=255)
                                    cgcCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCCCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGc                                                                                      >  1:152676/1‑149 (MQ=255)
                                                 caaCAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCa                                                                         <  2:145562/147‑1 (MQ=255)
                                                   aCAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCAcg                                                                       >  1:173138/1‑148 (MQ=255)
                                                          cacaAAACGTTATGCCGCATCCGCCAGAACAAACCTGCCGTACGGATGGATTTCAAGGTAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGGCCTGCCATAATATCGCCTGCCCGTGCCCTTCCACCGTCACTTcctcct                                                                >  1:219581/1‑144 (MQ=255)
                                                                            aTCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGGCCCATATAGGCGa                                              <  2:220259/149‑1 (MQ=255)
                                                                            aTCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGa                                              >  1:62059/1‑149 (MQ=255)
                                                                             tCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCCCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGAc                                             <  2:152676/149‑1 (MQ=255)
                                                                             tCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCCCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGAc                                             <  2:152689/149‑1 (MQ=255)
                                                                                gCCAGAACAAACATGCCGTACGGACGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATg                                          >  2:232093/1‑149 (MQ=255)
                                                                                 ccAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCCTATAGGCGACATgg                                         >  1:282288/1‑149 (MQ=255)
                                                                                 ccAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATgg                                         >  1:282315/1‑149 (MQ=255)
                                                                                                               cAAGATAATACTGCTCGCCGACGGCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTcgccg           <  1:101958/149‑1 (MQ=255)
                                                                                                                aaGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTcgccgc          <  1:153146/149‑1 (MQ=255)
                                                                                                                        aCTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATAAGGGCGACAGGTGGGATCAGGCGGCGCTGGCTCTCTTCGCcgcgctgactg  >  2:135605/1‑142 (MQ=255)
                                                                                                                             |                                                                                                                                               
CGCAGCGCATCAGGCATTGCCGGATGGCGGCACAAGCGCCTTATCCGGCCTACAAAACCACAAAACGTTATGCCGCATCCGCCAGAACAAACATGCCGTACGGATGGATTTCAAGATAATACTGCTCGCCGACGTCCGGTTGCAGACGCGTAGCGTTGACCTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTG  >  CP009273/273400‑273668

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