Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,643,336 C→A 100% 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,3360CA98.0% 788.4 / 5.3 245D771Y (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 (4/1);  new base A (113/127);  total (117/128)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 245 total aligned reads displayed.

ATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGA  >  CP009273/2643190‑2643481
                                                                                                                                                  |                                                                                                                                                 
aTTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGc                                                                                                                                                 >  1:704877/1‑149 (MQ=255)
aTTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGc                                                                                                                                                 <  1:3745462/149‑1 (MQ=255)
 ttCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGcc                                                                                                                                                <  1:935287/149‑1 (MQ=255)
  tCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCa                                                                                                                                               >  1:3993746/1‑149 (MQ=255)
      ggTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACAACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTAACCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGc                                                                                                                                           >  1:385596/1‑149 (MQ=255)
        ttttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGa                                                                                                                                         >  2:982103/1‑149 (MQ=255)
         tttATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGAt                                                                                                                                        >  1:3504996/1‑149 (MQ=255)
                 ccGGAATCGTCAGATCCCGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCt                                                                                                                                <  2:1440237/149‑1 (MQ=255)
                 ccGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCt                                                                                                                                <  2:1440713/149‑1 (MQ=255)
                   ggAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCttt                                                                                                                              >  2:3298891/1‑149 (MQ=255)
                     aaTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTAt                                                                                                                            >  1:1566903/1‑149 (MQ=255)
                     aaTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCACTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTAt                                                                                                                            >  1:1941007/1‑149 (MQ=255)
                      aTCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATc                                                                                                                           >  1:1996993/1‑149 (MQ=255)
                       tCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATcc                                                                                                                          >  2:3677083/1‑149 (MQ=255)
                         gTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAg                                                                                                                        >  1:4107875/1‑149 (MQ=255)
                          tCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGt                                                                                                                       >  1:3661817/1‑149 (MQ=255)
                           cAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGtt                                                                                                                      <  2:25139/149‑1 (MQ=255)
                              aTCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCa                                                                                                                   <  2:2425207/149‑1 (MQ=255)
                               tCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATACAGTTTCAg                                                                                                                  >  1:2945888/1‑149 (MQ=255)
                                ccAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGc                                                                                                                 <  2:35040/149‑1 (MQ=255)
                                  aGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGt                                                                                                               <  1:3943856/149‑1 (MQ=255)
                                    cccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGa                                                                                                             >  2:463977/1‑149 (MQ=255)
                                    cccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGa                                                                                                             <  2:846376/149‑1 (MQ=255)
                                     ccTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGac                                                                                                            <  1:2286658/149‑1 (MQ=255)
                                      cTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGaca                                                                                                           >  1:579217/1‑149 (MQ=255)
                                      cTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGaca                                                                                                           >  1:2645701/1‑149 (MQ=255)
                                          aGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggt                                                                                                       <  1:1226204/149‑1 (MQ=255)
                                             ccGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTGTCTGCGTGACACggtcgg                                                                                                    >  2:2464142/1‑149 (MQ=255)
                                               gAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggt                                                                                                  <  2:1561279/149‑1 (MQ=255)
                                                aaCATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACggtcgggtc                                                                                                 >  2:693736/1‑149 (MQ=255)
                                                  cATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGc                                                                                               <  1:1451295/149‑1 (MQ=255)
                                                   aTCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCa                                                                                              >  1:1896128/1‑149 (MQ=255)
                                                   aTCAATCTCTAGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCTACTTACTGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCa                                                                                              <  1:3216388/149‑1 (MQ=255)
                                                    tCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCAc                                                                                             <  2:1662053/149‑1 (MQ=255)
                                                    tCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCAc                                                                                             >  2:2694219/1‑149 (MQ=255)
                                                    tCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCAc                                                                                             >  1:2512751/1‑149 (MQ=255)
                                                      aaTCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTg                                                                                           >  1:1069679/1‑149 (MQ=255)
                                                         ctctTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTgcgc                                                                                        <  1:452704/149‑1 (MQ=255)
                                                          tctTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGcgcc                                                                                       <  1:1135245/149‑1 (MQ=255)
                                                            ttGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGcgccgc                                                                                     <  2:2334535/149‑1 (MQ=255)
                                                             tGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAGCTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCt                                                                                    <  1:1836894/149‑1 (MQ=255)
                                                               ccaccaCAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTAc                                                                                  <  2:2237727/149‑1 (MQ=255)
                                                                    acaGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCg                                                                             >  2:2842081/1‑149 (MQ=255)
                                                                     caGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGc                                                                            <  1:1758387/149‑1 (MQ=255)
                                                                     caGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGc                                                                            <  1:3527054/149‑1 (MQ=255)
                                                                       gCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtg                                                                          >  2:3747792/1‑149 (MQ=255)
                                                                        cGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgt                                                                         <  1:3444299/149‑1 (MQ=255)
                                                                          gCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCtgttg                                                                       >  2:3538591/1‑149 (MQ=255)
                                                                             cTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTcc                                                                    >  2:3658486/1‑149 (MQ=255)
                                                                              ttCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCTCACTGCGCCGCTACCGTCGCTGTTGTCCt                                                                   <  2:2136823/149‑1 (MQ=255)
                                                                                cACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGc                                                                 >  2:2953727/1‑149 (MQ=255)
                                                                                cACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGc                                                                 >  1:4075598/1‑149 (MQ=255)
                                                                                cACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCGGTTGTCCTGc                                                                 >  1:3737553/1‑149 (MQ=255)
                                                                                   tGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGcgccgc                                                                                     <  2:1287615/126‑1 (MQ=255)
                                                                                    ggACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGc                                                             <  2:2009505/149‑1 (MQ=255)
                                                                                       cTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCtat                                                          >  2:1991230/1‑149 (MQ=255)
                                                                                        tCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCCGCtata                                                         >  1:2193003/1‑149 (MQ=255)
                                                                                        tCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCtata                                                         >  2:2612706/1‑149 (MQ=255)
                                                                                        tCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCtata                                                         >  2:518817/1‑149 (MQ=255)
                                                                                          gACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagc                                                       >  1:2672045/1‑149 (MQ=255)
                                                                                          gACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagc                                                       >  1:2672015/1‑149 (MQ=255)
                                                                                           aCCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagcc                                                      >  2:3713626/1‑149 (MQ=255)
                                                                                             cATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCCGCTATagccag                                                    >  2:1406387/1‑149 (MQ=255)
                                                                                             cATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCTCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATagccag                                                    <  2:1070697/149‑1 (MQ=255)
                                                                                              aTCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGGCGCTGGTGTCCTGCCAGCTATagccagc                                                   >  2:2279146/1‑149 (MQ=255)
                                                                                                   aTAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGa                                              >  2:3569137/1‑149 (MQ=255)
                                                                                                   aTAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGa                                              >  2:320885/1‑149 (MQ=255)
                                                                                                    tAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGaa                                             >  2:2740042/1‑149 (MQ=255)
                                                                                                     aaCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGaaa                                            >  2:4054730/1‑149 (MQ=255)
                                                                                                      aCCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAAc                                           >  2:1693952/1‑149 (MQ=255)
                                                                                                         tttACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGaa                                        <  1:644358/149‑1 (MQ=255)
                                                                                                         tttACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGaa                                        <  2:549646/149‑1 (MQ=255)
                                                                                                          ttACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAac                                       <  2:2645701/149‑1 (MQ=255)
                                                                                                           tACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAaca                                      <  2:3477200/149‑1 (MQ=255)
                                                                                                            aCCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAacac                                     >  2:3594454/1‑149 (MQ=255)
                                                                                                              ccgccgTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACacta                                   <  2:443326/149‑1 (MQ=255)
                                                                                                               cgccgTTGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACactac                                  <  1:1878629/149‑1 (MQ=255)
                                                                                                               cgccgTTGGCGCGGCGATATGAAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACactac                                  <  1:1299020/149‑1 (MQ=255)
                                                                                                                    ttGGCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGAcc                             >  2:3065350/1‑149 (MQ=255)
                                                                                                                      ggCGCGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTAGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGc                           <  1:2887763/149‑1 (MQ=255)
                                                                                                                       gcgcGGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCt                          >  2:1293747/1‑149 (MQ=255)
                                                                                                                        cgcgGCGATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCtt                         <  1:839265/149‑1 (MQ=255)
                                                                                                                           ggcgATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCAt                      <  2:2454111/149‑1 (MQ=255)
                                                                                                                            gcgATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCAtt                     <  1:522360/149‑1 (MQ=255)
                                                                                                                              gATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCg                   >  2:1819308/1‑149 (MQ=255)
                                                                                                                              gATATGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCg                   <  1:2545013/149‑1 (MQ=255)
                                                                                                                               atatGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCgg                  >  2:3510753/1‑149 (MQ=255)
                                                                                                                               atatGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTACTGCCAGCAATAGCCAGCCCACAAACGAAAACTACTGAACCCTTAATTAgg                  >  2:3650747/1‑149 (MQ=255)
                                                                                                                                  tGCAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCGTCATTCGgcgc               >  1:1999161/1‑149 (MQ=255)
                                                                                                                                    cAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCtt             <  2:1083632/149‑1 (MQ=255)
                                                                                                                                    cAACTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCtt             <  2:91947/149‑1 (MQ=255)
                                                                                                                                      aCTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTg           >  2:3444803/1‑149 (MQ=255)
                                                                                                                                      aCTTAATGGTGTAGCCAGGGCGATAACAGGCTTTATCCATTTTCAGCGTGACAAGGTAGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGACAGCGATATCCAGACCACAAACCAAAAATAATGAAACCTTCATTCTGCGCtttt           >  1:186343/1‑148 (MQ=255)
                                                                                                                                       cTTAATGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGa          <  2:2308219/149‑1 (MQ=255)
                                                                                                                                          aaTGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGGCCTGCCAGCTATAGCCAGCCCAGAAACGAACACGACAGACCGCTTCATTAGGCGCTTTGACtt       >  2:16737/1‑149 (MQ=255)
                                                                                                                                            tGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATATCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTcc     >  1:1691664/1‑149 (MQ=255)
                                                                                                                                            tGGTGTAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTcc     <  2:1451384/149‑1 (MQ=255)
                                                                                                                                               tgtAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGa  <  2:4024588/149‑1 (MQ=255)
                                                                                                                                               tgtAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGa  <  1:897641/149‑1 (MQ=255)
                                                                                                                                               tgtAGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGa  <  1:1653335/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATTCCAGGTTTTATCGACCGGAATCGTCAGATCCAGCCCTTGAGCCCGAACATCAATCTCTTGCCACCACAGCGGCCCTTCACTGGACTCGACCATCGCATAACCTTTACCCGCCGTTGGCGCGGCGATATGCAACTTAATGGTGTCGCCAGGGCGATAACTGGCTTTATCCAGTTTCAGCGTGACACGGTCGGGTCGCACTGCGCCGCTACCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGA  >  CP009273/2643190‑2643481

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: