Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,067,908 T→A 10.8% V963V (GTT→GTA 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
*CP0092732,067,9080TA10.8% 545.3 / 47.6 191V963V (GTT→GTAfluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (83/86);  new base A (11/9);  total (95/96)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 194 total aligned reads displayed.

TGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGG  >  CP009273/2067761‑2068053
                                                                                                                                                   |                                                                                                                                                 
tGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGtt                                                                                                                                                   >  2:1224129/1‑148 (MQ=255)
tGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGtt                                                                                                                                                   <  2:126272/148‑1 (MQ=255)
 gCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTa                                                                                                                                                  >  1:3202083/1‑148 (MQ=255)
 gCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTa                                                                                                                                                  <  1:3311585/148‑1 (MQ=255)
  cacaACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTAt                                                                                                                                                 <  2:2901331/148‑1 (MQ=255)
  cacaACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTAt                                                                                                                                                 <  2:2899233/148‑1 (MQ=255)
  cacaACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTAt                                                                                                                                                 <  1:1401586/148‑1 (MQ=255)
   acaacaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTAATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATc                                                                                                                                                <  2:209890/148‑1 (MQ=255)
      acaTGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGc                                                                                                                                             <  1:894482/148‑1 (MQ=255)
          gtgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCt                                                                                                                                         <  2:2948401/148‑1 (MQ=255)
          gtgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCt                                                                                                                                         >  1:3334867/1‑148 (MQ=255)
           tgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCtt                                                                                                                                        <  1:1013843/148‑1 (MQ=255)
           tgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCtt                                                                                                                                        <  1:2628319/148‑1 (MQ=255)
              gtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAg                                                                                                                                     <  1:311451/148‑1 (MQ=255)
               tgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGc                                                                                                                                    >  1:907178/1‑148 (MQ=255)
               tgcCGGTTTCCGTCTGGGCAGCCACAACGATAGGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGc                                                                                                                                    <  2:2405975/148‑1 (MQ=255)
                   ggTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTccc                                                                                                                                >  2:140064/1‑148 (MQ=255)
                    gTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCg                                                                                                                               >  1:3014877/1‑148 (MQ=255)
                     tttCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGACTTCTGTTATCCGCACCTTCAGCTCCCgg                                                                                                                              >  2:3546907/1‑148 (MQ=255)
                     tttCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTAGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCgg                                                                                                                              >  2:2358944/1‑148 (MQ=255)
                      ttCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCggg                                                                                                                             <  1:1137351/148‑1 (MQ=255)
                           tCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGat                                                                                                                        >  1:1072703/1‑148 (MQ=255)
                                gCAGCCACAACGATATGACCTTTGGCTAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCg                                                                                                                   <  2:275250/148‑1 (MQ=255)
                                gCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCg                                                                                                                   <  2:389677/148‑1 (MQ=255)
                                    ccACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGtgtg                                                                                                               >  1:2899295/1‑148 (MQ=255)
                                      acaACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTggg                                                                                                             <  2:2459233/148‑1 (MQ=255)
                                        aaCGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGa                                                                                                           <  1:2283180/148‑1 (MQ=255)
                                         aCGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGAc                                                                                                          <  1:1774674/148‑1 (MQ=255)
                                          cGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACt                                                                                                         <  1:1633918/148‑1 (MQ=255)
                                          cGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACt                                                                                                         >  2:14870/1‑148 (MQ=255)
                                          cGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACt                                                                                                         >  2:1176856/1‑148 (MQ=255)
                                            atatGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTc                                                                                                       >  1:1688790/1‑148 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                     >  2:967418/1‑148 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                     >  1:3504100/1‑148 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                     <  2:1841816/148‑1 (MQ=255)
                                               tGACCTTTGGCGAAGGCTCCTCATCCCGTGCCCCCCGGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGACCTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCAc                                                                                                    <  2:559040/148‑1 (MQ=255)
                                               tGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCAc                                                                                                    >  1:2361877/1‑148 (MQ=255)
                                                   cTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGCGATATGCGTGTGGGGACTTCCACTgca                                                                                                >  1:3095640/1‑148 (MQ=255)
                                                   cTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgca                                                                                                <  2:3508143/148‑1 (MQ=255)
                                                    tttGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcag                                                                                               >  2:2834803/1‑148 (MQ=255)
                                                          gAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCg                                                                                         >  1:2790700/1‑148 (MQ=255)
                                                             ggCACCTCCTCCCGTGCCCCCCTGCGTGACAGTGCAAAACCCAGTGTGAGTGAAGGACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGa                                                                                      <  2:2688070/148‑1 (MQ=255)
                                                             ggCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGa                                                                                      >  1:1312870/1‑148 (MQ=255)
                                                             ggCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGa                                                                                      >  1:517786/1‑148 (MQ=255)
                                                               cACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATg                                                                                    <  2:2032345/148‑1 (MQ=255)
                                                               cACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATg                                                                                    >  2:2079959/1‑148 (MQ=255)
                                                               cACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATg                                                                                    <  1:1460629/148‑1 (MQ=255)
                                                                aCCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGa                                                                                   >  2:3020923/1‑148 (MQ=255)
                                                                  cTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACg                                                                                 <  1:1010413/148‑1 (MQ=255)
                                                                  cTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACg                                                                                 <  1:3023837/148‑1 (MQ=255)
                                                                  cTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACg                                                                                 <  2:2155637/148‑1 (MQ=255)
                                                                  cTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACg                                                                                 >  1:473510/1‑148 (MQ=255)
                                                                  cTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATGACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACg                                                                                 <  1:610509/148‑1 (MQ=255)
                                                                   tCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGt                                                                                <  1:1987587/148‑1 (MQ=255)
                                                                   tCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGt                                                                                <  1:2183467/148‑1 (MQ=255)
                                                                    cATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGtt                                                                               <  2:2427564/148‑1 (MQ=255)
                                                                       cccGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtctc                                                                            >  2:2712929/1‑148 (MQ=255)
                                                                        ccGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGNTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtctct                                                                           <  1:1806638/148‑1 (MQ=255)
                                                                               ccccTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTcac                                                                    >  2:2048408/1‑148 (MQ=255)
                                                                               ccccTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTcac                                                                    >  2:2050449/1‑148 (MQ=255)
                                                                                  cTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGa                                                                 <  2:537285/148‑1 (MQ=255)
                                                                                    gCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAAt                                                               >  1:1195696/1‑148 (MQ=255)
                                                                                       tGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGc                                                            >  1:2489759/1‑148 (MQ=255)
                                                                                        gACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGca                                                           <  1:769613/148‑1 (MQ=255)
                                                                                        gACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTAGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGca                                                           >  2:3275468/1‑148 (MQ=255)
                                                                                          cAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGcaca                                                         >  1:3221701/1‑148 (MQ=255)
                                                                                                aaaaCACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTg                                                   <  1:2796396/148‑1 (MQ=255)
                                                                                                aaaaCACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTg                                                   <  1:240766/148‑1 (MQ=255)
                                                                                                  aaCACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGa                                                 <  2:2569436/148‑1 (MQ=255)
                                                                                                  aaCACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGa                                                 >  1:179344/1‑148 (MQ=255)
                                                                                                   acacAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGAc                                                <  2:2214975/148‑1 (MQ=255)
                                                                                                        gtgtGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCa                                           >  1:2850715/1‑148 (MQ=255)
                                                                                                            gagtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGcc                                       >  1:1351006/1‑148 (MQ=255)
                                                                                                               tgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGa                                    >  2:1540195/1‑148 (MQ=255)
                                                                                                                gaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGAc                                   <  2:363990/148‑1 (MQ=255)
                                                                                                                   ttACCGGTGAACTGGTGGGAACAGCCTACTGTGATACGCACCTTCAGCTCCCGGGGAGATATGCGGTGGGGGACTTCCACTGCAGGCAGCGGGATGACCATCGCACCATCACAGAATCGCACAGCACTGGACCTGCAGGCCGGACTgg                                >  2:2086295/1‑148 (MQ=255)
                                                                                                                        ggTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGcc                           <  2:1615658/148‑1 (MQ=255)
                                                                                                                            aaCTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCgtg                       <  2:70982/148‑1 (MQ=255)
                                                                                                                             aCTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCgtgt                      <  2:1664526/148‑1 (MQ=255)
                                                                                                                               tggtggGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTcc                    >  1:878492/1‑148 (MQ=255)
                                                                                                                               tggtggGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTcc                    >  1:878486/1‑148 (MQ=255)
                                                                                                                               tggtggGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTcc                    >  1:1246336/1‑148 (MQ=255)
                                                                                                                                 gtggGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCgg                  >  1:2600337/1‑148 (MQ=255)
                                                                                                                                   gggTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGa                <  1:3124946/148‑1 (MQ=255)
                                                                                                                                    ggTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaa               >  1:2902044/1‑148 (MQ=255)
                                                                                                                                      tACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaaaa             <  2:593201/148‑1 (MQ=255)
                                                                                                                                       aCAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAat            <  1:1353736/148‑1 (MQ=255)
                                                                                                                                        cAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAata           >  1:1223358/1‑148 (MQ=255)
                                                                                                                                          gCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATc         >  2:3133387/1‑148 (MQ=255)
                                                                                                                                           ccTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCa        >  1:3537796/1‑148 (MQ=255)
                                                                                                                                            cTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAc       <  1:2290554/148‑1 (MQ=255)
                                                                                                                                            cTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAc       <  2:2288487/148‑1 (MQ=255)
                                                                                                                                            cTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAc       <  2:641202/148‑1 (MQ=255)
                                                                                                                                             ttCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGACAATATCAcc      >  1:3058783/1‑148 (MQ=255)
                                                                                                                                             ttCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAcc      >  1:2109060/1‑148 (MQ=255)
                                                                                                                                             ttCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAcc      <  2:1070263/148‑1 (MQ=255)
                                                                                                                                              tCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAccc     <  2:1001890/148‑1 (MQ=255)
                                                                                                                                              tCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAccc     >  2:3378456/1‑148 (MQ=255)
                                                                                                                                               cTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCt    <  1:1534024/148‑1 (MQ=255)
                                                                                                                                                 gTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTgg  >  1:2311158/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGG  >  CP009273/2067761‑2068053

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