Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,067,908 T→A 41.2% 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,9080TA41.2% 124.4 / 252.3 192V963V (GTT→GTAfluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (51/61);  new base A (32/47);  total (84/108)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

TGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGC  >  CP009273/2067761‑2068055
                                                                                                                                                   |                                                                                                                                                   
tGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGACCTGGTGGGTACAGCCTTCTGtt                                                                                                                                                     >  2:271237/1‑148 (MQ=255)
tGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGtt                                                                                                                                                     <  1:659028/148‑1 (MQ=255)
  cacaACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAAt                                                                                                                                                   >  1:1615856/1‑148 (MQ=255)
         tgtgCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCAcc                                                                                                                                            <  1:14368/148‑1 (MQ=255)
         tgtgCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCAcc                                                                                                                                            >  1:198804/1‑148 (MQ=255)
          gtgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCt                                                                                                                                           <  1:1381843/148‑1 (MQ=255)
          gtgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCt                                                                                                                                           >  2:661830/1‑148 (MQ=255)
           tgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCtt                                                                                                                                          <  1:2989109/148‑1 (MQ=255)
           tgcgtgcCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCtt                                                                                                                                          <  1:480121/148‑1 (MQ=255)
                    gTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACATTTTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCg                                                                                                                                 <  1:3497214/148‑1 (MQ=255)
                          gTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGgaga                                                                                                                           >  1:1208984/1‑148 (MQ=255)
                           tCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGat                                                                                                                          <  1:2029276/148‑1 (MQ=255)
                           tCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGat                                                                                                                          <  2:783673/148‑1 (MQ=255)
                            cTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGata                                                                                                                         <  1:1576932/148‑1 (MQ=255)
                               ggCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGc                                                                                                                      >  2:2200618/1‑148 (MQ=255)
                               ggCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGc                                                                                                                      >  1:3048897/1‑148 (MQ=255)
                               ggCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGc                                                                                                                      >  1:685656/1‑148 (MQ=255)
                                gCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCg                                                                                                                     <  2:513691/148‑1 (MQ=255)
                                 cAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCgt                                                                                                                    <  2:3285065/148‑1 (MQ=255)
                                    ccACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGtgtg                                                                                                                 <  2:417339/148‑1 (MQ=255)
                                    ccACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGtgtg                                                                                                                 >  1:45831/1‑148 (MQ=255)
                                    ccACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGtgtg                                                                                                                 >  2:1682850/1‑148 (MQ=255)
                                     cacaACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTgg                                                                                                                <  1:1705663/148‑1 (MQ=255)
                                        aaCGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGa                                                                                                             <  1:3161714/148‑1 (MQ=255)
                                         aCGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGAc                                                                                                            >  1:395883/1‑148 (MQ=255)
                                           gATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACtt                                                                                                          >  2:1659674/1‑148 (MQ=255)
                                            atatGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTc                                                                                                         >  2:817306/1‑148 (MQ=255)
                                            atatGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTc                                                                                                         <  1:3528065/148‑1 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                       <  2:950932/148‑1 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                       <  2:2375376/148‑1 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                       <  2:586436/148‑1 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                       <  2:484527/148‑1 (MQ=255)
                                              atGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCa                                                                                                       >  1:2774513/1‑148 (MQ=255)
                                               tGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCAc                                                                                                      <  2:1902450/148‑1 (MQ=255)
                                               tGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCAc                                                                                                      <  2:786478/148‑1 (MQ=255)
                                                   cTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgca                                                                                                  <  1:2152347/148‑1 (MQ=255)
                                                     ttGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcagg                                                                                                <  1:320851/148‑1 (MQ=255)
                                                     ttGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcagg                                                                                                <  1:5269/148‑1 (MQ=255)
                                                      tGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTgcaggc                                                                                               <  2:1732526/148‑1 (MQ=255)
                                                            aGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCggg                                                                                         <  1:2782518/148‑1 (MQ=255)
                                                            aGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCggg                                                                                         <  2:2536504/148‑1 (MQ=255)
                                                             ggCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGa                                                                                        <  1:2228884/148‑1 (MQ=255)
                                                              gCACCTCTTCCCGTGGCCCCCTGCGTGACAGTGCAAAACACAGGGTGAGTGAATTACCGGTGAACTGGTGGGTACAGACTTCTGTTATCCGCACCTTGAGCTCCCGGGGAGATATGCGTGTGGGGTCTTCCACTGCAGGCAGCGGGAt                                                                                       <  2:3246614/148‑1 (MQ=255)
                                                               cACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATg                                                                                      <  2:2313390/148‑1 (MQ=255)
                                                               cACCTCATCCCGGGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATg                                                                                      <  1:366784/148‑1 (MQ=255)
                                                                 ccTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGAt                                                                                       >  2:1275421/1‑145 (MQ=255)
                                                                 ccTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGAc                                                                                    >  1:1143284/1‑148 (MQ=255)
                                                                 ccTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGAc                                                                                    <  1:2307562/148‑1 (MQ=255)
                                                                  cTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACg                                                                                   <  2:512255/148‑1 (MQ=255)
                                                                    cATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGtt                                                                                 >  2:1330254/1‑148 (MQ=255)
                                                                     aTCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtc                                                                                >  2:2233339/1‑148 (MQ=255)
                                                                     aTCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTtc                                                                                >  1:2703353/1‑148 (MQ=255)
                                                                          gTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTcc                                                                           <  2:1259894/148‑1 (MQ=255)
                                                                              cccccTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTca                                                                       <  2:190101/148‑1 (MQ=255)
                                                                                    gCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAAt                                                                 >  2:2141647/1‑148 (MQ=255)
                                                                                    gCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAAt                                                                 >  1:1904981/1‑148 (MQ=255)
                                                                                     cGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATg                                                                <  2:708045/148‑1 (MQ=255)
                                                                                             tGCAAAACACAGTGTGAGTTAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCa                                                        >  1:568182/1‑148 (MQ=255)
                                                                                             tGCAAAACACAGTGTGAGTGACTTACCGGTGTACTGGTGGGTCCAGCCTTCTGTTCTCCGCACCTTCAGCTCCCGGGGAGTTATGCGTGTGGGGAGTTCCCCTGCCGGCAGCGGGCTGACGTTCTCTCCCTCACCGAATGGCACATCa                                                        >  1:2883934/1‑148 (MQ=255)
                                                                                             tGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCa                                                        >  2:1489072/1‑148 (MQ=255)
                                                                                              gCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCAc                                                       >  2:2694350/1‑148 (MQ=255)
                                                                                                   acacAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGAc                                                  <  1:661830/148‑1 (MQ=255)
                                                                                                    cacaGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGAcc                                                 >  2:3219216/1‑148 (MQ=255)
                                                                                                      caGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTg                                               >  1:2276941/1‑148 (MQ=255)
                                                                                                       aGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGc                                              <  1:249222/148‑1 (MQ=255)
                                                                                                         tgtgAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAg                                            <  2:1024544/148‑1 (MQ=255)
                                                                                                            gagtgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCCGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGcc                                         >  2:1818988/1‑148 (MQ=255)
                                                                                                               tgaATTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGa                                      <  2:200559/148‑1 (MQ=255)
                                                                                                                  aTTACCGGTGAACTGGTGGGTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTg                                   <  2:3187477/148‑1 (MQ=255)
                                                                                                                   ttACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTgg                                  <  2:2076125/148‑1 (MQ=255)
                                                                                                                       cGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGc                              <  2:1208984/148‑1 (MQ=255)
                                                                                                                       cGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGc                              <  1:2633776/148‑1 (MQ=255)
                                                                                                                         gTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGccc                            >  1:1014461/1‑148 (MQ=255)
                                                                                                                         gTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGccc                            >  1:2558584/1‑148 (MQ=255)
                                                                                                                            aaCTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGGGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACCTCACTGGACCTGCAGGCCGGACTGGAAGCCCgtg                         >  2:2691368/1‑148 (MQ=255)
                                                                                                                                ggtggGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTACCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCg                     >  1:1609107/1‑148 (MQ=255)
                                                                                                                                 gtggGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCgg                    >  1:751492/1‑148 (MQ=255)
                                                                                                                                    ggTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaa                 >  1:1680557/1‑148 (MQ=255)
                                                                                                                                     gTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaaa                <  2:395883/148‑1 (MQ=255)
                                                                                                                                     gTACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaaa                <  2:396908/148‑1 (MQ=255)
                                                                                                                                      tACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaaaa               >  1:1603392/1‑148 (MQ=255)
                                                                                                                                      tACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaaaa               >  1:355706/1‑148 (MQ=255)
                                                                                                                                      tACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaaaa               <  1:1492946/148‑1 (MQ=255)
                                                                                                                                      tACAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGaaaa               >  1:1257794/1‑148 (MQ=255)
                                                                                                                                       aCAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAat              <  2:2958511/148‑1 (MQ=255)
                                                                                                                                        cAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGACGCCCGTGTCCGGGactata             >  1:959809/1‑143 (MQ=255)
                                                                                                                                        cAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAata             >  2:867984/1‑148 (MQ=255)
                                                                                                                                        cAGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAata             <  1:1349581/148‑1 (MQ=255)
                                                                                                                                         aGCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAatat            >  2:2666496/1‑148 (MQ=255)
                                                                                                                                          gCCTTCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATc           <  2:3449863/148‑1 (MQ=255)
                                                                                                                                           ccTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCa          >  2:521144/1‑148 (MQ=255)
                                                                                                                                             ttCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAcc        >  2:2008410/1‑148 (MQ=255)
                                                                                                                                             ttCTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAcc        <  2:2830484/148‑1 (MQ=255)
                                                                                                                                             ttCTGGGATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTTGGGACATCCACTGCAGGCAGCGGGATGTCGTTCTCTCCCACACAGAATGGCACTTCACTGGACCTGCAGTCCTGACTGGAAGCACGTGTCGGGGACCATATCAcc        >  2:2089493/1‑148 (MQ=255)
                                                                                                                                              tCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAccc       <  1:828995/148‑1 (MQ=255)
                                                                                                                                              tCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAccc       >  1:3351913/1‑148 (MQ=255)
                                                                                                                                              tCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCAccc       >  1:3168999/1‑148 (MQ=255)
                                                                                                                                               cTGTAATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCt      <  1:3266251/148‑1 (MQ=255)
                                                                                                                                                 gTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTgg    >  2:428819/1‑148 (MQ=255)
                                                                                                                                                  ttATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTggg   >  1:2822157/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCACAACATGTGCGTGCCGGTTTCCGTCTGGGCAGCCACAACGATATGACCTTTGGCGAAGGCACCTCATCCCGTGCCCCCCTGCGTGACAGTGCAAAACACAGTGTGAGTGAATTACCGGTGAACTGGTGGGTACAGCCTTCTGTTATCCGCACCTTCAGCTCCCGGGGAGATATGCGTGTGGGGACTTCCACTGCAGGCAGCGGGATGACGTTCTCTCCCTCACAGAATGGCACATCACTGGACCTGCAGGCCGGACTGGAAGCCCGTGTCCGGGAAAATATCACCCTGGGC  >  CP009273/2067761‑2068055

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: