Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,072,718 G→A 7.8% A394T (GCT→ACT)  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,072,7180GA 7.8% 160.0 / 3.5 51A394T (GCT→ACT) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base G (19/28);  new base A (2/2);  total (21/30)
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.97e-01

GATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCA  >  NC_000913/2072586‑2072862
                                                                                                                                    |                                                                                                                                                
gATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCg                                                                                                                                  >  1:512116/1‑149 (MQ=255)
 aTGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCgg                                                                                                                                 >  1:213022/1‑149 (MQ=255)
       tGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGAc                                                                                                                           <  1:570269/149‑1 (MQ=255)
        gACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTACCACTGACAGCGGCAcc                                                                                                                                    <  1:259055/139‑1 (MQ=255)
        gACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTACCACTGACAGCGGCAcc                                                                                                                                    >  2:259055/1‑139 (MQ=255)
            ggACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGaaa                                                                                                                      >  1:439413/1‑149 (MQ=255)
             gACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAg                                                                                                                     <  2:146241/149‑1 (MQ=255)
                  acaGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAg                                                                                                                     <  4:56200/144‑1 (MQ=255)
                  acaGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAg                                                                                                                     >  3:56200/1‑144 (MQ=255)
                  acaGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCAtt                                                                                                                >  3:60781/1‑149 (MQ=255)
                     gCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAg                                                                                                             >  2:451051/1‑149 (MQ=255)
                        aCTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTAt                                                                                                          >  1:58884/1‑149 (MQ=255)
                        aCTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTAt                                                                                                          >  1:79499/1‑149 (MQ=255)
                         cTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATc                                                                                                         <  2:480691/149‑1 (MQ=255)
                              aCCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAgg                                                                                                    <  2:58884/149‑1 (MQ=255)
                               cccGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGc                                                                                                   <  2:419485/149‑1 (MQ=255)
                               cccGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGc                                                                                                   <  3:365840/149‑1 (MQ=255)
                               cccGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGc                                                                                                   <  2:241233/149‑1 (MQ=255)
                                       gatgatGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGc                                                                                           >  1:248326/1‑149 (MQ=255)
                                           atgGCCGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGAt                                                                                       <  1:218102/149‑1 (MQ=255)
                                            tgGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATg                                                                                      >  2:28836/1‑149 (MQ=255)
                                             ggCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGc                                                                                     >  2:400011/1‑149 (MQ=255)
                                               cGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGccc                                                                                   <  1:34121/149‑1 (MQ=255)
                                                       tGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTg                                                                           <  2:213022/149‑1 (MQ=255)
                                                       tGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTg                                                                           <  1:400011/149‑1 (MQ=255)
                                                             tCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGaaaaa                                                                     <  2:493787/149‑1 (MQ=255)
                                                              ccGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAg                                                                    <  4:60781/149‑1 (MQ=255)
                                                               cGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAgg                                                                   <  1:143565/149‑1 (MQ=255)
                                                                gCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGc                                                                  <  3:65264/149‑1 (MQ=255)
                                                                 cAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCa                                                                 <  3:78111/149‑1 (MQ=255)
                                                                 cAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCa                                                                 <  1:147725/149‑1 (MQ=255)
                                                                          cACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGttcattca                                                        <  3:76460/149‑1 (MQ=255)
                                                                             cGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGc                                                     >  3:131024/1‑149 (MQ=255)
                                                                                 accaccGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGaa                                                 >  1:108878/1‑149 (MQ=255)
                                                                                 accaccGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGaa                                                 >  4:471952/1‑149 (MQ=255)
                                                                                    accGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGccc                                                                                   >  4:26942/1‑112 (MQ=255)
                                                                                    accGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGccc                                                                                   <  3:26942/112‑1 (MQ=255)
                                                                                    accGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGttca                                                            >  1:118805/1‑135 (MQ=255)
                                                                                    accGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGttca                                                            <  2:118805/135‑1 (MQ=255)
                                                                                          tCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGa                                        >  2:553624/1‑149 (MQ=255)
                                                                                            cATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCACTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATa                                      >  2:330058/1‑149 (MQ=255)
                                                                                                              aCTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGaa                                                 <  2:236658/120‑1 (MQ=255)
                                                                                                              aCTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGaa                                                 >  1:236658/1‑120 (MQ=255)
                                                                                                                     gCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATgg             <  2:429980/149‑1 (MQ=255)
                                                                                                                      ccGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACgg                    >  3:59676/1‑141 (MQ=255)
                                                                                                                      ccGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACgg                    <  4:59676/141‑1 (MQ=255)
                                                                                                                      ccGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGc            <  4:187102/149‑1 (MQ=255)
                                                                                                                       cGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCg           <  2:29902/149‑1 (MQ=255)
                                                                                                                           tCCGGTGCCACTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACt       <  1:330058/149‑1 (MQ=255)
                                                                                                                              ggTGCCGGTGGCAGTGGTATCCGGAGCTACGGATTGGCAGTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCTCGCTGAACGACGGTGATATGGCCACGGATACCACGGTAAATGGCGGACTGtt    <  2:345880/149‑1 (MQ=255)
                                                                                                                                tGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCa  <  2:79499/149‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                
GATGTGCTGACCGGACACACAGCCACTAATACCCGCGTGGATGATGGCGGAACGCTGGATGTCCGCAACGGTGGCACCGCCACCACCGTATCCATGGGAAATGGCGGTGTACTGCTGGCCGATTCCGGTGCCGCTGTCAGTGGTACCCGGAGCGACGGAAAGGCATTCAGTATCGGAGGCGGTCAGGCGGATGCCCTGATGCTGGAAAAAGGCAGTTCATTCACGCTGAACGCCGGTGATACGGCCACGGATACCACGGTAAATGGCGGACTGTTCA  >  NC_000913/2072586‑2072862

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: ‑