Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,643,286 C→A 17.4% A787A (GCG→GCT 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,643,2860CA17.4% 1042.5 / 158.8 343A787A (GCG→GCTyfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (166/117);  new base A (30/30);  total (196/147)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 347 total aligned reads displayed.

GAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCA  >  CP009273/2643138‑2643435
                                                                                                                                                    |                                                                                                                                                     
 agaTTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGc                                                                                                                                                     >  2:3281832/1‑150 (MQ=255)
  gaTTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCa                                                                                                                                                    >  1:964750/1‑150 (MQ=255)
    tttATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATa                                                                                                                                                  <  1:3374521/150‑1 (MQ=255)
     ttATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACAGCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATaa                                                                                                                                                 >  1:4377415/1‑150 (MQ=255)
       aTCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAAcc                                                                                                                                               <  2:3316736/150‑1 (MQ=255)
        tCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCt                                                                                                                                              >  2:597186/1‑150 (MQ=255)
        tCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCt                                                                                                                                              >  1:3296383/1‑150 (MQ=255)
          gCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCttt                                                                                                                                            >  1:2992930/1‑150 (MQ=255)
               gACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACccg                                                                                                                                       <  1:1730445/150‑1 (MQ=255)
                   tACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGt                                                                                                                                   <  2:82515/150‑1 (MQ=255)
                   tACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGt                                                                                                                                   <  1:608518/150‑1 (MQ=255)
                    accaccAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGtt                                                                                                                                  <  1:1018224/150‑1 (MQ=255)
                    accaccAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGtt                                                                                                                                  >  1:2924537/1‑150 (MQ=255)
                    accaccAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGtt                                                                                                                                  >  1:1499291/1‑150 (MQ=255)
                     ccaccaGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTg                                                                                                                                 <  1:4788305/150‑1 (MQ=255)
                     ccaccaGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTg                                                                                                                                 >  2:2028646/1‑150 (MQ=255)
                        ccaGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGgcg                                                                                                                              >  1:2778480/1‑150 (MQ=255)
                          aGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGcgcg                                                                                                                            <  1:3744436/150‑1 (MQ=255)
                          aGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGcgcg                                                                                                                            >  1:3840951/1‑150 (MQ=255)
                          aGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGcgcg                                                                                                                            <  2:4416909/150‑1 (MQ=255)
                            cGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCgcggc                                                                                                                          <  2:2425344/150‑1 (MQ=255)
                            cGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCgcggc                                                                                                                          <  2:3416303/150‑1 (MQ=255)
                              tACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGa                                                                                                                        >  1:5109/1‑150 (MQ=255)
                               aCTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGat                                                                                                                       <  2:2509930/150‑1 (MQ=255)
                                cTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGata                                                                                                                      <  2:2965523/150‑1 (MQ=255)
                                    aaTACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCa                                                                                                                  <  1:1868036/150‑1 (MQ=255)
                                      tACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAAc                                                                                                                <  1:178184/150‑1 (MQ=255)
                                       aCAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACt                                                                                                               <  1:101044/150‑1 (MQ=255)
                                        cAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACtt                                                                                                              >  1:352871/1‑150 (MQ=255)
                                        cAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACtt                                                                                                              >  1:3607281/1‑150 (MQ=255)
                                        cAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACtt                                                                                                              >  1:2125924/1‑150 (MQ=255)
                                        cAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACtt                                                                                                              >  1:1385882/1‑150 (MQ=255)
                                         aGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACttc                                                                                                             >  2:90558/1‑149 (MQ=255)
                                         aGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTa                                                                                                             <  2:1308162/150‑1 (MQ=255)
                                         aGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTa                                                                                                             >  2:555548/1‑150 (MQ=255)
                                           aTCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAAt                                                                                                           <  1:2003186/150‑1 (MQ=255)
                                            tCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTCCCCGCCGTTGGCGCGGCGATATGCAACTTAATg                                                                                                          >  2:306238/1‑150 (MQ=255)
                                                gACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGgtgt                                                                                                      <  2:1254314/150‑1 (MQ=255)
                                                gACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGgtgt                                                                                                      >  1:748508/1‑150 (MQ=255)
                                                 aCGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTc                                                                                                     >  2:3914486/1‑150 (MQ=255)
                                                 aCGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTc                                                                                                     <  2:1826350/150‑1 (MQ=255)
                                                  cGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCg                                                                                                    >  1:2625478/1‑150 (MQ=255)
                                                     ttCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCa                                                                                                 <  2:4467124/150‑1 (MQ=255)
                                                     ttCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCa                                                                                                 <  1:2428280/150‑1 (MQ=255)
                                                     ttCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCa                                                                                                 <  1:4692237/150‑1 (MQ=255)
                                                         aGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGc                                                                                             >  1:2553996/1‑150 (MQ=255)
                                                         aGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGc                                                                                             <  2:1500714/150‑1 (MQ=255)
                                                            ttttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGAt                                                                                          >  1:2180090/1‑150 (MQ=255)
                                                             tttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATa                                                                                         >  2:1713586/1‑150 (MQ=255)
                                                                   gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGc                                                                                   <  1:2232325/150‑1 (MQ=255)
                                                                   gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGc                                                                                   <  2:48550/150‑1 (MQ=255)
                                                                   gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGc                                                                                   <  2:2793150/150‑1 (MQ=255)
                                                                    aCCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCt                                                                                  >  2:2117900/1‑150 (MQ=255)
                                                                    aCCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGAGATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCt                                                                                  >  2:4305743/1‑150 (MQ=255)
                                                                          aTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATcc                                                                            >  2:4060164/1‑150 (MQ=255)
                                                                          aTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATcc                                                                            <  2:4600907/150‑1 (MQ=255)
                                                                          aTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATcc                                                                            <  2:3513447/150‑1 (MQ=255)
                                                                          aTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATcc                                                                            <  1:4242239/150‑1 (MQ=255)
                                                                            cGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGAGATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAg                                                                          >  2:2719950/1‑150 (MQ=255)
                                                                             gTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGt                                                                         <  1:3194916/150‑1 (MQ=255)
                                                                             gTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGt                                                                         >  1:2299749/1‑150 (MQ=255)
                                                                               cAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCACATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGttt                                                                       >  1:4195368/1‑150 (MQ=255)
                                                                               cAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGttt                                                                       >  1:4615909/1‑150 (MQ=255)
                                                                                aGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTc                                                                      >  1:2564177/1‑150 (MQ=255)
                                                                                aGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTc                                                                      >  2:3796462/1‑150 (MQ=255)
                                                                                      aGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTg                                                                <  1:4294388/150‑1 (MQ=255)
                                                                                         ccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGaca                                                             <  1:2066765/150‑1 (MQ=255)
                                                                                          cTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGacac                                                            >  2:3165674/1‑150 (MQ=255)
                                                                                              aGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtc                                                        <  1:1910576/150‑1 (MQ=255)
                                                                                                 ccGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcggg                                                     <  2:383659/150‑1 (MQ=255)
                                                                                                 ccGAACATCAATCTCTTGCCACCACAGCGGCCCTGCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcggg                                                     >  1:2291250/1‑150 (MQ=255)
                                                                                                   gAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggtc                                                   <  1:1074754/150‑1 (MQ=255)
                                                                                                   gAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggtc                                                   <  2:1150613/150‑1 (MQ=255)
                                                                                                      cATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCa                                                >  2:2318051/1‑150 (MQ=255)
                                                                                                       aTCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGGAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCAc                                               >  2:4082421/1‑150 (MQ=255)
                                                                                                          aaTCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTgc                                            >  2:4223391/1‑150 (MQ=255)
                                                                                                           aTCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTgcg                                           >  2:1453216/1‑150 (MQ=255)
                                                                                                               ctTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTTAGCGTGACACGGTCGGGTCGCACTGcgccgc                                       <  2:1095932/150‑1 (MQ=255)
                                                                                                                   ccaccaCAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTAcc                                   >  1:3129766/1‑150 (MQ=255)
                                                                                                                     accacAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTGTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGt                                 >  2:3099783/1‑150 (MQ=255)
                                                                                                                      ccacAGCGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTc                                <  2:748508/150‑1 (MQ=255)
                                                                                                                       cacaGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCg                               <  1:975956/150‑1 (MQ=255)
                                                                                                                            cGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgtt                          <  1:1953275/150‑1 (MQ=255)
                                                                                                                            cGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgtt                          >  2:3376408/1‑150 (MQ=255)
                                                                                                                            cGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgtt                          <  2:364613/150‑1 (MQ=255)
                                                                                                                            cGGCCCTTCACTGGACTCGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgtt                          >  2:3210601/1‑150 (MQ=255)
                                                                                                                               cccTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgttgta                       >  1:3677890/1‑149 (MQ=255)
                                                                                                                                 cTTCACTGGACTCGACCATCGCATCACCTTGACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTGATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCt                     <  2:28553/150‑1 (MQ=255)
                                                                                                                                 cTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCt                     <  1:4024379/150‑1 (MQ=255)
                                                                                                                                           cTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCtata           <  2:4697096/150‑1 (MQ=255)
                                                                                                                                            tCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCGGTTGTGCTGCCCGCaatat          >  1:1599667/1‑149 (MQ=255)
                                                                                                                                            tCGACCATCGCATAACCTTGACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGAGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATag          <  1:2321340/150‑1 (MQ=255)
                                                                                                                                             cGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagc         <  1:4131638/150‑1 (MQ=255)
                                                                                                                                             cGACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagc         <  2:3585515/150‑1 (MQ=255)
                                                                                                                                              gACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagcc        <  1:1167969/150‑1 (MQ=255)
                                                                                                                                              gACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagcc        >  1:397298/1‑150 (MQ=255)
                                                                                                                                              gACCATAGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCAGCCAGCTATatcc        >  1:3429923/1‑150 (MQ=255)
                                                                                                                                                  aTCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagccagcc    >  2:387953/1‑150 (MQ=255)
                                                                                                                                                    cGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCa  >  2:2624317/1‑150 (MQ=255)
                                                                                                                                                    cGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCa  <  2:1648032/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAGATTTATCGCCAGGACGTACCACCAGCGTACTTAAATACAGATCATGACGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCA  >  CP009273/2643138‑2643435

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: