Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,644,397 A→T 5.9% L417* (TTG→TAG)  yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,644,3970AT 5.9% 1010.6 / 37.3 286L417* (TTG→TAG) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (112/157);  new base T (6/11);  total (118/168)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 295 total aligned reads displayed.

TGTTAGCGCGAATATGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAGCGCCGGAAGT  >  CP009273/2644248‑2644546
                                                                                                                                                     |                                                                                                                                                     
tgtTAGCGCGAATATGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCa                                                                                                                                                       <  1:1261857/150‑1 (MQ=255)
tgtTAGCGCGAATATGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCa                                                                                                                                                       >  1:4258766/1‑150 (MQ=255)
     gcgcGAATATGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCt                                                                                                                                                  <  1:4321144/150‑1 (MQ=255)
       gcgAATATGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCttt                                                                                                                                                <  1:4923539/150‑1 (MQ=255)
         gAATATGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTAc                                                                                                                                              <  2:4509737/150‑1 (MQ=255)
         gAATATGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTAc                                                                                                                                              >  2:1516602/1‑150 (MQ=255)
           atatGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCg                                                                                                                                            <  1:86260/150‑1 (MQ=255)
            tatGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGt                                                                                                                                           >  2:3802537/1‑150 (MQ=255)
            tatGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGt                                                                                                                                           <  2:2555337/150‑1 (MQ=255)
              tGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCt                                                                                                                                         >  2:3514497/1‑150 (MQ=255)
               gCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTg                                                                                                                                        >  1:1823636/1‑150 (MQ=255)
                ccACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGc                                                                                                                                       >  1:2314205/1‑150 (MQ=255)
                  acaTACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCAt                                                                                                                                     >  1:1930605/1‑150 (MQ=255)
                    aTACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCa                                                                                                                                   <  1:384713/150‑1 (MQ=255)
                     tACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCAc                                                                                                                                  <  1:5171122/150‑1 (MQ=255)
                         ggTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTAACCGTCTGCATCACgcag                                                                                                                              >  2:2859930/1‑150 (MQ=255)
                           ttGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACgcagca                                                                                                                            >  2:4428149/1‑150 (MQ=255)
                            tGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCaa                                                                                                                           >  2:2034592/1‑150 (MQ=255)
                            tGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCaa                                                                                                                           >  1:665676/1‑150 (MQ=255)
                             gCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCaaa                                                                                                                          <  1:1589568/150‑1 (MQ=255)
                              ccGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAAc                                                                                                                         <  2:1611270/150‑1 (MQ=255)
                               cGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAAcc                                                                                                                        <  1:483561/150‑1 (MQ=255)
                                 cATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCTACGCTTTACCGTCTGCATCACACAGCAAACCAt                                                                                                                      <  1:3748379/150‑1 (MQ=255)
                                 cATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCAt                                                                                                                      <  1:1483275/150‑1 (MQ=255)
                                  aTTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCAtt                                                                                                                     >  1:3059232/1‑150 (MQ=255)
                                  aTTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCAtt                                                                                                                     >  1:4625428/1‑150 (MQ=255)
                                   ttGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTg                                                                                                                    <  2:891134/150‑1 (MQ=255)
                                     gCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAg                                                                                                                  <  2:3160766/150‑1 (MQ=255)
                                       tATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGa                                                                                                                >  1:3893406/1‑150 (MQ=255)
                                        aTCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGAt                                                                                                               <  2:1134661/150‑1 (MQ=255)
                                        aTCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGAt                                                                                                               >  2:2318161/1‑150 (MQ=255)
                                        aTCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGAt                                                                                                               <  2:1965589/150‑1 (MQ=255)
                                           gAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTAc                                                                                                            >  1:1446071/1‑150 (MQ=255)
                                                gCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGttt                                                                                                       >  1:1339745/1‑150 (MQ=255)
                                                   aGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCAc                                                                                                    >  2:74432/1‑150 (MQ=255)
                                                     gTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCt                                                                                                  <  1:2763143/150‑1 (MQ=255)
                                                       aaaGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTGTCACCTgg                                                                                                >  2:781359/1‑150 (MQ=255)
                                                           tGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGa                                                                                            >  2:3122965/1‑150 (MQ=255)
                                                              tAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCTACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATaa                                                                                         >  1:4962477/1‑150 (MQ=255)
                                                              tAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATaa                                                                                         <  1:3033784/150‑1 (MQ=255)
                                                                 agGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCTACGCTTTACAGTCTGCAGCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGa                                                                                      >  2:3872466/1‑150 (MQ=255)
                                                                 agGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGa                                                                                      <  2:4862240/150‑1 (MQ=255)
                                                                    ccATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATcg                                                                                   <  2:5200390/150‑1 (MQ=255)
                                                                    ccATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATcg                                                                                   >  2:3090512/1‑150 (MQ=255)
                                                                     cATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATcgc                                                                                  >  2:3414228/1‑150 (MQ=255)
                                                                      aTTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATcgcg                                                                                 >  2:3016042/1‑150 (MQ=255)
                                                                      aTTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATcgcg                                                                                 >  1:3983169/1‑150 (MQ=255)
                                                                          tcCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTggg                                                                             <  1:2978305/150‑1 (MQ=255)
                                                                          tcCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTggg                                                                             <  2:3131512/150‑1 (MQ=255)
                                                                            cGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGcc                                                                           <  1:2794729/150‑1 (MQ=255)
                                                                              gCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCTACGCTTTACAGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCaa                                                                         <  1:357937/150‑1 (MQ=255)
                                                                                 gACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACa                                                                      <  1:4723641/150‑1 (MQ=255)
                                                                                   cTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATg                                                                    <  2:3258240/150‑1 (MQ=255)
                                                                                         acgCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAAc                                                              <  2:4320050/150‑1 (MQ=255)
                                                                                          cgCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACt                                                             >  1:5002704/1‑150 (MQ=255)
                                                                                            cTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGt                                                           >  2:3504587/1‑150 (MQ=255)
                                                                                             tCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGtt                                                          >  1:3087483/1‑150 (MQ=255)
                                                                                              ccTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGttt                                                         <  1:3337733/150‑1 (MQ=255)
                                                                                              ccTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGttt                                                         <  1:235166/150‑1 (MQ=255)
                                                                                                 gAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCTACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCt                                                      >  2:4512673/1‑150 (MQ=255)
                                                                                                 gAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCt                                                      >  1:2916976/1‑150 (MQ=255)
                                                                                                   gTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCtat                                                    <  2:4079004/150‑1 (MQ=255)
                                                                                                       cTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGcc                                                >  2:3108985/1‑150 (MQ=255)
                                                                                                       cTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGcc                                                >  2:3109089/1‑150 (MQ=255)
                                                                                                         gCCCATCGGGTTTACTCACGTCTTACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACTCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTg                                              <  2:1271258/150‑1 (MQ=255)
                                                                                                         gCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTg                                              <  2:639517/150‑1 (MQ=255)
                                                                                                         gCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTg                                              >  2:4720912/1‑150 (MQ=255)
                                                                                                          cccATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTgg                                             <  2:1199871/150‑1 (MQ=255)
                                                                                                            cATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGgcg                                           >  1:55515/1‑150 (MQ=255)
                                                                                                            cATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGgcg                                           <  2:673724/150‑1 (MQ=255)
                                                                                                               cGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGgcgcgc                                        >  2:2004434/1‑150 (MQ=255)
                                                                                                                 ggTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCa                                      <  1:747453/150‑1 (MQ=255)
                                                                                                                 ggTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCa                                      <  1:319092/150‑1 (MQ=255)
                                                                                                                 ggTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCa                                      >  1:4382636/1‑150 (MQ=255)
                                                                                                                 ggTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCa                                      <  2:4043451/150‑1 (MQ=255)
                                                                                                                  gTTTAATCACGTCTAACTTGATGGGTTGATTGGGCTACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCAACTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACGGTTTGCTATAGCCTGGCGCGCCAg                                     >  2:2475867/1‑150 (MQ=255)
                                                                                                                   tttAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGc                                    >  1:3741411/1‑150 (MQ=255)
                                                                                                                   tttAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGc                                    <  1:1305658/150‑1 (MQ=255)
                                                                                                                    ttAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCa                                   <  1:2687815/150‑1 (MQ=255)
                                                                                                                       aTCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATg                                <  1:3454460/150‑1 (MQ=255)
                                                                                                                        tCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGt                               >  2:3739534/1‑150 (MQ=255)
                                                                                                                           cGTCTAACGTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTaa                            <  2:58355/150‑1 (MQ=255)
                                                                                                                              cTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAAtt                         <  2:2582633/150‑1 (MQ=255)
                                                                                                                                aaCTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCt                       <  2:3364147/150‑1 (MQ=255)
                                                                                                                                 aCTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTg                      <  2:1039133/150‑1 (MQ=255)
                                                                                                                                 aCTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCAAGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTg                      >  1:681461/1‑150 (MQ=255)
                                                                                                                                     gATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTaa                  >  2:3809818/1‑150 (MQ=255)
                                                                                                                                       tGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGt                <  2:1660075/150‑1 (MQ=255)
                                                                                                                                        gggTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTc               <  1:261739/150‑1 (MQ=255)
                                                                                                                                          gTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCa             <  1:1034652/150‑1 (MQ=255)
                                                                                                                                          gTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCa             <  2:5035044/150‑1 (MQ=255)
                                                                                                                                           ttGATTGGGCTACGCTTTACAGTCTGCATCACGCAGCAAACCATTGAGGATTACGGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAATAACTGTTTGCTAGAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAg            <  2:4218173/150‑1 (MQ=255)
                                                                                                                                            tGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAgc           <  2:1823636/150‑1 (MQ=255)
                                                                                                                                            tGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAgc           <  1:619179/150‑1 (MQ=255)
                                                                                                                                             gATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGGAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAgcg          >  1:4566606/1‑150 (MQ=255)
                                                                                                                                             gATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAgcg          >  1:420175/1‑150 (MQ=255)
                                                                                                                                                 gggCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAGCGCCgg      <  1:514003/150‑1 (MQ=255)
                                                                                                                                                  ggCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAGCGCCGGa     <  1:2333092/150‑1 (MQ=255)
                                                                                                                                                  ggCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAGCGCCGGa     <  1:2805693/150‑1 (MQ=255)
                                                                                                                                                     aaCGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAGCGCCGGAAGt  >  1:1008321/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTAGCGCGAATATGCCACATACCGGTTGCCGCATTGCTATCGAGTGGCCAGGTAAAGTGGTAGAGGCCATTCTCCGGCTGACTAACGACGCTCCTGAGTACCTGCCCATCGGGTTTAATCACGTCTAACTTGATGGGTTGATTGGGCAACGCTTTACCGTCTGCATCACGCAGCAAACCATTGAGGATTACCGTTTCACCTGGGCGATAAAGATCGCGTGGGCCAAACATGAAAAACTGTTTGCTATAGCCTGGCGCGCCAGCAATGTTAAATTCTGCTAAGTCCAGCGCCGGAAGT  >  CP009273/2644248‑2644546

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: