Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,073,225	T→C	W563R (TGG→CGG)	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_000913	2,073,225	0	T	C	100.0%	211.8 / NA	56	W563R (TGG→CGG)	flu	CP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/26);  total (30/26)

ACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGA  >  NC_000913/2073090‑2073360                                                                                                                                        |                                                                                                                                        acgctgacgctgAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                      <  2:930466/139‑1 (MQ=255) acgctgacgctgAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGa                                                                                                                                      <  2:1127752/139‑1 (MQ=255)  cgctgacgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGaa                                                                                                                                     >  1:715907/1‑139 (MQ=255)       acgctgaACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATcc                                                                                                                                <  2:1049434/139‑1 (MQ=255)                   aGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCa                                                                                                                    <  2:715907/139‑1 (MQ=255)                    gTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCAc                                                                                                                   <  2:1264917/139‑1 (MQ=255)                            accacGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGt                                                                                                           >  2:12251/1‑139 (MQ=255)                               acGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCa                                                                                                             <  1:1029378/134‑1 (MQ=255)                               acGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCa                                                                                                             >  2:1029378/1‑134 (MQ=255)                               acGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCgg                                                                                                        <  1:146027/139‑1 (MQ=255)                               acGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCgg                                                                                                        >  1:502462/1‑139 (MQ=255)                                      tCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                                                 <  1:1263054/139‑1 (MQ=255)                                        aTTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATg                                                                                               >  2:13898/1‑139 (MQ=255)                                            cTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCt                                                                                           >  1:2824/1‑139 (MQ=255)                                               aGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAg                                                                                        >  1:907436/1‑139 (MQ=255)                                                     gTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGc                                                                                  <  1:13898/139‑1 (MQ=255)                                                        cAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCgg                                                                               <  1:654149/139‑1 (MQ=255)                                                           cccTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACa                                                                            >  2:1227737/1‑139 (MQ=255)                                                             cTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGa                                                                          <  2:170688/139‑1 (MQ=255)                                                              tGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGAt                                                                         >  2:504561/1‑139 (MQ=255)                                                              tGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCAGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGCTGACCTCAGCCATGCCGGACAGAt                                                                         >  1:385806/1‑139 (MQ=255)                                                               gAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGAtt                                                                        <  1:1135593/139‑1 (MQ=255)                                                                aaGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTc                                                                       >  1:454905/1‑139 (MQ=255)                                                                  gCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCAt                                                                     >  1:523510/1‑139 (MQ=255)                                                                   cTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCAtt                                                                    <  1:245651/139‑1 (MQ=255)                                                                            agcaCTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCa                                                           <  2:907436/139‑1 (MQ=255)                                                                                          cGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAg                                             >  1:855523/1‑139 (MQ=255)                                                                                             tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAg                                                                                        >  1:367886/1‑93 (MQ=255)                                                                                             tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAg                                                                                        <  2:367886/93‑1 (MQ=255)                                                                                             tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                                          <  1:323079/123‑1 (MQ=255)                                                                                             tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                                          >  2:323079/1‑123 (MQ=255)                                                                                             tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGcaca                                                   <  2:1192781/130‑1 (MQ=255)                                                                                             tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGcaca                                                   >  1:1192781/1‑130 (MQ=255)                                                                                             tGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTc                                          >  2:802204/1‑139 (MQ=255)                                                                                                 aTTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTAc                                      <  1:504561/139‑1 (MQ=255)                                                                                                  ttGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                                                 >  1:831871/1‑79 (MQ=255)                                                                                                  ttGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGa                                                                                                 <  2:831871/79‑1 (MQ=255)                                                                                                     aCCCCACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCATACCGGc                                  <  2:916095/139‑1 (MQ=255)                                                                                                         cACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                                          <  1:424277/111‑1 (MQ=255)                                                                                                         cACGAATGTCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                                          >  2:424277/1‑111 (MQ=255)                                                                                                                 tCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGAc                                                                             >  2:939555/1‑84 (MQ=255)                                                                                                                 tCACTCTCGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGAc                                                                             <  1:939555/84‑1 (MQ=255)                                                                                                                       tcGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTGATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGaaaaa                >  2:700456/1‑139 (MQ=255)                                                                                                                       tcGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCg                                                       >  1:904354/1‑100 (MQ=255)                                                                                                                       tcGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCg                                                       <  2:904354/100‑1 (MQ=255)                                                                                                                       tcGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGaaaaa                >  2:722159/1‑139 (MQ=255)                                                                                                                       tcGCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGaaaaa                >  1:12073/1‑139 (MQ=255)                                                                                                                         gCCTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAAcc              >  2:567484/1‑139 (MQ=255)                                                                                                                          ccTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                                          <  2:827522/94‑1 (MQ=255)                                                                                                                          ccTCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCAc                                                          >  1:827522/1‑94 (MQ=255)                                                                                                                            tCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGAtt                                                                        >  2:629154/1‑78 (MQ=255)                                                                                                                            tCCGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGAtt                                                                        <  1:629154/78‑1 (MQ=255)                                                                                                                             ccGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCAtt                                                                    <  2:1006577/81‑1 (MQ=255)                                                                                                                             ccGGTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCAtt                                                                    >  1:1006577/1‑81 (MQ=255)                                                                                                                               ggTGCCACCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACg        >  2:865375/1‑139 (MQ=255)                                                                                                                                     aCCCGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGa  >  1:202846/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        ACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGAATATCCCCGATAACGCCACGGTGCAGTCGGTGGTGGATGACCTCAGCCATGCCGGACAGATTCATTTCACCTCCACCCGCACAGGGAAGTTCGTACCGGCAACCCTGAAAGTGAAAAACCTGAACGGACAGA  >  NC_000913/2073090‑2073360

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