Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,098,628 G→T 16.0% H9N (CAT→AAT)  csgD ← csgBAC operon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,098,6280GT16.0% 77.4 / 2.9 25H9N (CAT→AAT) csgDcsgBAC operon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base T (2/2);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

CAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTAAAATATAAAACTAATGGATTACATCTGATTTCAATCTAGCCATTACAAATCTTAA  >  CP009273/1098398‑1098828
                                                                                                                                                                                                                                      |                                                                                                                                                                                                        
caaTATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCtttttt                                                                                                                                                               >  1:41208/1‑274 (MQ=255)
                                          ccAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTAt                                                                                                                                                            <  1:257155/235‑1 (MQ=255)
                                               aatGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGAAGAAACACAGCTTTTTTAATTGATGGCAAACCTGACAACTGCCTCTAAAATAANAGCACCAGAAg                                                                                                               >  2:49859/1‑275 (MQ=255)
                                                                tctgctttgattatatGCAAAAGAATTATAGAGACTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCTTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGGGAGCTGTTGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGCCAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTCCa                                                                                              <  1:466853/263‑1 (MQ=255)
                                                                          tcatATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtgtg                                                                                    <  1:49861/275‑1 (MQ=255)
                                                                                    aagaaTAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATTAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTaaa                                                                                                                               >  2:395792/1‑222 (MQ=255)
                                                                                        ataataGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATTAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAg                                                                                                                           <  1:395805/222‑1 (MQ=255)
                                                                                                 ccTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATTAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATg                                                             >  1:139059/1‑275 (MQ=255)
                                                                                                      ctcgaTATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATTAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTa                                                         <  2:139053/271‑1 (MQ=255)
                                                                                                      agagATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGTGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCTTAAATGGGTa                                                         <  1:481264/274‑1 (MQ=255)
                                                                                                           tatCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTAAAatat                                                   <  2:41207/275‑1 (MQ=255)
                                                                                                                                                            cGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTaaaa                                                                                                                              >  2:225802/1‑151 (MQ=255)
                                                                                                                                                            cGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTACAATATAAAAATAATGGATTACATCTGATTTAAATCTAGCCATTACAAATCATaa  >  1:50517/1‑275 (MQ=255)
                                                                                                                                                                tGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGa                                                                                                                          <  1:225809/151‑1 (MQ=255)
                                                                                                                                                                  tttAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACtt                                                                                                                                                                                            >  2:51053/1‑83 (MQ=255)
                                                                                                                                                                      aGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCAtt                                                                                                                                                                                        <  1:51055/83‑1 (MQ=255)
                                                                                                                                                                              cAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACa                                                                                                       >  1:378030/1‑156 (MQ=255)
                                                                                                                                                                                  agagCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATg                                                                                                   <  2:378017/156‑1 (MQ=255)
                                                                                                                                                                                         tCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGt                                                                                                  >  1:366658/1‑150 (MQ=255)
                                                                                                                                                                                             cTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCa                                                                                              <  2:366645/150‑1 (MQ=255)
                                                                                                                                                                                                cAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCACGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGt                                                                                                  >  2:321993/1‑143 (MQ=255)
                                                                                                                                                                                                    gaagaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAATCACCAGAAGTACTGACAGATGTTGCa                                                                                              <  1:322003/143‑1 (MQ=255)
                                                                                                                                                                                                                               tGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTaaaa                                                                                                                              >  2:159615/1‑84 (MQ=255)
                                                                                                                                                                                                                                  ccATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTAAAATATaaaa                                               >  2:186294/1‑160 (MQ=255)
                                                                                                                                                                                                                                   cATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGa                                                                                                                          <  1:159621/84‑1 (MQ=255)
                                                                                                                                                                                                                                      gAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTAAAATATAAAACTaa                                           <  1:186300/160‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                        
CAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCAGCGATTGTTTAAGGTGCTGCAAGAGAGCTGTCGCCTGCAAAGAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTAAAATATAAAACTAATGGATTACATCTGATTTCAATCTAGCCATTACAAATCTTAA  >  CP009273/1098398‑1098828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: