Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,073,094 T→C L519P (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,0940TC100.0% 212.8 / NA 57L519P (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 (30/27);  total (30/27)

CGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGA  >  NC_000913/2072969‑2073228
                                                                                                                             |                                                                                                                                      
cGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGAc                                                                                                                           >  2:31782/1‑139 (MQ=255)
 gATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCAGAACGACa                                                                                                                          >  1:453833/1‑139 (MQ=255)
  aTGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGa                                                                                                                            >  1:422256/1‑136 (MQ=255)
  aTGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGa                                                                                                                            <  2:422256/136‑1 (MQ=255)
    gCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTa                                                                                                                       >  1:49396/1‑139 (MQ=255)
          cTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTca                                                                                                                 <  2:427892/139‑1 (MQ=255)
               gggAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACg                                                                                                            <  2:22060/139‑1 (MQ=255)
                  aGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGAt                                                                                                         >  1:135024/1‑139 (MQ=255)
                  aGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGAt                                                                                                         >  1:318193/1‑139 (MQ=255)
                    gCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGt                                                                                                       >  1:391663/1‑139 (MQ=255)
                      tctctcACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCa                                                                                                     >  2:338232/1‑139 (MQ=255)
                           cACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCt                                                                                                <  2:204124/139‑1 (MQ=255)
                                gTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           <  1:388045/139‑1 (MQ=255)
                                        aGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAg                                                                                   >  1:361758/1‑139 (MQ=255)
                                          cGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGcc                                                                                 >  2:332174/1‑139 (MQ=255)
                                           gTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGccc                                                                                <  2:64154/139‑1 (MQ=255)
                                              gAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCGCACTCACCCAGAAAGCCGTCAACCTGAATGAAGTCACGCCGACGCTGAACGACAGTACCGTCACCACGGACGTCATTGCTCAGCGCGGTACAGCCCTGa                                                                             <  2:49396/139‑1 (MQ=255)
                                                   aTCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTg                                                                        <  2:291809/139‑1 (MQ=255)
                                                   aTCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTg                                                                        <  2:135024/139‑1 (MQ=255)
                                                             ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  2:391663/139‑1 (MQ=255)
                                                             ggCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  1:31782/139‑1 (MQ=255)
                                                                    tcactGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCt                                                       >  1:301137/1‑139 (MQ=255)
                                                                      actGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGa                                                     >  1:143067/1‑139 (MQ=255)
                                                                         gTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACg                                                  >  1:203350/1‑139 (MQ=255)
                                                                            aGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgca                                                                 <  1:189735/121‑1 (MQ=255)
                                                                            aGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgca                                                                 >  2:189735/1‑121 (MQ=255)
                                                                            aGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGACAGCACTGTGCTGAACGGTg                                               >  2:332885/1‑139 (MQ=255)
                                                                             gCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACg                                                  >  1:195624/1‑135 (MQ=255)
                                                                             gCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACg                                                  <  2:195624/135‑1 (MQ=255)
                                                                              cAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              >  2:10745/1‑122 (MQ=255)
                                                                              cAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  1:10745/122‑1 (MQ=255)
                                                                                      cacTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTAc                                                                                     >  2:320376/1‑91 (MQ=255)
                                                                                      cacTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTAc                                                                                     <  1:320376/91‑1 (MQ=255)
                                                                                        cTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCATCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGCCCGGCAGCACTGTGCTGAACGGTGCCCTTTACCCCa                                   >  2:323510/1‑139 (MQ=255)
                                                                                         tCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCAc                                  <  2:318193/139‑1 (MQ=255)
                                                                                           aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTcacc                                                                                                               <  2:150871/60‑1 (MQ=255)
                                                                                           aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTcacc                                                                                                               >  1:150871/1‑60 (MQ=255)
                                                                                           aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                                                                           <  1:310610/80‑1 (MQ=255)
                                                                                           aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              >  2:148770/1‑109 (MQ=255)
                                                                                           aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  1:148770/109‑1 (MQ=255)
                                                                                           aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGa                                >  2:313660/1‑139 (MQ=255)
                                                                                           aCCCAGAAAGCCGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGAAGTCATTGCTCAgcg                                                                                           >  2:310610/1‑80 (MQ=255)
                                                                                                     ccGTCAACCTGAATGAAGGCACGCCGACGCTGAACGa                                                                                                                            >  2:274931/1‑37 (MQ=255)
                                                                                                     ccGTCAACCTGAATGAAGGCACGCCGACGCTGAACGa                                                                                                                            <  1:274931/37‑1 (MQ=255)
                                                                                                      cGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGaa                                                    >  2:32132/1‑108 (MQ=255)
                                                                                                      cGTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGaa                                                    <  1:32132/108‑1 (MQ=255)
                                                                                                       gTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCAc                                  >  2:154577/1‑125 (MQ=255)
                                                                                                       gTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCAc                                  <  1:154577/125‑1 (MQ=255)
                                                                                                       gTCAACCTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCg                    >  1:141393/1‑139 (MQ=255)
                                                                                                            ccTGAATGAAGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTcc               >  2:75030/1‑139 (MQ=255)
                                                                                                               gaatgaaGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGt            <  1:75030/139‑1 (MQ=255)
                                                                                                                aatgaaGGCACGCCGACGCTGAATGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGACAGCACTGTGCTGAACTGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTg           <  2:301137/139‑1 (MQ=255)
                                                                                                                    aaGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCAc       <  1:332174/139‑1 (MQ=255)
                                                                                                                     aGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGa                                                                                                          <  2:385631/39‑1 (MQ=255)
                                                                                                                     aGGCACGCCGACGCTGAACGACAGTACCGTCACCACGGa                                                                                                          >  1:385631/1‑39 (MQ=255)
                                                                                                                      ggCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              >  1:201176/1‑82 (MQ=255)
                                                                                                                      ggCACGCCGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagca                                                              <  2:201176/82‑1 (MQ=255)
                                                                                                                         acgccgacgctgAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGa  <  2:361758/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                                      
CGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTGTGCTGAACGGTGCCATTGACCCCACGAATGTCACTCTCGCCTCCGGTGCCACCTGGA  >  NC_000913/2072969‑2073228

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: