Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,067,205 G→A 13.4% 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,2050GA13.4% 1174.0 / 140.6 381G729D (GGC→GAC) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base G (155/174);  new base A (18/33);  total (173/208)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 390 total aligned reads displayed.

CAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGG  >  CP009273/2067056‑2067354
                                                                                                                                                     |                                                                                                                                                     
cAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCACCCGGGACAGTGATGAGAGCTGGTATCTGCACAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAgg                                                                                                                                                       <  2:4847363/150‑1 (MQ=255)
 aGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGc                                                                                                                                                      <  1:4108271/150‑1 (MQ=255)
 aGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGc                                                                                                                                                      >  1:785280/1‑150 (MQ=255)
  ggCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACt                                                                                                                                                     >  1:2851370/1‑150 (MQ=255)
    cTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTcc                                                                                                                                                   <  2:1418747/150‑1 (MQ=255)
    cTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTcc                                                                                                                                                   >  1:4299862/1‑150 (MQ=255)
     tGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTccc                                                                                                                                                  >  2:1350627/1‑150 (MQ=255)
     tGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTccc                                                                                                                                                  >  1:2292967/1‑150 (MQ=255)
      gCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCg                                                                                                                                                 >  1:3110203/1‑150 (MQ=255)
       cAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGc                                                                                                                                                >  2:803601/1‑150 (MQ=255)
          gCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGc                                                                                                                                             >  1:580461/1‑150 (MQ=255)
            cGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCa                                                                                                                                           >  1:3874646/1‑150 (MQ=255)
              gTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATc                                                                                                                                         <  1:2404142/150‑1 (MQ=255)
               tgacTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCACGCAGCCATCa                                                                                                                                        <  1:4914080/147‑1 (MQ=255)
                gCCTTGACCGACTCCCTCCCCCGGGACCGGGATGAGAGCTGGTATCTGCGCAGTGCAATTGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCAGGCTGACACAGGCAATGGACGATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAg                                                                                                                                       <  2:3287433/150‑1 (MQ=255)
                   tttAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGAcc                                                                                                                                    >  2:1665769/1‑150 (MQ=255)
                    ttAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAAAGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCg                                                                                                                                   >  1:4714466/1‑150 (MQ=255)
                      aaCTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgt                                                                                                                                 <  1:3417739/150‑1 (MQ=255)
                       actactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCACAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtg                                                                                                                                <  1:2952286/150‑1 (MQ=255)
                       actactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCACAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtg                                                                                                                                <  1:4613007/150‑1 (MQ=255)
                        ctactCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGgtgt                                                                                                                               >  2:1451860/1‑150 (MQ=255)
                          actCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTaa                                                                                                                             <  1:4485555/150‑1 (MQ=255)
                           ctCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTaaa                                                                                                                            <  2:1761002/150‑1 (MQ=255)
                               cTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGt                                                                                                                        >  2:4326048/1‑150 (MQ=255)
                                 cAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGa                                                                                                                      >  1:82980/1‑150 (MQ=255)
                                 cAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGa                                                                                                                      <  2:4062020/150‑1 (MQ=255)
                                 cAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGa                                                                                                                      <  1:3146721/150‑1 (MQ=255)
                                    ccGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGaaaa                                                                                                                   <  1:4323961/150‑1 (MQ=255)
                                     cGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAaac                                                                                                                  <  2:706151/150‑1 (MQ=255)
                                     cGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAaac                                                                                                                  <  2:2145087/150‑1 (MQ=255)
                                         aCAGTGATGAGAGCTGGGATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAacaaca                                                                                                              <  2:1231127/150‑1 (MQ=255)
                                          cAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAg                                                                                                             >  1:3377098/1‑150 (MQ=255)
                                          cAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAg                                                                                                             >  2:4646420/1‑150 (MQ=255)
                                           aGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGc                                                                                                            >  1:218465/1‑150 (MQ=255)
                                            gTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcg                                                                                                           <  1:5046848/150‑1 (MQ=255)
                                             tgatgaGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgt                                                                                                          >  2:2010255/1‑150 (MQ=255)
                                              gatgaGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtc                                                                                                         <  2:1138713/150‑1 (MQ=255)
                                                  agagCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGcgtccgtc                                                                                                     <  1:3047940/150‑1 (MQ=255)
                                                     gCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCa                                                                                                  <  2:2931474/150‑1 (MQ=255)
                                                       tGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGc                                                                                                >  1:1252564/1‑150 (MQ=255)
                                                         gTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCAt                                                                                              >  1:5101868/1‑150 (MQ=255)
                                                         gTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCAt                                                                                              >  1:3587319/1‑150 (MQ=255)
                                                          tATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCAtt                                                                                             <  1:1927252/150‑1 (MQ=255)
                                                           aTCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTc                                                                                            <  2:4500514/150‑1 (MQ=255)
                                                           aTCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTc                                                                                            >  1:4010540/1‑150 (MQ=255)
                                                             cTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGAGTCCGTCTCAGCATTCAg                                                                                          <  2:3790482/150‑1 (MQ=255)
                                                               gcgcAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAggg                                                                                        <  1:633822/150‑1 (MQ=255)
                                                                      gAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCAt                                                                                 <  1:4963768/150‑1 (MQ=255)
                                                                         aaTGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCAtctc                                                                              <  2:3043257/150‑1 (MQ=255)
                                                                          aTGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCg                                                                             >  1:4856466/1‑150 (MQ=255)
                                                                          aTGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCg                                                                             >  2:1477730/1‑150 (MQ=255)
                                                                           tGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCgg                                                                            <  1:758559/150‑1 (MQ=255)
                                                                           tGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACGATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAACTGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCgc                                                                            <  1:308108/150‑2 (MQ=255)
                                                                             cTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTc                                                                          >  2:370453/1‑150 (MQ=255)
                                                                             cTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTc                                                                          <  1:1014734/150‑1 (MQ=255)
                                                                             cTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATGGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAACTGGTGAAAACAACAGCGGCCGTCTCAGCATTCAGGGCGGTCATCTCgggc                                                                          >  2:217447/1‑148 (MQ=255)
                                                                                       aGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACaa                                                                >  2:3428668/1‑150 (MQ=255)
                                                                                         aaGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATg                                                              <  1:1288282/150‑1 (MQ=255)
                                                                                         aaGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATg                                                              <  2:2344203/150‑1 (MQ=255)
                                                                                            tCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCg                                                           >  2:542287/1‑150 (MQ=255)
                                                                                              ccccTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGt                                                         <  1:4465913/150‑1 (MQ=255)
                                                                                               cccTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGTCTCCCGCAGCCATCAGTACGGTGTAAATGTTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTa                                                        <  2:4012853/150‑1 (MQ=255)
                                                                                               cccTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTa                                                        >  2:2661578/1‑150 (MQ=255)
                                                                                                 cTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTAtt                                                      <  1:2172933/150‑1 (MQ=255)
                                                                                                 cTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTAtt                                                      <  1:4416120/150‑1 (MQ=255)
                                                                                                 cTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTAtt                                                      <  2:4197537/150‑1 (MQ=255)
                                                                                                  tGTATGCCTCCATGCTGTCACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGTGCGGTCATCTCGGTCACGATAACAATGGCTGTATTg                                                     <  1:3918627/150‑1 (MQ=255)
                                                                                                  tGTATGCCTCCATGCTGACACAGGCAATGGACAATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGAATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTg                                                     >  1:87935/1‑150 (MQ=255)
                                                                                                   gTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGc                                                    >  2:498746/1‑150 (MQ=255)
                                                                                                    tATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTACATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGcc                                                   >  2:3231397/1‑150 (MQ=255)
                                                                                                    tATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGcc                                                   >  2:2323994/1‑150 (MQ=255)
                                                                                                     aTGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCTGCCATCAGACCGGTGTAGATGGTGAAAACAACAGCGGCCGTCTCAGCCGTCAGGGCGGTCATCGCGGTCACGGTCCCTATGGCGGTTGTGccc                                                  >  2:2717605/1‑150 (MQ=255)
                                                                                                     aTGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGccc                                                  <  1:534396/150‑1 (MQ=255)
                                                                                                      tGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTCTTGCCCg                                                 >  2:386150/1‑150 (MQ=255)
                                                                                                      tGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCg                                                 >  1:1750215/1‑150 (MQ=255)
                                                                                                       gCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGt                                                <  1:2885724/150‑1 (MQ=255)
                                                                                                           ccATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTgggg                                            <  1:1634542/150‑1 (MQ=255)
                                                                                                                cTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACg                                       >  2:1412538/1‑150 (MQ=255)
                                                                                                                cTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACg                                       >  1:1666895/1‑150 (MQ=255)
                                                                                                                 tGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGc                                      >  1:2627160/1‑150 (MQ=255)
                                                                                                                   acacAGGCAATGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCg                                    >  2:5116147/1‑150 (MQ=255)
                                                                                                                        ggCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAag                               >  1:2351572/1‑150 (MQ=255)
                                                                                                                            aTGGACTATGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACATCATCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACATTGGCGGTATTGCCCGTGGGGCCAAGCCGGAAagcagc                           <  1:4227624/150‑1 (MQ=255)
                                                                                                                             tGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCg                          >  2:3120204/1‑150 (MQ=255)
                                                                                                                               gACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGc                        <  2:4302810/150‑1 (MQ=255)
                                                                                                                               gACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGc                        <  1:1821356/150‑1 (MQ=255)
                                                                                                                                   aTGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCt                    >  2:128514/1‑150 (MQ=255)
                                                                                                                                    tGACCGGATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTa                   <  2:255628/150‑1 (MQ=255)
                                                                                                                                       ccGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATgg                <  1:1471798/150‑1 (MQ=255)
                                                                                                                                        cGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGAGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGa               <  2:940950/150‑1 (MQ=255)
                                                                                                                                         ggATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAt              <  2:1080376/150‑1 (MQ=255)
                                                                                                                                         ggATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAt              <  1:3266044/150‑1 (MQ=255)
                                                                                                                                          gATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAtt             >  2:3577730/1‑150 (MQ=255)
                                                                                                                                          gATTGTGGCAGACTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGAtt             <  1:1201263/150‑1 (MQ=255)
                                                                                                                                           aTTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTc            >  2:3669520/1‑150 (MQ=255)
                                                                                                                                              gtgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtc         >  2:1551461/1‑150 (MQ=255)
                                                                                                                                              gtgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtc         >  1:3529198/1‑150 (MQ=255)
                                                                                                                                              gtgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtc         <  1:2010255/150‑1 (MQ=255)
                                                                                                                                              gtgGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtc         <  1:1349386/150‑1 (MQ=255)
                                                                                                                                                 gCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTcgtccgt      <  2:2220996/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGCTGCAGGCCGGTGCCTTTAACTACTCCCTCAACCGGGACAGTGATGAGAGCTGGTATCTGCGCAGTGAAAATGCTTATCGTGCAGAAGTCCCCCTGTATGCCTCCATGCTGACACAGGCAATGGACTATGACCGGATTGTGGCAGGCTCCCGCAGCCATCAGACCGGTGTAAATGGTGAAAACAACAGCGTCCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGG  >  CP009273/2067056‑2067354

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