Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,218,022 C→A 100% E110* (GAA→TAA)  pbpG ← D‑alanyl‑D‑alanine endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,218,0220CA97.9% 594.1 / 3.6 188E110* (GAA→TAA) pbpGD‑alanyl‑D‑alanine endopeptidase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (99/85);  total (101/87)
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.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 188 total aligned reads displayed.

TTGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCAC  >  CP009273/2217877‑2218166
                                                                                                                                                 |                                                                                                                                                
ttgttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGc                                                                                                                                                >  2:2381853/1‑148 (MQ=255)
ttgttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGc                                                                                                                                                <  1:2563017/148‑1 (MQ=255)
 tgttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCt                                                                                                                                               >  2:1984615/1‑148 (MQ=255)
 tgttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCt                                                                                                                                               >  1:1656694/1‑148 (MQ=255)
   ttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTAt                                                                                                                                             <  2:2982407/148‑1 (MQ=255)
    tCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTAtt                                                                                                                                            >  2:2877285/1‑148 (MQ=255)
     cATTCCGAGCGATTTCGCTTTCGCATTCAGTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTc                                                                                                                                           <  2:3460716/148‑1 (MQ=255)
      aTTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCa                                                                                                                                          <  1:374941/148‑1 (MQ=255)
        tCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGt                                                                                                                                        >  2:905999/1‑148 (MQ=255)
        tCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGt                                                                                                                                        >  1:3195703/1‑148 (MQ=255)
         ccGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTc                                                                                                                                       <  1:1869617/148‑1 (MQ=255)
         ccGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTc                                                                                                                                       <  1:1869320/148‑1 (MQ=255)
           gAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGt                                                                                                                                     <  1:1014518/148‑1 (MQ=255)
           gAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGt                                                                                                                                     <  1:1014457/148‑1 (MQ=255)
            aGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTa                                                                                                                                    >  2:2138686/1‑148 (MQ=255)
             gCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAc                                                                                                                                   >  2:770343/1‑148 (MQ=255)
                    cgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGaa                                                                                                                            <  1:2735879/148‑1 (MQ=255)
                    cgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGaa                                                                                                                            >  2:1743350/1‑148 (MQ=255)
                      ctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATa                                                                                                                          <  1:916213/148‑1 (MQ=255)
                        ttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGa                                                                                                                        >  1:643025/1‑148 (MQ=255)
                            cattcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGAcccc                                                                                                                    >  1:1833154/1‑148 (MQ=255)
                             attcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTGAGCTATTCAGTCGTACTCGCGAATAGACCCCt                                                                                                                   >  2:2275678/1‑148 (MQ=255)
                                   tGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtc                                                                                                             <  1:3326043/148‑1 (MQ=255)
                                       ttAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCg                                                                                                         >  2:2687027/1‑148 (MQ=255)
                                        taataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCgg                                                                                                        >  2:1411577/1‑148 (MQ=255)
                                          ataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCg                                                                                                      <  1:435140/148‑1 (MQ=255)
                                           taaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGt                                                                                                     >  2:2113076/1‑148 (MQ=255)
                                            aaaGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTc                                                                                                    <  1:2195820/148‑1 (MQ=255)
                                            aaaGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTc                                                                                                    <  2:3495071/148‑1 (MQ=255)
                                             aaGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCt                                                                                                   >  1:1517195/1‑148 (MQ=255)
                                                gCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGc                                                                                                >  1:3391015/1‑148 (MQ=255)
                                                gCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGc                                                                                                >  2:2245176/1‑148 (MQ=255)
                                                   ttGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGa                                                                                             <  2:2330293/148‑1 (MQ=255)
                                                     gTAACCACCGGTATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGTAGACATCAGCGCCAGCTACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGata                                                                                           <  2:1432787/148‑1 (MQ=255)
                                                     gTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGata                                                                                           <  2:1385242/148‑1 (MQ=255)
                                                      tAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGatat                                                                                          <  2:2290803/148‑1 (MQ=255)
                                                       aaCCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATc                                                                                         >  1:387302/1‑148 (MQ=255)
                                                        accaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATcc                                                                                        >  2:3040651/1‑148 (MQ=255)
                                                        accaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATcc                                                                                        <  2:2538258/148‑1 (MQ=255)
                                                        accaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATcc                                                                                        >  2:3040640/1‑148 (MQ=255)
                                                         ccaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCa                                                                                       >  1:3150714/1‑148 (MQ=255)
                                                         ccaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCa                                                                                       >  1:2704233/1‑148 (MQ=255)
                                                            ccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTCTTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACtt                                                                                    >  2:1392740/1‑148 (MQ=255)
                                                            ccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACtt                                                                                    <  1:2781971/148‑1 (MQ=255)
                                                            ccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACtt                                                                                    >  2:807827/1‑148 (MQ=255)
                                                               ggATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGACTAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTa                                                                                 >  2:2617063/1‑148 (MQ=255)
                                                               ggATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTa                                                                                 >  2:2726135/1‑148 (MQ=255)
                                                                gATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAg                                                                                <  1:1594479/148‑1 (MQ=255)
                                                                gATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAg                                                                                >  2:378407/1‑148 (MQ=255)
                                                                 aTAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGt                                                                               >  1:1927907/1‑148 (MQ=255)
                                                                 aTAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGt                                                                               <  1:1976283/148‑1 (MQ=255)
                                                                   aaTGGTGCGAAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttt                                                                             <  2:51792/148‑1 (MQ=255)
                                                                    aTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGtttt                                                                            <  1:1067410/148‑1 (MQ=255)
                                                                     tGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                           <  2:3239103/148‑1 (MQ=255)
                                                                     tGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                           >  2:1180315/1‑148 (MQ=255)
                                                                     tGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                           >  1:3359161/1‑148 (MQ=255)
                                                                       gTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCa                                                                         >  1:1576161/1‑148 (MQ=255)
                                                                       gTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCa                                                                         <  2:1027894/148‑1 (MQ=255)
                                                                             aaGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTGTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAga                                                                   >  2:2045529/1‑147 (MQ=255)
                                                                                   tGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGt                                                             <  2:3386828/148‑1 (MQ=255)
                                                                                   tGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGt                                                             <  1:770343/148‑1 (MQ=255)
                                                                                     ccgccgCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCg                                                           >  1:1480511/1‑148 (MQ=255)
                                                                                       gccgcGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTg                                                         <  1:3034198/148‑1 (MQ=255)
                                                                                           cgcgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATc                                                     <  2:1561179/148‑1 (MQ=255)
                                                                                                ttttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGca                                                >  1:2969141/1‑148 (MQ=255)
                                                                                                ttttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGca                                                >  1:2969223/1‑148 (MQ=255)
                                                                                                  ttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGcaca                                              >  1:2440159/1‑148 (MQ=255)
                                                                                                    cTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAAc                                            <  1:991247/148‑1 (MQ=255)
                                                                                                      gAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCa                                          <  2:2572975/148‑1 (MQ=255)
                                                                                                       aaGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCAt                                         <  2:2490787/148‑1 (MQ=255)
                                                                                                       aaGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCAt                                         <  2:3218340/148‑1 (MQ=255)
                                                                                                         gACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATcg                                       <  1:2643168/148‑1 (MQ=255)
                                                                                                          aCATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATcgc                                      >  1:2641301/1‑148 (MQ=255)
                                                                                                           cATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATcgcg                                     <  2:965629/148‑1 (MQ=255)
                                                                                                             tCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGt                                   >  1:3146794/1‑148 (MQ=255)
                                                                                                              cAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTc                                  >  2:2168968/1‑148 (MQ=255)
                                                                                                                gcgcCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCAt                                <  2:2781340/148‑1 (MQ=255)
                                                                                                                   ccAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTaa                             <  2:1765665/148‑1 (MQ=255)
                                                                                                                          cAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTg                      <  1:3316803/148‑1 (MQ=255)
                                                                                                                           aGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAAACATCGCGGTCATTAATTTGCTGa                     <  1:2113076/148‑1 (MQ=255)
                                                                                                                             cATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATa                   >  2:3135841/1‑148 (MQ=255)
                                                                                                                             cATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATa                   >  1:2155813/1‑148 (MQ=255)
                                                                                                                             cATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATAACCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAAATTGCTGATa                   >  2:1663846/1‑148 (MQ=255)
                                                                                                                              atatCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAg                  <  1:905999/148‑1 (MQ=255)
                                                                                                                                  cTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGc              >  1:1157181/1‑148 (MQ=255)
                                                                                                                                  cTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGc              >  1:3281588/1‑148 (MQ=255)
                                                                                                                                  cTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGc              >  1:1153577/1‑148 (MQ=255)
                                                                                                                                  cTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGc              >  1:2903740/1‑148 (MQ=255)
                                                                                                                                   tttACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCa             <  1:1342818/148‑1 (MQ=255)
                                                                                                                                   tttACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCa             <  1:148670/148‑1 (MQ=255)
                                                                                                                                    ttACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCaa            <  1:2610936/148‑1 (MQ=255)
                                                                                                                                    ttACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCaa            <  1:260086/148‑1 (MQ=255)
                                                                                                                                     tACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAAt           <  2:3104475/148‑1 (MQ=255)
                                                                                                                                       cGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCg         >  1:1098874/1‑148 (MQ=255)
                                                                                                                                          cTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGAc      >  1:1777040/1‑148 (MQ=255)
                                                                                                                                          cTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGAc      >  1:693905/1‑148 (MQ=255)
                                                                                                                                          cTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGAc      >  2:1610371/1‑148 (MQ=255)
                                                                                                                                            gATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGc    <  1:1043090/148‑1 (MQ=255)
                                                                                                                                             aTTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCa   >  2:3508708/1‑148 (MQ=255)
                                                                                                                                              tttAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCAc  >  2:469896/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCAC  >  CP009273/2217877‑2218166

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