Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,645,592 C→A 85.6% 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,5920CA85.6% 355.1 / 65.6 160D19Y (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 (14/9);  new base A (59/78);  total (73/87)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 160 total aligned reads displayed.

AGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACA  >  CP009273/2645447‑2645738
                                                                                                                                                 |                                                                                                                                                  
aGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGc                                                                                                                                                  >  1:540111/1‑148 (MQ=255)
 gCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATAGTTGTTGTAGCa                                                                                                                                                 >  1:1747689/1‑148 (MQ=255)
   aCCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCAcc                                                                                                                                               <  1:2166499/148‑1 (MQ=255)
   aCCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCAcc                                                                                                                                               <  1:2165245/148‑1 (MQ=255)
   aCCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAcc                                                                                                                                               <  1:561566/148‑1 (MQ=255)
   aCCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAcc                                                                                                                                               <  2:1116290/148‑1 (MQ=255)
   aCCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAcc                                                                                                                                               <  2:2685750/148‑1 (MQ=255)
    ccTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAccc                                                                                                                                              <  2:1459827/148‑1 (MQ=255)
    ccTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAccc                                                                                                                                              <  1:2144654/148‑1 (MQ=255)
       tACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgc                                                                                                                                           >  2:2532625/1‑148 (MQ=255)
           ggCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagc                                                                                                                                       >  1:881461/1‑148 (MQ=255)
                 ctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgccagc                                                                                                                                 >  1:3020527/1‑148 (MQ=255)
                       ggCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGc                                                                                                                           >  2:2499528/1‑148 (MQ=255)
                       ggCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGc                                                                                                                           >  1:2988263/1‑148 (MQ=255)
                         ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCAt                                                                                                                         >  1:1044636/1‑148 (MQ=255)
                          cAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATg                                                                                                                        >  2:1913655/1‑148 (MQ=255)
                             tttttGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAg                                                                                                                     <  2:2858997/148‑1 (MQ=255)
                                ttGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGACGGCGCATCTTTTTTCACCGGTGTTGGCGCGTTATCGTTGGTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcg                                                                                                                  >  1:1942891/1‑148 (MQ=255)
                                ttGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcg                                                                                                                  >  1:273482/1‑148 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                                 <  1:2746867/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                                 <  2:314890/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                                 <  2:2539185/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                                 <  2:2539436/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                                 <  1:1371213/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                                 <  2:1255970/148‑1 (MQ=255)
                                 tGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                                 <  1:3047238/148‑1 (MQ=255)
                                  gTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcggc                                                                                                                <  1:769921/148‑1 (MQ=255)
                                  gTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACGTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcggc                                                                                                                <  2:73424/148‑1 (MQ=255)
                                   tctctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCt                                                                                                               >  2:1442753/1‑148 (MQ=255)
                                   tctctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCt                                                                                                               >  2:2819650/1‑148 (MQ=255)
                                   tctctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGGTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCt                                                                                                               >  2:1380115/1‑148 (MQ=255)
                                   tctctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTGTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGGTGTTGTAGCACCCTGCCAGCGCCCGCATTAGCATGCAGGCGGCt                                                                                                               >  2:3084645/1‑148 (MQ=255)
                                        ccGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcgcg                                                                                                          >  1:1406335/1‑148 (MQ=255)
                                         cGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGt                                                                                                         <  1:3141260/148‑1 (MQ=255)
                                          ggCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTa                                                                                                        <  2:2631082/148‑1 (MQ=255)
                                            cACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAAc                                                                                                      >  1:2944586/1‑148 (MQ=255)
                                            cACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAAc                                                                                                      >  1:1293793/1‑148 (MQ=255)
                                               ggAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACttt                                                                                                   <  2:2825016/148‑1 (MQ=255)
                                                gagaGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg                                                                                                                 <  1:104003/133‑1 (MQ=255)
                                                gagaGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttt                                                                                                  >  1:2879615/1‑148 (MQ=255)
                                                gagaGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttt                                                                                                  >  1:1286671/1‑148 (MQ=255)
                                                 agagTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACttttt                                                                                                 >  2:87283/1‑148 (MQ=255)
                                                          cTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatttttca                                                                                        <  1:883997/148‑1 (MQ=255)
                                                            gAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATc                                                                                      <  2:1825397/148‑1 (MQ=255)
                                                             aaCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATcc                                                                                     <  1:401027/148‑1 (MQ=255)
                                                             aaCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATcc                                                                                     <  1:1162967/148‑1 (MQ=255)
                                                             aaCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATcc                                                                                     <  1:505192/148‑1 (MQ=255)
                                                             aaCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATcc                                                                                     <  1:505194/148‑1 (MQ=255)
                                                                tCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTg                                                                                  <  2:881461/148‑1 (MQ=255)
                                                                       ttCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAAt                                                                           <  2:2585015/148‑1 (MQ=255)
                                                                        tCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATg                                                                          <  2:2703762/148‑1 (MQ=255)
                                                                            aaGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGccc                                                                      >  1:2864618/1‑148 (MQ=255)
                                                                             agagGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCg                                                                     >  1:790420/1‑148 (MQ=255)
                                                                              gagGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGt                                                                    >  1:1062145/1‑148 (MQ=255)
                                                                                   cGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCa                                                               <  1:1778419/148‑1 (MQ=255)
                                                                                     cTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAAc                                                             >  1:1918586/1‑148 (MQ=255)
                                                                                     cTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAAc                                                             >  2:2481430/1‑148 (MQ=255)
                                                                                     cTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAAc                                                             >  2:2142216/1‑148 (MQ=255)
                                                                                       ttAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTTTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGt                                                           <  1:3041958/148‑1 (MQ=255)
                                                                                       ttAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGt                                                           >  2:767286/1‑148 (MQ=255)
                                                                                             aCTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaata                                                     >  1:349725/1‑148 (MQ=255)
                                                                                             aCTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaata                                                     >  1:1412832/1‑148 (MQ=255)
                                                                                              cTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaataa                                                    <  2:334978/148‑1 (MQ=255)
                                                                                              cTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaataa                                                    >  2:2684932/1‑148 (MQ=255)
                                                                                                tCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCACCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAatt                                                  <  1:545417/148‑1 (MQ=255)
                                                                                                   gAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAattatt                                               <  1:773384/148‑1 (MQ=255)
                                                                                                   gAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGATAGCAACGTCGAAGAattatt                                               <  2:3190628/148‑1 (MQ=255)
                                                                                                      ggCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATg                                            >  2:2938632/1‑148 (MQ=255)
                                                                                                      ggCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATg                                            >  2:2938627/1‑148 (MQ=255)
                                                                                                       gcgcATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGc                                           >  2:1978361/1‑148 (MQ=255)
                                                                                                          cATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTCCGCGTCACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGACTAATTATTATGCGTg                                        >  1:2629979/1‑148 (MQ=255)
                                                                                                          cATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTg                                        >  1:1499605/1‑148 (MQ=255)
                                                                                                           aTCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTga                                       >  1:2546073/1‑148 (MQ=255)
                                                                                                             cTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTgaga                                     >  1:1064144/1‑148 (MQ=255)
                                                                                                             cTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTgaga                                     >  1:42778/1‑148 (MQ=255)
                                                                                                              ttttttCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGaa                                    <  2:539875/148‑1 (MQ=255)
                                                                                                              ttttttCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGaa                                    <  2:540111/148‑1 (MQ=255)
                                                                                                                ttttCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAAtt                                  <  1:2520213/148‑1 (MQ=255)
                                                                                                                 tttCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAAttt                                 <  1:634984/148‑1 (MQ=255)
                                                                                                                 tttCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAAttt                                 <  1:2297669/148‑1 (MQ=255)
                                                                                                                    cACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACg                              <  1:3059526/148‑1 (MQ=255)
                                                                                                                      ccGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGtt                            <  1:700733/148‑1 (MQ=255)
                                                                                                                       cGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTa                           >  1:1907744/1‑148 (MQ=255)
                                                                                                                          tGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGt                        >  1:2513525/1‑148 (MQ=255)
                                                                                                                            ttGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTaa                      <  2:3307207/148‑1 (MQ=255)
                                                                                                                                cgcgTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCa                  >  2:853834/1‑148 (MQ=255)
                                                                                                                                 gcgTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAg                 <  1:2004010/148‑1 (MQ=255)
                                                                                                                                 gcgTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTATTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAg                 <  1:2003853/148‑1 (MQ=255)
                                                                                                                                    ttATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTAt              >  1:223175/1‑148 (MQ=255)
                                                                                                                                     tATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATc             >  2:1533248/1‑148 (MQ=255)
                                                                                                                                     tATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAAATATTATGCGTGAGAATTTACGTTATGTAATACAGTATc             >  2:1002690/1‑148 (MQ=255)
                                                                                                                                      aTCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCa            <  1:2532625/148‑1 (MQ=255)
                                                                                                                                        cGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACg          >  2:2034448/1‑148 (MQ=255)
                                                                                                                                           tgttgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGccc       <  1:311665/148‑1 (MQ=255)
                                                                                                                                            gttgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCa      <  2:619340/148‑1 (MQ=255)
                                                                                                                                             ttgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACTCGTAACTTTTTCCGTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTAGGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCAt     <  1:1380115/148‑1 (MQ=255)
                                                                                                                                             ttgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCAt     >  2:2454328/1‑148 (MQ=255)
                                                                                                                                             ttgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCAt     >  2:2454617/1‑148 (MQ=255)
                                                                                                                                             ttgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCAt     <  2:1027547/148‑1 (MQ=255)
                                                                                                                                              tgtAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATa    <  1:2500722/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACA  >  CP009273/2645447‑2645738

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: