Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,073,151 T→C L538P (CTG→CCG)  flu → CP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,073,1510TC100.0% 89.3 / NA 25L538P (CTG→CCG) fluCP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/12);  total (13/12)

GGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGA  >  NC_000913/2073030‑2073284
                                                                                                                         |                                                                                                                                     
ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGGCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGgcagca                                                                                                                      <  2:450544/139‑1 (MQ=255)
ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGgcagca                                                                                                                      <  1:286310/139‑1 (MQ=255)
            gttaGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACg                                                                                                          <  2:94103/136‑1 (MQ=255)
             tCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACgg                                                                                                         <  2:184051/139‑1 (MQ=255)
                 cAACACCACACTCACCCAGAAAGCCGTCAACCTGGATGAAGGCACGCTGCCGCTGAACGACCGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGcc                                                                                                     >  2:462705/1‑139 (MQ=255)
                                                        aGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCAcc                                                              <  2:603945/139‑1 (MQ=255)
                                                              gctgacgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACg                                                                                         <  1:389545/106‑1 (MQ=255)
                                                              gctgacgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACg                                                                                         >  2:389545/1‑106 (MQ=255)
                                                              gctgacgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAat                                                        <  1:687522/139‑1 (MQ=255)
                                                                              aGTACCGTCACCTCGGATGTCATTGCCCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCa                                        >  2:171610/1‑139 (MQ=255)
                                                                                  ccGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCa                                        <  2:324675/135‑1 (MQ=255)
                                                                                  ccGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCa                                        >  1:324675/1‑135 (MQ=255)
                                                                                       accacGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGt                               >  2:625857/1‑139 (MQ=255)
                                                                                        ccacGGATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTc                              >  1:362118/1‑139 (MQ=255)
                                                                                            ggATGTCATTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCggtg                          >  1:93307/1‑139 (MQ=255)
                                                                                                   aTTGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATg                   >  2:401761/1‑139 (MQ=255)
                                                                                                     tGCTCAGCGCGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGAc                 >  1:589655/1‑139 (MQ=255)
                                                                                                           gcgcGGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGc           >  2:708377/1‑139 (MQ=255)
                                                                                                            cgcgGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCAtt                                                                                                  <  2:44617/51‑1 (MQ=255)
                                                                                                            cgcgGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCAtt                                                                                                  >  1:44617/1‑51 (MQ=255)
                                                                                                            cgcgGTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGcc          >  2:181152/1‑139 (MQ=255)
                                                                                                               ggTACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATg       <  2:589655/139‑1 (MQ=255)
                                                                                                                 tACAGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGcc     >  2:335746/1‑139 (MQ=255)
                                                                                                                    aGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGa  <  2:362118/139‑1 (MQ=255)
                                                                                                                    aGCCCCGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGa  <  1:335746/139‑1 (MQ=255)
                                                                                                                         |                                                                                                                                     
GGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGA  >  NC_000913/2073030‑2073284

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