Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,074,691 C→G 23.5% intergenic (+33/‑88) flu → / → yeeR CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter/CP4‑44 prophage; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,074,6910CG23.5% 74.7 / 21.1 34intergenic (+33/‑88)flu/yeeRCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter/CP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (14/12);  new base G (4/4);  total (18/16)
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.46e-01

GGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCT  >  NC_000913/2074556‑2074837
                                                                                                                                       |                                                                                                                                                  
ggAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCg                                                                                                                                             <  1:309728/143‑1 (MQ=255)
ggAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCg                                                                                                                                             >  2:309728/1‑143 (MQ=255)
 gAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCatgatg                                                                                                                                                 <  3:54823/138‑1 (MQ=255)
 gAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCatgatg                                                                                                                                                 >  4:54823/1‑138 (MQ=255)
       atCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatc                                                                                                                                                    <  3:804798/129‑1 (MQ=255)
       atCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatc                                                                                                                                                    >  4:804798/1‑129 (MQ=255)
                         gCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCgg                                                                                                                                            <  2:957187/119‑1 (MQ=255)
                         gCCCGTTATGCCCACAGCGTCAGTGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCgg                                                                                                                                            >  1:957187/1‑119 (MQ=255)
                          ccGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATGATGACCggg                                                                                                                                           <  3:376867/119‑1 (MQ=255)
                          caggTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGCCCCGGTCATGATGACCggg                                                                                                                                           >  4:376867/3‑119 (MQ=255)
                                             cAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatc                                                                                                                                                    >  2:676868/1‑91 (MQ=255)
                                             cAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatc                                                                                                                                                    <  1:676868/91‑1 (MQ=255)
                                                                   tataACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCAt                                                                                              >  4:10009/1‑123 (MQ=255)
                                                                   tataACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCAt                                                                                              <  3:10009/123‑1 (MQ=255)
                                                                                ccACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCggg                                                                                                                                           >  2:149270/1‑65 (MQ=255)
                                                                                ccACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCggg                                                                                                                                           <  1:149270/65‑1 (MQ=255)
                                                                                      tGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatc                                                                                                                                                    >  4:49965/1‑50 (MQ=255)
                                                                                      tGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGtcatc                                                                                                                                                    <  3:49965/50‑1 (MQ=255)
                                                                                      tGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATc                                                                                                                     >  4:1104452/1‑81 (MQ=255)
                                                                                      tGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATc                                                                                                                     <  3:1104452/81‑1 (MQ=255)
                                                                                       gAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGa                                                                                                                       >  2:467838/1‑78 (MQ=255)
                                                                                       gAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGa                                                                                                                       <  1:467838/78‑1 (MQ=255)
                                                                                        aaTGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGAcc                                                                                                                                        >  2:738789/1‑60 (MQ=255)
                                                                                        aaTGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGAcc                                                                                                                                        <  1:738789/60‑1 (MQ=255)
                                                                                              aCCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                                                                                                     <  4:952773/89‑1 (MQ=255)
                                                                                              aCCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                                                                                                     >  3:952773/1‑89 (MQ=255)
                                                                                                       cAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTAc                                                      >  3:46200/1‑127 (MQ=255)
                                                                                                       cAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTAc                                                      <  4:46200/127‑1 (MQ=255)
                                                                                                                    cctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTa                                                            >  4:942844/2‑108 (MQ=255)
                                                                                                                    cctctGTGATCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTa                                                            <  3:942844/107‑1 (MQ=255)
                                                                                                                     ctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGa                                    >  1:731125/1‑131 (MQ=255)
                                                                                                                     ctctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGa                                    <  2:731125/131‑1 (MQ=255)
                                                                                                                              cccGGTCATGATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACgg                                                                    <  1:681150/90‑1 (MQ=255)
                                                                                                                              cccGGTCATGATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACgg                                                                    >  2:681150/1‑90 (MQ=255)
                                                                                                                                cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATa                                                                                      >  3:545891/1‑70 (MQ=255)
                                                                                                                                cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATa                                                                                      <  4:545891/70‑1 (MQ=255)
                                                                                                                                cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGat                                                                               <  3:750575/77‑1 (MQ=255)
                                                                                                                                cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGat                                                                               >  4:750575/1‑77 (MQ=255)
                                                                                                                                cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcgcg                                        >  2:529682/1‑116 (MQ=255)
                                                                                                                                cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcgcg                                        <  1:529682/116‑1 (MQ=255)
                                                                                                                                cGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct        >  3:67046/1‑148 (MQ=255)
                                                                                                                                 ggTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttctt       >  3:105838/1‑148 (MQ=255)
                                                                                                                                   tcatcaTGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAAc                                                                                    >  4:260243/1‑69 (MQ=255)
                                                                                                                                   tcatcaTGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAAc                                                                                    <  3:260243/69‑1 (MQ=255)
                                                                                                                                      tcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCt  >  1:869447/1‑148 (MQ=255)
                                                                                                                                       catGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgc                              >  4:103672/1‑119 (MQ=255)
                                                                                                                                       catGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgc                              <  3:103672/119‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                  
GGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCT  >  NC_000913/2074556‑2074837

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: