Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,645,592 C→A 100% D19Y (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,645,5920CA98.0% 800.3 / 4.1 253D19Y (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 (3/2);  new base A (128/120);  total (131/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 253 total aligned reads displayed.

AGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATA  >  CP009273/2645447‑2645736
                                                                                                                                                 |                                                                                                                                                
 gCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCa                                                                                                                                               <  1:2378554/148‑1 (MQ=255)
 gCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCa                                                                                                                                               <  1:589471/148‑1 (MQ=255)
   aCCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAcc                                                                                                                                             <  1:1890496/148‑1 (MQ=255)
    ccTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAccc                                                                                                                                            >  2:2735489/1‑148 (MQ=255)
       tACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgc                                                                                                                                         <  2:342014/148‑1 (MQ=255)
          cGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccag                                                                                                                                      >  2:1317305/1‑148 (MQ=255)
          cGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccag                                                                                                                                      >  1:2011448/1‑148 (MQ=255)
             cACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagcgc                                                                                                                                   <  2:223278/148‑1 (MQ=255)
               ctctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagcgcca                                                                                                                                 >  1:2753905/1‑148 (MQ=255)
               ctctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagcgcca                                                                                                                                 >  1:971470/1‑148 (MQ=255)
                tctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagcgccag                                                                                                                                >  2:1646361/1‑148 (MQ=255)
                    ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCAtt                                                                                                                            <  2:2969168/148‑1 (MQ=255)
                    ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCAtt                                                                                                                            >  1:3079267/1‑148 (MQ=255)
                    ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCAtt                                                                                                                            <  1:3360515/148‑1 (MQ=255)
                    ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCAtt                                                                                                                            >  2:2437429/1‑148 (MQ=255)
                     tGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTa                                                                                                                           >  2:466750/1‑148 (MQ=255)
                     tGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTa                                                                                                                           >  2:463787/1‑148 (MQ=255)
                      gggCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAg                                                                                                                          <  2:304644/148‑1 (MQ=255)
                         ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCAt                                                                                                                       >  1:1658264/1‑148 (MQ=255)
                         ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCAt                                                                                                                       >  1:2938996/1‑148 (MQ=255)
                          cAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATg                                                                                                                      >  2:2604930/1‑148 (MQ=255)
                             tttttGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAg                                                                                                                   <  2:1754594/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                               <  1:483020/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                               <  2:1180747/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                               <  1:2274491/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                               >  2:809229/1‑148 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                               >  1:1035297/1‑148 (MQ=255)
                                       tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCACATCTTTTTTCACCGCTGTTGGCGCGTGATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcgc                                                                                                         >  1:3295833/1‑148 (MQ=255)
                                        ccGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcgcg                                                                                                        <  2:2001700/148‑1 (MQ=255)
                                         cGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGt                                                                                                       <  2:1053994/148‑1 (MQ=255)
                                         cGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGt                                                                                                       <  1:2320370/148‑1 (MQ=255)
                                           gCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTaa                                                                                                     >  1:855085/1‑148 (MQ=255)
                                            cACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAAc                                                                                                    <  2:1876914/148‑1 (MQ=255)
                                            cACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAAc                                                                                                    >  1:1755040/1‑148 (MQ=255)
                                              cGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtt                                                                                                  >  2:2841928/1‑148 (MQ=255)
                                              cGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtt                                                                                                  <  1:2430534/148‑1 (MQ=255)
                                              cGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtt                                                                                                  >  1:656239/1‑148 (MQ=255)
                                               ggAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACttt                                                                                                 <  2:2490288/148‑1 (MQ=255)
                                               ggAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACttt                                                                                                 >  2:122206/1‑148 (MQ=255)
                                               ggAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACttt                                                                                                 <  1:2925472/148‑1 (MQ=255)
                                                gagaGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTGTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTAGCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCGACGCGTAACtttt                                                                                                >  2:55589/1‑148 (MQ=255)
                                                 agagTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACttttt                                                                                               <  1:2815402/148‑1 (MQ=255)
                                                     ttttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatt                                                                                           >  2:1599003/1‑148 (MQ=255)
                                                       ttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatttt                                                                                         <  1:526340/148‑1 (MQ=255)
                                                         gCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatttttc                                                                                       <  2:2137250/148‑1 (MQ=255)
                                                          cTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatttttca                                                                                      <  2:3502072/148‑1 (MQ=255)
                                                             aaCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATcc                                                                                   <  1:3204918/148‑1 (MQ=255)
                                                              aCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTCCGCGGAACTTTTTCATTTTTCATccc                                                                                  >  1:1132723/1‑148 (MQ=255)
                                                                  gcgTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGcc                                                                              <  1:157693/148‑1 (MQ=255)
                                                                   cgTTTTCTGAAGAGGCCGCTTAAGTAACTTCAGAAGGCGCATCTATTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATTCAGGCGGCTCCGCGTAACTTTATCATTTTTAATCCCTGAca                                                                             >  2:1701060/1‑148 (MQ=255)
                                                                      tttCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACaa                                                                          <  1:1974399/148‑1 (MQ=255)
                                                                          tGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGc                                                                      <  1:2593851/148‑1 (MQ=255)
                                                                          tGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGc                                                                      <  1:805379/148‑1 (MQ=255)
                                                                            aaGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGccc                                                                    >  1:367343/1‑148 (MQ=255)
                                                                              gagGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGt                                                                  <  2:1084350/148‑1 (MQ=255)
                                                                                 gCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAg                                                               <  2:251188/148‑1 (MQ=255)
                                                                                   cGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCa                                                             <  1:3346067/148‑1 (MQ=255)
                                                                                   cGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCa                                                             >  2:1750297/1‑148 (MQ=255)
                                                                                    gCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCaa                                                            <  1:1372555/148‑1 (MQ=255)
                                                                                      tttAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACg                                                          >  2:1796158/1‑148 (MQ=255)
                                                                                        tAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCa                                                             <  2:509280/143‑1 (MQ=255)
                                                                                          gTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGa                                                      >  1:2989145/1‑148 (MQ=255)
                                                                                           tAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaa                                                     <  1:1468868/148‑1 (MQ=255)
                                                                                           tAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaa                                                     <  2:1657479/148‑1 (MQ=255)
                                                                                            aaCTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaat                                                    >  2:1303020/1‑148 (MQ=255)
                                                                                              cTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaataa                                                  >  1:2524638/1‑148 (MQ=255)
                                                                                              cTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaataa                                                  <  2:33949/148‑1 (MQ=255)
                                                                                              cTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaataa                                                  >  1:1311365/1‑148 (MQ=255)
                                                                                                tCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCATGCCACAATGGCCCGTTAGCAACGTCGAATAatt                                                >  2:642675/1‑148 (MQ=255)
                                                                                                 cAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAatta                                               >  2:1363005/1‑148 (MQ=255)
                                                                                                 cAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAatta                                               <  2:3529785/148‑1 (MQ=255)
                                                                                                   gAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCTGCTACTCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAattatt                                             <  1:1931776/148‑1 (MQ=255)
                                                                                                   gAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAattatt                                             >  1:1235632/1‑148 (MQ=255)
                                                                                                   gAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAattatt                                             <  1:130582/148‑1 (MQ=255)
                                                                                                      ggCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATg                                          <  2:874970/148‑1 (MQ=255)
                                                                                                        cgcATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCg                                        >  1:3399219/1‑148 (MQ=255)
                                                                                                           aTCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTga                                     <  1:1193139/148‑1 (MQ=255)
                                                                                                           aTCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTga                                     <  1:508901/148‑1 (MQ=255)
                                                                                                             cTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTgaga                                   >  1:106031/1‑148 (MQ=255)
                                                                                                             cTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTgaga                                   <  1:1843359/148‑1 (MQ=255)
                                                                                                             cTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTgaga                                   >  2:3124457/1‑148 (MQ=255)
                                                                                                              ttttttCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGaa                                  <  2:1760938/148‑1 (MQ=255)
                                                                                                                 tttCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAAttt                               >  2:3187533/1‑148 (MQ=255)
                                                                                                                     aCCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGt                           <  1:1646361/148‑1 (MQ=255)
                                                                                                                         cTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATg                       >  1:3388113/1‑148 (MQ=255)
                                                                                                                         cTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATg                       >  2:3081560/1‑148 (MQ=255)
                                                                                                                         cTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCCCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATACCTGACACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATg                       >  1:1047496/1‑148 (MQ=255)
                                                                                                                           gTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTa                     >  2:1521782/1‑148 (MQ=255)
                                                                                                                            ttGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTaa                    <  2:3079267/148‑1 (MQ=255)
                                                                                                                             tGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAAt                   >  2:3377897/1‑148 (MQ=255)
                                                                                                                             tGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAAt                   >  2:1611150/1‑148 (MQ=255)
                                                                                                                                cgcgTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCa                <  1:2437429/148‑1 (MQ=255)
                                                                                                                                  cgTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACATTATGTAATTCAGt              <  2:3262129/148‑1 (MQ=255)
                                                                                                                                   gTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACACTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTa             <  2:697778/148‑1 (MQ=255)
                                                                                                                                     tATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATc           <  2:2082120/148‑1 (MQ=255)
                                                                                                                                     tATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATc           <  2:855085/148‑1 (MQ=255)
                                                                                                                                           tgttgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGccc     <  1:2091444/148‑1 (MQ=255)
                                                                                                                                           tgttgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGccc     <  1:1763229/148‑1 (MQ=255)
                                                                                                                                             ttgtCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACACTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCAt   >  1:2884102/1‑148 (MQ=255)
                                                                                                                                              tgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATa  >  1:2203843/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATA  >  CP009273/2645447‑2645736

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: ‑