Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,067,205 G→A 5.2% G729D (GGC→GAC)  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,2050GA 5.2% 1555.1 / 44.3 422G729D (GGC→GAC) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base G (235/165);  new base A (16/6);  total (251/171)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 425 total aligned reads displayed.

CAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTG  >  CP009273/2067056‑2067353
                                                                                                                                                     |                                                                                                                                                    
cAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAgg                                                                                                                                                      <  2:2764965/150‑1 (MQ=255)
 aGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACAGAGGCAATGGACTATGACAGGATTGTGGCAGGc                                                                                                                                                     >  1:756240/1‑150 (MQ=255)
   gCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTc                                                                                                                                                   <  2:2385653/150‑1 (MQ=255)
     tGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTccc                                                                                                                                                 >  1:4572274/1‑150 (MQ=255)
     tGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTccc                                                                                                                                                 <  1:1411272/150‑1 (MQ=255)
       cAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGc                                                                                                                                               >  1:4037910/1‑150 (MQ=255)
       cAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGc                                                                                                                                               >  2:1783650/1‑150 (MQ=255)
         ggCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAg                                                                                                                                             >  2:1170283/1‑150 (MQ=255)
           ccGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGcc                                                                                                                                           >  2:3933135/1‑150 (MQ=255)
           ccGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGcc                                                                                                                                           <  1:3526413/150‑1 (MQ=255)
            cGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCa                                                                                                                                          >  2:4516934/1‑150 (MQ=255)
             ggTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCAt                                                                                                                                         <  2:1450836/150‑1 (MQ=255)
             ggTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCAt                                                                                                                                         >  2:3131678/1‑150 (MQ=255)
             ggTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCAt                                                                                                                                         >  2:705857/1‑150 (MQ=255)
               tGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCa                                                                                                                                       >  2:4724554/1‑150 (MQ=255)
                gCCTTTAACTACTCCCTCACCCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAg                                                                                                                                      <  2:1379549/150‑1 (MQ=255)
                gCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAg                                                                                                                                      >  2:836697/1‑150 (MQ=255)
                gCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAg                                                                                                                                      >  1:682685/1‑150 (MQ=255)
                gCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAg                                                                                                                                      >  2:3719231/1‑150 (MQ=255)
                  cTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGAc                                                                                                                                    <  1:2195199/150‑1 (MQ=255)
                       actcctCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtg                                                                                                                               <  1:2525575/150‑1 (MQ=255)
                       actactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtg                                                                                                                               <  1:2885102/150‑1 (MQ=255)
                        ctactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtgt                                                                                                                              <  1:3678593/150‑1 (MQ=255)
                             cccTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATg                                                                                                                         >  1:807101/1‑150 (MQ=255)
                             cccTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATg                                                                                                                         >  1:4497604/1‑150 (MQ=255)
                             cccTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATg                                                                                                                         >  1:3190826/1‑150 (MQ=255)
                                 cAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGa                                                                                                                     <  1:1840471/150‑1 (MQ=255)
                                  aaCCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGaa                                                                                                                    <  2:547063/150‑1 (MQ=255)
                                   aCCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGaaa                                                                                                                   >  1:1800681/1‑150 (MQ=255)
                                      gggACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAaaca                                                                                                                <  2:320135/150‑1 (MQ=255)
                                         aCAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAacaaca                                                                                                             <  2:4487465/150‑1 (MQ=255)
                                          cAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAg                                                                                                            <  1:2357955/150‑1 (MQ=255)
                                           aGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGc                                                                                                           <  2:3514170/150‑1 (MQ=255)
                                            gTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcg                                                                                                          <  1:3789313/150‑1 (MQ=255)
                                            gTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcg                                                                                                          <  2:2103413/150‑1 (MQ=255)
                                               atgaGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtcc                                                                                                       >  2:1317608/1‑150 (MQ=255)
                                                tgaGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtccg                                                                                                      <  2:2330147/150‑1 (MQ=255)
                                                  agagCTCGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtccgtc                                                                                                    >  2:4006152/1‑150 (MQ=255)
                                                   gagCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGtct                                                                                                   >  1:1177629/1‑150 (MQ=255)
                                                     gCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCa                                                                                                 >  2:1264980/1‑150 (MQ=255)
                                                      cTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAg                                                                                                >  2:597239/1‑150 (MQ=255)
                                                      cTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAg                                                                                                >  1:4661431/1‑150 (MQ=255)
                                                      cTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAg                                                                                                >  1:2280699/1‑150 (MQ=255)
                                                       tGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGc                                                                                               >  1:1273615/1‑150 (MQ=255)
                                                       tGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGc                                                                                               >  1:3313575/1‑150 (MQ=255)
                                                         gTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCAt                                                                                             >  1:4241306/1‑150 (MQ=255)
                                                          tATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCAtt                                                                                            >  1:300012/1‑150 (MQ=255)
                                                          tATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCAtt                                                                                            <  2:2156395/150‑1 (MQ=255)
                                                              tGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAgg                                                                                        <  2:1658913/150‑1 (MQ=255)
                                                              tGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAgg                                                                                        >  2:812457/1‑150 (MQ=255)
                                                                cgcAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGc                                                                                      >  2:1860052/1‑150 (MQ=255)
                                                                  cAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCgg                                                                                    >  2:3621614/1‑150 (MQ=255)
                                                                   aGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGt                                                                                   <  1:608045/150‑1 (MQ=255)
                                                                    gTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTc                                                                                  >  2:1224730/1‑150 (MQ=255)
                                                                       aaaaTGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTACAGGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCAtc                                                                               >  2:2228445/1‑150 (MQ=255)
                                                                       aaaaTGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCAtc                                                                               >  1:1048704/1‑150 (MQ=255)
                                                                             cTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTc                                                                         >  1:2448552/1‑150 (MQ=255)
                                                                              ttATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCa                                                                        <  1:3013334/150‑1 (MQ=255)
                                                                                     gCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAAc                                                                 <  1:3578046/150‑1 (MQ=255)
                                                                                      cAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACa                                                                <  2:3081635/150‑1 (MQ=255)
                                                                                      cAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACa                                                                >  1:956871/1‑150 (MQ=255)
                                                                                       aGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACaa                                                               >  1:39995/1‑150 (MQ=255)
                                                                                           gTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGc                                                           >  1:1735528/1‑150 (MQ=255)
                                                                                            tCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCg                                                          <  1:2575405/150‑1 (MQ=255)
                                                                                            tCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCg                                                          >  1:1684432/1‑150 (MQ=255)
                                                                                             cccccTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCgg                                                         >  1:2834491/1‑150 (MQ=255)
                                                                                              ccccTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGt                                                        >  1:2592269/1‑150 (MQ=255)
                                                                                                ccTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTAt                                                      >  1:2625131/1‑150 (MQ=255)
                                                                                                 cTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGAGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTAtt                                                     <  1:2743243/150‑1 (MQ=255)
                                                                                                  tGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTg                                                    >  2:2865148/1‑150 (MQ=255)
                                                                                                        ccTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATCACAATGGCGGTATTGCCCGTg                                              >  1:2066097/1‑150 (MQ=255)
                                                                                                        ccTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTg                                              >  2:1016272/1‑150 (MQ=255)
                                                                                                            cATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGc                                          >  1:3403683/1‑150 (MQ=255)
                                                                                                              tGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCa                                        >  1:2895169/1‑150 (MQ=255)
                                                                                                              tGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCa                                        <  2:4046665/150‑1 (MQ=255)
                                                                                                                 tGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGc                                     >  2:2460308/1‑150 (MQ=255)
                                                                                                                 tGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGc                                     >  2:444446/1‑150 (MQ=255)
                                                                                                                  gACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGcc                                    >  1:3253276/1‑150 (MQ=255)
                                                                                                                  gACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGcc                                    <  2:4304865/150‑1 (MQ=255)
                                                                                                                   acacAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCg                                   <  1:3131678/150‑1 (MQ=255)
                                                                                                                       aGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGaaa                               >  1:3118286/1‑150 (MQ=255)
                                                                                                                         gCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAagc                             >  2:1588154/1‑150 (MQ=255)
                                                                                                                         gCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAagc                             <  2:2431970/150‑1 (MQ=255)
                                                                                                                          cAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAagca                            <  2:4741903/150‑1 (MQ=255)
                                                                                                                           aaTGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAagcag                           >  2:2684197/1‑150 (MQ=255)
                                                                                                                                aCTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCa                      >  2:3439782/1‑150 (MQ=255)
                                                                                                                                 cTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGTAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAGTGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAg                     <  1:622877/150‑1 (MQ=255)
                                                                                                                                  tATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGc                    >  1:3480247/1‑150 (MQ=255)
                                                                                                                                  tATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGc                    <  2:682685/150‑1 (MQ=255)
                                                                                                                                   aTGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCt                   <  2:149637/150‑1 (MQ=255)
                                                                                                                                     gACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTAt                 <  2:3562258/150‑1 (MQ=255)
                                                                                                                                     cacaggATTGTGGCAGGCTCACGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGACCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTAt                 >  2:413574/5‑150 (MQ=255)
                                                                                                                                         ggATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAt             >  2:2023405/1‑150 (MQ=255)
                                                                                                                                         ggATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAt             >  2:174905/1‑150 (MQ=255)
                                                                                                                                             tgtgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgt         >  2:2558962/1‑150 (MQ=255)
                                                                                                                                              gtgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtc        >  1:2016410/1‑150 (MQ=255)
                                                                                                                                              gtgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtc        >  1:122389/1‑150 (MQ=255)
                                                                                                                                               tgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtcc       >  1:1618482/1‑150 (MQ=255)
                                                                                                                                               tgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtcc       >  1:413017/1‑150 (MQ=255)
                                                                                                                                                    ggCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTg  <  1:4342183/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTG  >  CP009273/2067056‑2067353

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: