Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,066,503 C→A 7.8% S495* (TCA→TAA)  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,066,5030CA 7.8% 1282.6 / 53.2 340S495* (TCA→TAA) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base C (158/155);  new base A (11/16);  total (169/171)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

GGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTG  >  CP009273/2066355‑2066652
                                                                                                                                                    |                                                                                                                                                     
ggcACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCa                                                                                                                                                      >  1:5077658/1‑150 (MQ=255)
 gcACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATAAg                                                                                                                                                     <  1:4421658/150‑1 (MQ=255)
      cTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGt                                                                                                                                                <  1:899290/150‑1 (MQ=255)
      cTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGt                                                                                                                                                >  2:3097603/1‑150 (MQ=255)
         gCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATAAGGAAGTGGc                                                                                                                                             <  1:758340/150‑1 (MQ=255)
          cGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGca                                                                                                                                            <  1:293102/150‑1 (MQ=255)
           gggCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGcac                                                                                                                                           >  2:4443517/1‑150 (MQ=255)
            ggCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTAACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGcaca                                                                                                                                          <  1:1176725/150‑1 (MQ=255)
              caccaccacGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCAcact                                                                                                                                        <  1:4520089/150‑1 (MQ=255)
              caccaccacGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCAcact                                                                                                                                        <  1:4061971/150‑1 (MQ=255)
              caccaccacGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCAcact                                                                                                                                        >  1:812866/1‑150 (MQ=255)
               accaccacGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATAAGGAAGTGGCAcactc                                                                                                                                       <  2:576802/150‑1 (MQ=255)
                ccaccacGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCAcactca                                                                                                                                      >  1:2327306/1‑150 (MQ=255)
                    cacGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGt                                                                                                                                  <  1:3901975/150‑1 (MQ=255)
                    cacGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATAAGGAAGTGGCACACTCACTGt                                                                                                                                  <  1:1060223/150‑1 (MQ=255)
                     acGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATAAGGAAGTGGCACACTCACTGTc                                                                                                                                 <  2:5073392/150‑1 (MQ=255)
                       gCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAg                                                                                                                               <  2:50340/150‑1 (MQ=255)
                         tGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCa                                                                                                                             >  2:3864638/1‑150 (MQ=255)
                          gAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCaa                                                                                                                            >  2:3279379/1‑150 (MQ=255)
                          gAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCaa                                                                                                                            >  1:4199061/1‑150 (MQ=255)
                          gAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCaa                                                                                                                            <  1:4172151/150‑1 (MQ=255)
                              aaCGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAAcacc                                                                                                                        <  2:2480253/150‑1 (MQ=255)
                              aaCGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAAcacc                                                                                                                        <  2:1652059/150‑1 (MQ=255)
                                cGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAAcaccac                                                                                                                      >  1:4751606/1‑150 (MQ=255)
                                 ggCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACcaca                                                                                                                     <  2:1824560/150‑1 (MQ=255)
                                   cgcCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACt                                                                                                                   >  2:1287095/1‑150 (MQ=255)
                                   cgcCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACt                                                                                                                   >  2:4062758/1‑150 (MQ=255)
                                   cgcCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACt                                                                                                                   >  2:2537982/1‑150 (MQ=255)
                                   cgcCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACt                                                                                                                   <  1:2116711/150‑1 (MQ=255)
                                       aTACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCAcc                                                                                                               >  2:4356572/1‑150 (MQ=255)
                                        tACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCAccc                                                                                                              <  2:4257802/150‑1 (MQ=255)
                                         aCTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCa                                                                                                             >  2:3526221/1‑150 (MQ=255)
                                         aCTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATAAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCa                                                                                                             >  1:1608523/1‑150 (MQ=255)
                                           ttACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGa                                                                                                           <  1:1805043/150‑1 (MQ=255)
                                           ttACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGa                                                                                                           <  1:1294724/150‑1 (MQ=255)
                                             aCCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGaaa                                                                                                         >  1:4332505/1‑150 (MQ=255)
                                               ccTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGc                                                                                                       >  2:1345249/1‑150 (MQ=255)
                                                cTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCTTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGcc                                                                                                      >  1:1007775/1‑150 (MQ=255)
                                                  ttCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGAGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTACCGGCAGCGTGGAAAAATCAGGACGTGGCCCACTCACTGGCAGCACCACCACACTCACCCAGAAAGCCgg                                                                                                    >  2:2705443/1‑149 (MQ=255)
                                                      gggAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAAc                                                                                                <  2:1411526/150‑1 (MQ=255)
                                                       ggAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAAcc                                                                                               >  2:1937746/1‑150 (MQ=255)
                                                         aaGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCtg                                                                                             >  1:2509156/1‑150 (MQ=255)
                                                         aaGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCtg                                                                                             >  1:1953541/1‑150 (MQ=255)
                                                         aaGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCtg                                                                                             <  2:3967601/150‑1 (MQ=255)
                                                          aGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCtga                                                                                            >  2:4415528/1‑150 (MQ=255)
                                                              ggTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCtgaatga                                                                                        >  2:4625926/1‑150 (MQ=255)
                                                              ggTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCtgaatga                                                                                        <  2:2784160/150‑1 (MQ=255)
                                                              ggTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCtgaatga                                                                                        <  2:215141/150‑1 (MQ=255)
                                                                         aTACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCacgctga                                                                             >  2:4242249/1‑150 (MQ=255)
                                                                          tACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCacgctgac                                                                            <  1:4182507/150‑1 (MQ=255)
                                                                            cccTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCacgctgacgc                                                                          >  1:285291/1‑150 (MQ=255)
                                                                            cccTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCacgctgacgc                                                                          >  1:3259895/1‑150 (MQ=255)
                                                                             ccTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCCCCCAGAAAGCCGTCTACCTGAATGAAGGCacgctgacgct                                                                         <  2:2233424/150‑1 (MQ=255)
                                                                              cTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCacgctgacgctg                                                                        <  1:1998604/150‑1 (MQ=255)
                                                                                gACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCCGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGaa                                                                      >  2:3139895/1‑150 (MQ=255)
                                                                                gACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATAAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGaa                                                                      <  1:2238419/150‑1 (MQ=255)
                                                                                        gTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTa                                                              >  1:1598472/1‑150 (MQ=255)
                                                                                        gTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATAAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTa                                                              >  2:491670/1‑150 (MQ=255)
                                                                                         tGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAACGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTAc                                                             <  1:924832/150‑1 (MQ=255)
                                                                                         tGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTAc                                                             >  2:1748582/1‑150 (MQ=255)
                                                                                            aGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGt                                                          <  1:1710463/150‑1 (MQ=255)
                                                                                                gATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTcacc                                                      >  2:706948/1‑150 (MQ=255)
                                                                                                gATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTcacc                                                      >  2:1079202/1‑150 (MQ=255)
                                                                                                 aTGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATAAGGAAGTGGCACAGTCACTGTCAGCCACTCCCCACTCACGCAGATAGCGGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTcacca                                                     <  2:542798/150‑1 (MQ=255)
                                                                                                  tGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTcaccac                                                    <  1:3052658/150‑1 (MQ=255)
                                                                                                   gCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACg                                                   <  2:3394807/150‑1 (MQ=255)
                                                                                                   gCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACg                                                   >  2:665172/1‑150 (MQ=255)
                                                                                                     aCTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGa                                                 >  2:4863718/1‑150 (MQ=255)
                                                                                                          tGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCa                                            <  1:794056/150‑1 (MQ=255)
                                                                                                            cAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCAtt                                          <  1:3664933/150‑1 (MQ=255)
                                                                                                             aGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTg                                         <  1:4381398/150‑1 (MQ=255)
                                                                                                             aGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTg                                         <  2:4384260/150‑1 (MQ=255)
                                                                                                              gggAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGc                                        <  2:4670351/150‑1 (MQ=255)
                                                                                                               ggAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCt                                       >  2:2269963/1‑150 (MQ=255)
                                                                                                                gAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTc                                      <  1:2279206/150‑1 (MQ=255)
                                                                                                                   gCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgc                                   <  2:3197425/150‑1 (MQ=255)
                                                                                                                    ctctctCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcg                                  <  1:3097603/150‑1 (MQ=255)
                                                                                                                     tctctcACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAgcgc                                 >  2:1945252/1‑150 (MQ=255)
                                                                                                                      ctctcACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGcgcg                                <  2:3571738/150‑1 (MQ=255)
                                                                                                                      ctctcACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGcgcg                                <  2:1717834/150‑1 (MQ=255)
                                                                                                                        ctcACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGt                              <  2:4698819/150‑1 (MQ=255)
                                                                                                                         tcACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTa                             >  1:4709906/1‑150 (MQ=255)
                                                                                                                          cACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTAc                            >  2:1114687/1‑150 (MQ=255)
                                                                                                                          cACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTAc                            <  1:666185/150‑1 (MQ=255)
                                                                                                                           aCCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACa                           <  2:4507938/150‑1 (MQ=255)
                                                                                                                             cGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGc                         >  2:4759770/1‑150 (MQ=255)
                                                                                                                             cGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGc                         >  2:2882206/1‑150 (MQ=255)
                                                                                                                               gTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGccc                       >  2:1554338/1‑150 (MQ=255)
                                                                                                                                tAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCt                      >  1:1077955/1‑150 (MQ=255)
                                                                                                                                 aaCGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTg                     >  1:3361364/1‑150 (MQ=255)
                                                                                                                                    ggCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAg                  >  2:4142719/1‑150 (MQ=255)
                                                                                                                                      ctgCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCCGCACCACCACACTCACCCAGACAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCt                <  1:1761839/148‑1 (MQ=255)
                                                                                                                                          gTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGAcc            >  1:4827311/1‑150 (MQ=255)
                                                                                                                                          gTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGAcc            >  2:3677675/1‑150 (MQ=255)
                                                                                                                                          gTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGAcc            <  1:629267/150‑1 (MQ=255)
                                                                                                                                            ggAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCgg          >  2:4309905/1‑150 (MQ=255)
                                                                                                                                                aaaTCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGgcagc      <  2:2146134/150‑1 (MQ=255)
                                                                                                                                                 aaTCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGAGCGGTACAGCCCTGAAGCTGACCGgcagca     <  1:2538095/150‑1 (MQ=255)
                                                                                                                                                  aTCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCAc    <  1:3293815/150‑1 (MQ=255)
                                                                                                                                                    cAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACTCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACtg  <  1:1937746/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCACACTGGCGGGCACCACCACGCTGAATAACGGCGCCATACTTACCCTTTCCGGGAAGACGGTGAACAACGATACCCTGACCATCCGTGAAGGCGATGCACTCCTGCAGGGAGGCTCTCTCACCGGTAACGGCAGCGTGGAAAAATCAGGAAGTGGCACACTCACTGTCAGCAACACCACACTCACCCAGAAAGCCGTCAACCTGAATGAAGGCACGCTGACGCTGAACGACAGTACCGTCACCACGGATGTCATTGCTCAGCGCGGTACAGCCCTGAAGCTGACCGGCAGCACTG  >  CP009273/2066355‑2066652

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: