Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,678,030 C→A 11.2% A279A (GCG→GCT glyA ← serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,678,0300CA11.2% 661.5 / 39.7 195A279A (GCG→GCTglyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (99/74);  new base A (13/9);  total (112/83)
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.95e-01

Only 100 of 198 total aligned reads displayed.

CAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGG  >  CP009273/2677883‑2678177
                                                                                                                                                   |                                                                                                                                                   
cAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATc                                                                                                                                                     >  2:162375/1‑148 (MQ=255)
   gTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGct                                                                                                                                                  >  1:1526833/1‑148 (MQ=255)
    tttttATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGctt                                                                                                                                                 <  1:814559/148‑1 (MQ=255)
    tttttATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGctt                                                                                                                                                 >  2:460281/1‑148 (MQ=255)
    tttttATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGctt                                                                                                                                                 <  2:1403003/148‑1 (MQ=255)
    tttttATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGACGGAAACCACTTTGTAGCCGCGCTCGAGGAAAACTTCTACCATCGCGTTAGCGTTTTTAGAGACCTGCTGATGGTAAGTTTTGAACTCAGGCTCCATCGctt                                                                                                                                                 >  2:73549/1‑148 (MQ=255)
     ttttATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGcttc                                                                                                                                                <  2:709206/148‑1 (MQ=255)
        tATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCTCTTCttt                                                                                                                                             <  2:2970109/148‑1 (MQ=255)
           caaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCag                                                                                                                                          >  1:2498979/1‑148 (MQ=255)
               cagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCa                                                                                                                                      <  1:2907219/148‑1 (MQ=255)
               cagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCa                                                                                                                                      <  1:2906912/148‑1 (MQ=255)
                agatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCaa                                                                                                                                     <  2:43985/148‑1 (MQ=255)
                agatcaaccagGAACAGGTGGTTATAAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCaa                                                                                                                                     <  1:1623124/148‑1 (MQ=255)
                agatcaaccagGAACAGGTGGTTATAAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCaa                                                                                                                                     <  1:892684/148‑1 (MQ=255)
                agatcaaccagGAACAGGTGGTTATAAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCaa                                                                                                                                     <  1:892665/148‑1 (MQ=255)
                 gatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAAc                                                                                                                                    >  1:924633/1‑148 (MQ=255)
                   tcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCg                                                                                                                                  >  2:2039095/1‑148 (MQ=255)
                    caaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGc                                                                                                                                 <  2:2850023/148‑1 (MQ=255)
                    caaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTATAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGc                                                                                                                                 <  1:1393838/148‑1 (MQ=255)
                      accagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCagag                                                                                                                                        <  2:960834/139‑1 (MQ=255)
                        cagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTa                                                                                                                             <  2:564330/148‑1 (MQ=255)
                        cagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTa                                                                                                                             <  2:1661374/148‑1 (MQ=255)
                         agGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTAc                                                                                                                            <  2:818921/148‑1 (MQ=255)
                          ggAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTAcc                                                                                                                           >  1:1684108/1‑148 (MQ=255)
                          ggAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTAcc                                                                                                                           >  2:2691721/1‑148 (MQ=255)
                           gAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCg                                                                                                                          >  2:1016559/1‑148 (MQ=255)
                           gAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCg                                                                                                                          >  2:2063126/1‑148 (MQ=255)
                            aaCAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCgg                                                                                                                         >  1:2963148/1‑148 (MQ=255)
                            aaCAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCgg                                                                                                                         >  2:1794203/1‑148 (MQ=255)
                            aaCAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCgg                                                                                                                         >  2:1794130/1‑148 (MQ=255)
                             aCAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGc                                                                                                                        >  2:2885295/1‑148 (MQ=255)
                              cAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCg                                                                                                                       >  1:2382895/1‑148 (MQ=255)
                              cAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCg                                                                                                                       >  2:1067930/1‑148 (MQ=255)
                               aGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGa                                                                                                                      <  1:318625/148‑1 (MQ=255)
                                ggtggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAt                                                                                                                     >  2:2603087/1‑148 (MQ=255)
                                 gtggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAtt                                                                                                                    >  1:1455953/1‑148 (MQ=255)
                                  tggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTa                                                                                                                   >  1:1053068/1‑148 (MQ=255)
                                  tggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTa                                                                                                                   >  1:1053416/1‑148 (MQ=255)
                                      tATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTg                                                                                                               >  1:2445957/1‑148 (MQ=255)
                                       aTCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGc                                                                                                              <  1:113893/148‑1 (MQ=255)
                                       aTCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGc                                                                                                              <  1:113499/148‑1 (MQ=255)
                                        tCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCTGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCa                                                                                                             >  2:145644/1‑148 (MQ=255)
                                            tGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCaa                                                                                                         <  2:2377608/148‑1 (MQ=255)
                                            tGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCaa                                                                                                         >  1:2716332/1‑148 (MQ=255)
                                             gccgccGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTCTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAAc                                                                                                        <  1:2872878/148‑1 (MQ=255)
                                                gccgGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGa                                                                                                     <  1:2735125/148‑1 (MQ=255)
                                                gccgGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGa                                                                                                     <  1:1061643/148‑1 (MQ=255)
                                                     aaaCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGcc                                                                                                <  2:2458333/148‑1 (MQ=255)
                                                     aaaCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGcc                                                                                                <  1:450741/148‑1 (MQ=255)
                                                       aCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCt                                                                                              >  1:170194/1‑148 (MQ=255)
                                                        ccACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAAAGGACCGCCCTg                                                                                             >  2:351327/1‑148 (MQ=255)
                                                         cACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTGTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAGTACGTGCATCAACCGACCGCCCTGa                                                                                            >  1:1971985/1‑148 (MQ=255)
                                                             ttGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGaccac                                                                                        >  1:75868/1‑148 (MQ=255)
                                                             ttGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGaccac                                                                                        >  1:3267443/1‑148 (MQ=255)
                                                                  gCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCGGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCCTCCACGGCCCGCCCGGACCACCCggg                                                                                   >  1:585294/1‑148 (MQ=255)
                                                                     gcgcTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGaaa                                                                                <  2:1309682/148‑1 (MQ=255)
                                                                      cgcTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGaaaa                                                                               <  2:101935/148‑1 (MQ=255)
                                                                          cGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGc                                                                           >  1:3117148/1‑148 (MQ=255)
                                                                             ggAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCaga                                                                        >  1:48648/1‑148 (MQ=255)
                                                                              gAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCagag                                                                       >  2:2551265/1‑148 (MQ=255)
                                                                              gAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCagag                                                                       >  1:385458/1‑148 (MQ=255)
                                                                                    ttCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGt                                                                 >  2:309446/1‑148 (MQ=255)
                                                                                    ttCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGt                                                                 >  2:1817701/1‑148 (MQ=255)
                                                                                    ttCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGt                                                                 >  2:326941/1‑148 (MQ=255)
                                                                                     tCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGACCTCAGGCTCCATCGCTTCGTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtt                                                                <  2:2353497/148‑1 (MQ=255)
                                                                                     tCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCCTCAACGGACCGCCCTGACCACCAGGGAAAACGGCCGAGGTCAGtt                                                                >  1:1180676/1‑148 (MQ=255)
                                                                                      cTACCATCGCTTTCGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGttt                                                               <  1:2066833/148‑1 (MQ=255)
                                                                                       tACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTACGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtttt                                                              >  2:2364872/1‑148 (MQ=255)
                                                                                       tACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtttt                                                              <  1:2987871/148‑1 (MQ=255)
                                                                                       tACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCCGGGAAAACGGCAGAGTTCAGtttt                                                              >  2:2365634/1‑148 (MQ=255)
                                                                                               cTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAg                                                      <  1:975694/148‑1 (MQ=255)
                                                                                               cTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAg                                                      >  2:2205017/1‑148 (MQ=255)
                                                                                                  tAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGctc                                                   <  1:2220913/148‑1 (MQ=255)
                                                                                                     cGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTc                                                >  2:682465/1‑148 (MQ=255)
                                                                                                      gTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCg                                               <  1:1749215/148‑1 (MQ=255)
                                                                                                       tttttAGCTACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGTCCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGc                                              >  1:2671201/1‑148 (MQ=255)
                                                                                                       tttttAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGc                                              <  1:1236872/148‑1 (MQ=255)
                                                                                                       tttttAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGc                                              <  2:1744876/148‑1 (MQ=255)
                                                                                                         tttAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTa                                            >  1:1477763/1‑148 (MQ=255)
                                                                                                          ttAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTac                                           <  1:2603087/148‑1 (MQ=255)
                                                                                                             gCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTACTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCCGTGTTCAGTTTTTTGTACAGCTCTTCGATaccac                                        >  2:1560543/1‑148 (MQ=255)
                                                                                                              cGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTaccacc                                       >  1:2310418/1‑148 (MQ=255)
                                                                                                              cGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTATTCGCTaccacc                                       >  2:2594928/1‑148 (MQ=255)
                                                                                                                 ccTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCttt                                    >  1:2235185/1‑148 (MQ=255)
                                                                                                                  ctgctgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTc                                   >  1:2351708/1‑148 (MQ=255)
                                                                                                                  ctgctgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTc                                   >  1:960810/1‑148 (MQ=255)
                                                                                                                   tgctgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCg                                  >  2:480095/1‑148 (MQ=255)
                                                                                                                    gctgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGc                                 <  1:2039095/148‑1 (MQ=255)
                                                                                                                    gctgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGc                                 >  1:2430904/1‑148 (MQ=255)
                                                                                                                      tgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCa                               <  1:2063140/148‑1 (MQ=255)
                                                                                                                       gctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAg                              >  2:1279499/1‑148 (MQ=255)
                                                                                                                        ctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAgg                             <  1:1783386/148‑1 (MQ=255)
                                                                                                                           gTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATc                          <  1:1082251/148‑1 (MQ=255)
                                                                                                                            tAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCa                         >  1:2193423/1‑148 (MQ=255)
                                                                                                                              aGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAgg                       <  2:216527/148‑1 (MQ=255)
                                                                                                                              aGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAgg                       <  2:1053068/148‑1 (MQ=255)
                                                                                                                                         cAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGTCCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCgg            <  2:2671201/148‑1 (MQ=255)
                                                                                                                                         cAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCgg            >  2:2063468/1‑148 (MQ=255)
                                                                                                                                            gCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGAcc         <  1:1837743/148‑1 (MQ=255)
                                                                                                                                               ccATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGc      <  1:2205017/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGG  >  CP009273/2677883‑2678177

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: