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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,678,0300CA10.6% 855.4 / 51.9 244A279A (GCG→GCTglyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (103/115);  new base A (19/7);  total (122/122)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 248 total aligned reads displayed.

GTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGT  >  CP009273/2677881‑2678179
                                                                                                                                                     |                                                                                                                                                     
  cAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGc                                                                                                                                                     <  2:1604441/150‑1 (MQ=255)
   aGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGct                                                                                                                                                    <  1:2206981/150‑1 (MQ=255)
     gTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGcttc                                                                                                                                                  <  2:2569986/150‑1 (MQ=255)
      tttttATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGcttct                                                                                                                                                 <  2:2726099/150‑1 (MQ=255)
       ttttATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGcttctt                                                                                                                                                >  2:2802894/1‑150 (MQ=255)
          tATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCa                                                                                                                                             >  2:1929664/1‑150 (MQ=255)
           atcaaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCag                                                                                                                                            >  1:4187623/1‑150 (MQ=255)
              aaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGc                                                                                                                                         >  1:2856419/1‑150 (MQ=255)
               accagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCa                                                                                                                                        >  2:1973718/1‑150 (MQ=255)
                  agatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAAcc                                                                                                                                     >  1:3438493/1‑150 (MQ=255)
                    atcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGc                                                                                                                                   >  1:1772891/1‑150 (MQ=255)
                    atcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGc                                                                                                                                   >  1:4298554/1‑150 (MQ=255)
                     tcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCt                                                                                                                                  >  2:4077035/1‑150 (MQ=255)
                       aaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCttt                                                                                                                                <  1:3221065/150‑1 (MQ=255)
                       aaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCttt                                                                                                                                <  1:3363131/150‑1 (MQ=255)
                       aaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCttt                                                                                                                                >  1:2435603/1‑150 (MQ=255)
                        accagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTa                                                                                                                               >  1:752180/1‑150 (MQ=255)
                         ccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTAc                                                                                                                              >  1:4168917/1‑150 (MQ=255)
                           agGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCg                                                                                                                            >  1:1917951/1‑150 (MQ=255)
                           agGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCg                                                                                                                            >  1:1747455/1‑150 (MQ=255)
                            ggAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCgg                                                                                                                           <  1:1340667/150‑1 (MQ=255)
                             gAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTCGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGc                                                                                                                          <  2:2494214/150‑1 (MQ=255)
                               aCAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGa                                                                                                                        >  1:2236403/1‑150 (MQ=255)
                                cAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAt                                                                                                                       >  2:1305490/1‑150 (MQ=255)
                                cAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAt                                                                                                                       <  1:1281179/150‑1 (MQ=255)
                                  ggtggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTa                                                                                                                     <  2:271314/150‑1 (MQ=255)
                                   gtggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTAc                                                                                                                    <  1:1118463/150‑1 (MQ=255)
                                     ggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGt                                                                                                                  >  2:1788604/1‑150 (MQ=255)
                                      gtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTg                                                                                                                 >  1:639591/1‑150 (MQ=255)
                                         aTCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCAt                                                                                                              >  1:112609/1‑150 (MQ=255)
                                          tCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATc                                                                                                             <  1:2547396/150‑1 (MQ=255)
                                          tCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATc                                                                                                             <  1:1510286/150‑1 (MQ=255)
                                            aGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAGTACGTGCATCaa                                                                                                           >  1:584735/1‑150 (MQ=255)
                                                  gccgGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGAcc                                                                                                     <  1:2332049/150‑1 (MQ=255)
                                                      gAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGccc                                                                                                 <  1:2946894/150‑1 (MQ=255)
                                                       aaaCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCt                                                                                                <  1:3159301/150‑1 (MQ=255)
                                                         aCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGa                                                                                              <  2:3966198/150‑1 (MQ=255)
                                                             cTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGaccac                                                                                          >  2:3979745/1‑150 (MQ=255)
                                                               ttGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGAccacca                                                                                        <  2:875630/150‑1 (MQ=255)
                                                               tggTAGCCGCGCTCGAGGCACACTTCTACCCTCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCACAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGAccacca                                                                                        <  2:3518449/148‑1 (MQ=255)
                                                                  tAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAGTACGTGCATCAACGGACCGCCCTGACCACCAggg                                                                                     >  2:2428057/1‑150 (MQ=255)
                                                                   aGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGa                                                                                    >  2:2636355/1‑150 (MQ=255)
                                                                   aGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGa                                                                                    >  2:1894595/1‑150 (MQ=255)
                                                                      cgcgCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGaaaa                                                                                 >  1:2343386/1‑150 (MQ=255)
                                                                       gcgcTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAAc                                                                                <  2:1878044/150‑1 (MQ=255)
                                                                        cgcTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACg                                                                               <  1:3940259/150‑1 (MQ=255)
                                                                           tCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCa                                                                            <  1:4076668/150‑1 (MQ=255)
                                                                            cGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCag                                                                           >  1:2066658/1‑150 (MQ=255)
                                                                            cGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCag                                                                           >  1:2067453/1‑150 (MQ=255)
                                                                             gAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGACAACGGCaga                                                                          >  1:1195153/1‑150 (MQ=255)
                                                                                gAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTTCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGtt                                                                       <  2:2963263/150‑1 (MQ=255)
                                                                                gAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGtt                                                                       <  2:816173/150‑1 (MQ=255)
                                                                                gAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTCTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGtt                                                                       >  1:4450279/1‑150 (MQ=255)
                                                                                 aaCACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACAGCCCTGACCACCAGGGAAAACGGCAGAGTTc                                                                      <  2:2291849/150‑1 (MQ=255)
                                                                                  acacTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCa                                                                     >  2:3021592/1‑150 (MQ=255)
                                                                                   cacTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTAGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAg                                                                    >  1:2592676/1‑150 (MQ=255)
                                                                                     cTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtt                                                                  >  2:3656407/1‑150 (MQ=255)
                                                                                        cTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGttttt                                                               <  2:1422002/150‑1 (MQ=255)
                                                                                         tACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtttttt                                                              <  2:773074/150‑1 (MQ=255)
                                                                                          aCCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTg                                                             >  2:693532/1‑150 (MQ=255)
                                                                                           ccATCGCTTTTGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTGAGAGCAATCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGt                                                            <  1:2941671/150‑1 (MQ=255)
                                                                                             aTCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTAc                                                          >  2:4463900/1‑150 (MQ=255)
                                                                                                 cTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGct                                                      <  2:2848786/150‑1 (MQ=255)
                                                                                                    tAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCtt                                                   <  2:3852305/150‑1 (MQ=255)
                                                                                                     aGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTc                                                  <  2:847352/150‑1 (MQ=255)
                                                                                                     aGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTc                                                  >  2:3208481/1‑150 (MQ=255)
                                                                                                        gTTTTTAGCGACCCGCTGCTGGTAAGTATTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCt                                               <  2:4262318/150‑1 (MQ=255)
                                                                                                            ttAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTacca                                           >  1:3549931/1‑150 (MQ=255)
                                                                                                              aGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTaccacc                                         >  1:390587/1‑150 (MQ=255)
                                                                                                              aGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTaccacc                                         <  2:2078111/150‑1 (MQ=255)
                                                                                                              aGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTCACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTaccacc                                         <  2:2078173/150‑1 (MQ=255)
                                                                                                               gCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCt                                        <  1:3631072/150‑1 (MQ=255)
                                                                                                               gCGACCTGCTGCTGGTAAGTTATGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCt                                        <  2:3300270/150‑1 (MQ=255)
                                                                                                                 gACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCttt                                      >  1:49595/1‑150 (MQ=255)
                                                                                                                 gACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCttt                                      >  1:49674/1‑150 (MQ=255)
                                                                                                                 gACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCttt                                      <  2:1747498/150‑1 (MQ=255)
                                                                                                                     tgctgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGcc                                  <  2:3761841/150‑1 (MQ=255)
                                                                                                                      gctgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCa                                 <  2:1597408/150‑1 (MQ=255)
                                                                                                                        tgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAgg                               >  1:3730341/1‑150 (MQ=255)
                                                                                                                          cttgTAAGTTTTGAACTCAGGCTCGATCGCTTCTTTCAGAGCAACCGCTTTACAGGCGATTTAGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGTCATAGTTCAGTTTTTTTTACAGCTATTCGCTACCACCTTTAGCCAGGAt                             <  2:1130291/147‑1 (MQ=255)
                                                                                                                          ctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGAt                             <  1:1788620/150‑1 (MQ=255)
                                                                                                                           tgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACAAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATc                            >  1:4467157/1‑150 (MQ=255)
                                                                                                                             gTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAg                          >  2:2624365/1‑150 (MQ=255)
                                                                                                                              tAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAgg                         >  1:2154298/1‑150 (MQ=255)
                                                                                                                                aGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGgcc                       <  1:3716495/150‑1 (MQ=255)
                                                                                                                                aGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGgcc                       <  1:1451359/150‑1 (MQ=255)
                                                                                                                                     tGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGTATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCcgccgc                  >  1:952671/1‑150 (MQ=255)
                                                                                                                                       aaCTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTGTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCcgcgc                <  1:3718687/150‑1 (MQ=255)
                                                                                                                                        aCTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCcgcgcg               <  2:3190254/150‑1 (MQ=255)
                                                                                                                                         cTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCgg              >  2:4472199/1‑150 (MQ=255)
                                                                                                                                          tCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGa             >  1:3187135/1‑150 (MQ=255)
                                                                                                                                           cAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGAc            >  1:2917868/1‑150 (MQ=255)
                                                                                                                                           cAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGAc            <  2:977545/150‑1 (MQ=255)
                                                                                                                                                tCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGcc       <  2:2236403/150‑1 (MQ=255)
                                                                                                                                                tCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGcc       <  1:3979745/150‑1 (MQ=255)
                                                                                                                                                 ccATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCa      <  1:80884/150‑1 (MQ=255)
                                                                                                                                                  cATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAg     >  2:1090378/1‑150 (MQ=255)
                                                                                                                                                   aTCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAgg    <  2:2412118/150‑1 (MQ=255)
                                                                                                                                                    tCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAggg   >  2:3350478/1‑150 (MQ=255)
                                                                                                                                                     cGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGt  <  2:2815115/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGT  >  CP009273/2677881‑2678179

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: