Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,643,336 C→A 86.8% D771Y (GAC→TAC)  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,3360CA86.8% 678.0 / 110.8 282D771Y (GAC→TAC) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base A (138/107);  total (157/125)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 284 total aligned reads displayed.

CGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGA  >  CP009273/2643188‑2643484
                                                                                                                                                    |                                                                                                                                                    
  aTTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGc                                                                                                                                                    <  2:145843/149‑1 (MQ=255)
  aTTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGc                                                                                                                                                    >  1:3284010/1‑149 (MQ=255)
   ttCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGcc                                                                                                                                                   >  2:3909776/1‑149 (MQ=255)
   ttCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGcc                                                                                                                                                   >  2:1520570/1‑149 (MQ=255)
   ttCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGcc                                                                                                                                                   >  2:2854725/1‑149 (MQ=255)
   ttCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGcc                                                                                                                                                   >  2:4406860/1‑149 (MQ=255)
     ccAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAg                                                                                                                                                 >  1:2069365/1‑149 (MQ=255)
      cAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAgg                                                                                                                                                <  2:314197/149‑1 (MQ=255)
      cAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAgg                                                                                                                                                <  2:4559878/149‑1 (MQ=255)
        ggTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGc                                                                                                                                              >  1:693176/1‑149 (MQ=255)
        ggTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGc                                                                                                                                              >  2:1814164/1‑149 (MQ=255)
              aTCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAAc                                                                                                                                        <  1:1070777/149‑1 (MQ=255)
               tCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACt                                                                                                                                       >  2:2465434/1‑149 (MQ=255)
                 gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTgg                                                                                                                                     >  2:4207256/1‑149 (MQ=255)
                 gACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTgg                                                                                                                                     >  2:1936466/1‑149 (MQ=255)
                      gAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTa                                                                                                                                <  2:4470911/149‑1 (MQ=255)
                       aaTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTAt                                                                                                                               <  1:41836/149‑1 (MQ=255)
                       aaTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTAt                                                                                                                               <  2:400299/149‑1 (MQ=255)
                        atagTCAGATACAGCCATTGAGCACGAACATCAATCTCTTGCAACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACATTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATc                                                                                                                              >  2:69512/4‑149 (MQ=255)
                        aTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATc                                                                                                                              >  1:2136868/1‑149 (MQ=255)
                         tCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATcc                                                                                                                             <  1:3670026/149‑1 (MQ=255)
                            tCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGt                                                                                                                          >  2:1008105/1‑149 (MQ=255)
                             cAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGtt                                                                                                                         <  1:2975277/149‑1 (MQ=255)
                               gATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTc                                                                                                                       >  1:4007203/1‑149 (MQ=255)
                               gATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTc                                                                                                                       >  1:4317365/1‑149 (MQ=255)
                               gATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTc                                                                                                                       >  1:1398902/1‑149 (MQ=255)
                                 tCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAg                                                                                                                     <  1:4238061/149‑1 (MQ=255)
                                 tCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAg                                                                                                                     <  2:1069104/149‑1 (MQ=255)
                                      cccTTGAGCCCGACCATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGa                                                                                                                <  2:4258135/149‑1 (MQ=255)
                                      cccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGa                                                                                                                <  2:2195213/149‑1 (MQ=255)
                                       ccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGac                                                                                                               <  2:3458536/149‑1 (MQ=255)
                                             gCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtc                                                                                                         <  2:2663545/149‑1 (MQ=255)
                                                 gAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggt                                                                                                     <  1:750448/149‑1 (MQ=255)
                                                  aaCATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggtc                                                                                                    >  1:1921800/1‑149 (MQ=255)
                                                  aaCATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggtc                                                                                                    >  2:1510321/1‑149 (MQ=255)
                                                      tCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCAc                                                                                                <  1:4777774/149‑1 (MQ=255)
                                                       cAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACt                                                                                               <  1:1240618/149‑1 (MQ=255)
                                                        aaTCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTg                                                                                              >  2:3081595/1‑149 (MQ=255)
                                                         aTCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCCGTTTCAGCGTGACACGGTCGGGTCGCACTgc                                                                                             >  2:4521686/1‑149 (MQ=255)
                                                          tctcTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTgcg                                                                                            >  2:671600/1‑149 (MQ=255)
                                                           ctctTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTgcgc                                                                                           <  2:2469378/149‑1 (MQ=255)
                                                               tGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCt                                                                                       >  2:3613266/1‑149 (MQ=255)
                                                               tGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCt                                                                                       >  2:1240442/1‑149 (MQ=255)
                                                               tGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCt                                                                                       <  1:3073025/149‑1 (MQ=255)
                                                                gCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTa                                                                                      <  2:3933274/149‑1 (MQ=255)
                                                                 ccaccaCAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTAc                                                                                     <  2:3693360/149‑1 (MQ=255)
                                                                 ccaccaCAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTAc                                                                                     <  1:2711139/149‑1 (MQ=255)
                                                                   accacAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCg                                                                                   >  2:534298/1‑149 (MQ=255)
                                                                    ccacAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGt                                                                                  <  1:1100897/149‑1 (MQ=255)
                                                                    ccacAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGt                                                                                  >  1:4557290/1‑149 (MQ=255)
                                                                      acaGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCg                                                                                >  1:452725/1‑149 (MQ=255)
                                                                      acaGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCg                                                                                >  1:2378654/1‑149 (MQ=255)
                                                                         gCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtg                                                                             >  2:3844043/1‑149 (MQ=255)
                                                                          cGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgt                                                                            <  2:2572879/149‑1 (MQ=255)
                                                                             cccTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgttgt                                                                         >  2:3373227/1‑149 (MQ=255)
                                                                               cTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTcc                                                                       >  2:2269869/1‑149 (MQ=255)
                                                                               cTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTcc                                                                       >  2:3128150/1‑149 (MQ=255)
                                                                                ttCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCt                                                                      <  1:4399027/149‑1 (MQ=255)
                                                                                     tGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAg                                                                 >  2:1926588/1‑149 (MQ=255)
                                                                                      ggACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGc                                                                >  1:3557598/1‑149 (MQ=255)
                                                                                       gACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCt                                                               <  1:1682270/149‑1 (MQ=255)
                                                                                       gACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCt                                                               >  1:4116780/1‑149 (MQ=255)
                                                                                        aCTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCaccaa                                                              >  1:4415201/1‑145 (MQ=255)
                                                                                         cTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCtat                                                             >  2:2358989/1‑149 (MQ=255)
                                                                                             aCCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagcc                                                         <  2:141569/149‑1 (MQ=255)
                                                                                             aCCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCCGCTATagcc                                                         >  1:970013/1‑149 (MQ=255)
                                                                                             aCCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagcc                                                         <  2:862432/149‑1 (MQ=255)
                                                                                             aCCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagcc                                                         <  1:4074631/149‑1 (MQ=255)
                                                                                             aCCATCGCATAACCATGACCCGCCGTTGGCGCGGCGAGATGCAACTTAATGGTGTCGCCAGGGCGATAACGGGCTTTATCCAGTTTAAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagcc                                                         <  2:3386224/149‑1 (MQ=255)
                                                                                              ccATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagcca                                                        >  1:4698309/1‑149 (MQ=255)
                                                                                                       aaCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGaaa                                               >  1:2756758/1‑149 (MQ=255)
                                                                                                        aCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTGCTGCCAGCTATAGCCAGCCCAGAAAc                                              >  1:4583629/1‑149 (MQ=255)
                                                                                                        aCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAAc                                              <  1:3915467/149‑1 (MQ=255)
                                                                                                         ccTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACg                                             <  1:1850456/149‑1 (MQ=255)
                                                                                                         ccTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACg                                             >  2:2752067/1‑149 (MQ=255)
                                                                                                         ccTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACg                                             >  2:4046509/1‑149 (MQ=255)
                                                                                                          cTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGa                                            >  2:4746913/1‑149 (MQ=255)
                                                                                                             tACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCGGTTGTCCTGCCCGCTATAGCCACCCCAGAAACCaaaa                                         >  1:1635695/1‑147 (MQ=255)
                                                                                                              aCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAacac                                        <  1:4460688/149‑1 (MQ=255)
                                                                                                              aCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCAGCCCCCCAAAGCCACCCCACAAAACAacac                                        >  1:828139/1‑149 (MQ=255)
                                                                                                                 cgccgTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACactac                                     <  2:4516390/149‑1 (MQ=255)
                                                                                                                 cgccgTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACactac                                     <  1:3958221/149‑1 (MQ=255)
                                                                                                                 cgccgTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACactac                                     <  1:1568775/149‑1 (MQ=255)
                                                                                                                 cgccgTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACactac                                     <  1:770560/149‑1 (MQ=255)
                                                                                                                  gccgTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGACACGAACactact                                    >  2:4375471/1‑149 (MQ=255)
                                                                                                                   ccgTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTg                                   >  2:1763589/1‑149 (MQ=255)
                                                                                                                     gTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGAc                                 <  2:3284010/149‑1 (MQ=255)
                                                                                                                     gTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGAc                                 <  2:2112257/149‑1 (MQ=255)
                                                                                                                      ttGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCAGCCCCCTATAGCCACCCCAGAAACGAACACAACAGAcc                                >  1:20025/1‑149 (MQ=255)
                                                                                                                         gcgcGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCt                             <  2:3189288/149‑1 (MQ=255)
                                                                                                                              gcgATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCAtt                        <  2:2062015/149‑1 (MQ=255)
                                                                                                                               cgATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTc                       >  2:959106/1‑149 (MQ=255)
                                                                                                                                gATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCCGCGATAGCCAGCCCAGAAACGAACACTACAGACCGCGTCAATCg                      >  1:3714668/1‑149 (MQ=255)
                                                                                                                                   atGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTAGgcg                   >  2:2826425/1‑149 (MQ=255)
                                                                                                                                    tGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGgcgc                  <  1:526440/149‑1 (MQ=255)
                                                                                                                                             aTGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTc         >  2:854212/1‑149 (MQ=255)
                                                                                                                                             aTGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCAGCCCGCGATATCCCGCCCAGAAACGAACACAACTGACCGCTTCATTAGGCGCTTTGACTTc         >  1:1431128/1‑149 (MQ=255)
                                                                                                                                              tGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTcc        <  1:2597983/149‑1 (MQ=255)
                                                                                                                                               ggTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCGGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACAGAACGCTTAATTAGGCGCTTTGACTTCCa       >  1:2221557/1‑149 (MQ=255)
                                                                                                                                                gtgtCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAg      <  1:3452912/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGA  >  CP009273/2643188‑2643484

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: