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,5920CA99.0% 643.4 / ‑2.0 200D19Y (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 (2/0);  new base A (94/104);  total (96/104)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 200 total aligned reads displayed.

CACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCAT  >  CP009273/2645449‑2645735
                                                                                                                                               |                                                                                                                                               
cACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAc                                                                                                                                             >  2:244594/1‑148 (MQ=255)
 aCCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAcc                                                                                                                                            >  1:118292/1‑148 (MQ=255)
 aCCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCAcc                                                                                                                                            >  1:1199359/1‑148 (MQ=255)
    ttACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTg                                                                                                                                         >  1:2076184/1‑148 (MQ=255)
     tACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgc                                                                                                                                        <  1:594236/148‑1 (MQ=255)
     tACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgc                                                                                                                                        <  1:2358468/148‑1 (MQ=255)
          gCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagcg                                                                                                                                   >  2:1427579/1‑148 (MQ=255)
          gCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagcg                                                                                                                                   >  2:2160511/1‑148 (MQ=255)
           cACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagcgc                                                                                                                                  >  1:151166/1‑148 (MQ=255)
           cACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagcgc                                                                                                                                  >  2:2507169/1‑148 (MQ=255)
               ctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTgccagcgccagc                                                                                                                              >  1:1576974/1‑148 (MQ=255)
                tGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCa                                                                                                                             <  1:3326149/148‑1 (MQ=255)
                 gTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCAt                                                                                                                            >  1:1995080/1‑148 (MQ=255)
                  ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCAtt                                                                                                                           >  2:642869/1‑148 (MQ=255)
                  ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCAtt                                                                                                                           <  1:923884/148‑1 (MQ=255)
                   tGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTa                                                                                                                          <  2:2945515/148‑1 (MQ=255)
                     ggtcAGTTGTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGc                                                                                                                        >  2:3134798/4‑148 (MQ=255)
                     ggCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGc                                                                                                                        <  1:246561/148‑1 (MQ=255)
                           tttttGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAg                                                                                                                  >  2:1791485/1‑148 (MQ=255)
                                  ctctcCGGCACGGAGAGTTGTGCTGAACTCGCTTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCGGCCAGCGCCAGCATTAGCATGCAGGCGGCTa                                                                                                           <  2:918337/148‑1 (MQ=255)
                                    ctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcg                                                                                                         >  2:1810287/1‑148 (MQ=255)
                                    ctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcg                                                                                                         >  2:1270770/1‑148 (MQ=255)
                                         gCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTaa                                                                                                    <  2:184800/148‑1 (MQ=255)
                                         gCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTaa                                                                                                    <  2:788537/148‑1 (MQ=255)
                                         gCACGGAGAGTTTAGCAGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTaa                                                                                                    <  2:672541/148‑1 (MQ=255)
                                          cACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAAc                                                                                                   >  1:3438750/1‑148 (MQ=255)
                                           aCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACt                                                                                                  >  2:1839320/1‑148 (MQ=255)
                                           aCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACt                                                                                                  <  2:1097586/148‑1 (MQ=255)
                                            cGGAGAGTTTGGCTGAACTCGCGTGTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtt                                                                                                 <  2:3104367/148‑1 (MQ=255)
                                             ggAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACttt                                                                                                >  2:3179442/1‑148 (MQ=255)
                                               agagTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGACGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACttttt                                                                                              >  1:819634/1‑148 (MQ=255)
                                                gagTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttc                                                                                             <  1:2507181/148‑1 (MQ=255)
                                                 agTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttca                                                                                            <  1:499124/148‑1 (MQ=255)
                                                 agTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttca                                                                                            <  1:3172572/148‑1 (MQ=255)
                                                 agTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAAAtttttca                                                                                            >  2:2349787/1‑148 (MQ=255)
                                                   ttttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatt                                                                                          >  2:2274367/1‑148 (MQ=255)
                                                   ttttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatt                                                                                          >  1:3286043/1‑148 (MQ=255)
                                                   ttttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatt                                                                                          <  2:1664711/148‑1 (MQ=255)
                                                   ttttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatt                                                                                          <  2:2433992/148‑1 (MQ=255)
                                                   ttttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatt                                                                                          >  2:2274479/1‑148 (MQ=255)
                                                    tttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcattt                                                                                         <  2:149955/148‑1 (MQ=255)
                                                    tttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcattt                                                                                         <  1:639231/148‑1 (MQ=255)
                                                    tttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcattt                                                                                         <  2:921777/148‑1 (MQ=255)
                                                    tttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcattt                                                                                         <  1:3451383/148‑1 (MQ=255)
                                                     ttGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcatttt                                                                                        >  2:763660/1‑148 (MQ=255)
                                                      tGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACtttttcattttt                                                                                       >  1:2324901/1‑148 (MQ=255)
                                                          gAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATc                                                                                   >  2:3474078/1‑148 (MQ=255)
                                                           aaCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTTCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATcc                                                                                  <  2:2173763/148‑1 (MQ=255)
                                                           aaCTCGCGTTGTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATcc                                                                                  <  2:3212506/148‑1 (MQ=255)
                                                            aCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATccc                                                                                 <  2:2307051/148‑1 (MQ=255)
                                                            aCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATccc                                                                                 >  2:853918/1‑148 (MQ=255)
                                                            aCTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATccc                                                                                 >  2:733423/1‑148 (MQ=255)
                                                                   ttttCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCcaca                                                                          <  1:291055/148‑1 (MQ=255)
                                                                    tttCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACaa                                                                         >  1:3262982/1‑148 (MQ=255)
                                                                       cTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATgg                                                                      <  2:457254/148‑1 (MQ=255)
                                                                        tGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGc                                                                     <  1:2415763/148‑1 (MQ=255)
                                                                         gAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGcc                                                                    <  2:2273150/148‑1 (MQ=255)
                                                                         gAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGcc                                                                    >  1:253773/1‑148 (MQ=255)
                                                                          aaGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGccc                                                                   >  1:870626/1‑148 (MQ=255)
                                                                            gagGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGt                                                                 >  2:2612553/1‑148 (MQ=255)
                                                                            gagGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGt                                                                 >  1:2270344/1‑148 (MQ=255)
                                                                            gagGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGt                                                                 <  1:1223840/148‑1 (MQ=255)
                                                                            gagGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGt                                                                 <  1:108169/148‑1 (MQ=255)
                                                                             agGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGtt                                                                <  2:2009963/148‑1 (MQ=255)
                                                                             agGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGtt                                                                <  2:210586/148‑1 (MQ=255)
                                                                                 cGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCa                                                            <  1:2664297/148‑1 (MQ=255)
                                                                                        gTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGa                                                     <  1:2918908/148‑1 (MQ=255)
                                                                                        gTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGa                                                     >  1:1355360/1‑148 (MQ=255)
                                                                                           aCTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaata                                                  <  1:531034/148‑1 (MQ=255)
                                                                                           aCTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaata                                                  <  1:469301/148‑1 (MQ=255)
                                                                                           aCTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaata                                                  <  1:2565919/148‑1 (MQ=255)
                                                                                            cTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaataa                                                 <  2:2076184/148‑1 (MQ=255)
                                                                                             ttCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGaataat                                                >  2:2761933/1‑148 (MQ=255)
                                                                                              tCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAatt                                               <  2:2394496/148‑1 (MQ=255)
                                                                                               cAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAatta                                              <  1:100365/148‑1 (MQ=255)
                                                                                                  aaGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAAttatta                                           <  2:1995080/148‑1 (MQ=255)
                                                                                                    ggCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATg                                         >  2:297572/1‑148 (MQ=255)
                                                                                                    ggCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATg                                         >  1:2437703/1‑148 (MQ=255)
                                                                                                        cATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTg                                     >  1:2243906/1‑148 (MQ=255)
                                                                                                         aTCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTga                                    >  2:3188871/1‑148 (MQ=255)
                                                                                                          tCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTgag                                   >  1:612392/1‑148 (MQ=255)
                                                                                                               tttCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAAttt                              >  2:2242524/1‑148 (MQ=255)
                                                                                                               tttCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAAttt                              <  1:419273/148‑1 (MQ=255)
                                                                                                               tttCACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAAttt                              <  1:3107617/148‑1 (MQ=255)
                                                                                                                  cACCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACg                           <  2:1627311/148‑1 (MQ=255)
                                                                                                                   aCCGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGt                          <  1:120489/148‑1 (MQ=255)
                                                                                                                    ccGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGtt                         <  2:1529575/148‑1 (MQ=255)
                                                                                                                     cGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTa                        <  1:1909786/148‑1 (MQ=255)
                                                                                                                     cGCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTa                        <  1:3136173/148‑1 (MQ=255)
                                                                                                                      gCTGTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTAt                       >  1:1777403/1‑148 (MQ=255)
                                                                                                                         gTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTa                    <  2:1199359/148‑1 (MQ=255)
                                                                                                                         gTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTa                    <  2:118292/148‑1 (MQ=255)
                                                                                                                         gTTGGCGCGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTa                    >  1:2922218/1‑148 (MQ=255)
                                                                                                                             gcgcGTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTc                >  1:675169/1‑148 (MQ=255)
                                                                                                                                cgTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGt             >  1:2402345/1‑148 (MQ=255)
                                                                                                                                cgTTATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGt             >  2:2520522/1‑148 (MQ=255)
                                                                                                                                  ttATCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTAt           <  1:2027398/148‑1 (MQ=255)
                                                                                                                                    aTCGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCa         >  2:1019917/1‑148 (MQ=255)
                                                                                                                                      cGTTGTTGTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACg       >  2:1553815/1‑148 (MQ=255)
                                                                                                                                        ttgttgTAGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGcc     >  1:1638862/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAACTTTTTCATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCAT  >  CP009273/2645449‑2645735

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: