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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,678,0300CA 7.3% 932.2 / 31.6 244A279A (GCG→GCTglyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (113/113);  new base A (4/14);  total (117/127)
Fisher's exact test for biased strand distribution p-value = 2.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 252 total aligned reads displayed.

GTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGT  >  CP009273/2677881‑2678179
                                                                                                                                                     |                                                                                                                                                     
gTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATc                                                                                                                                                       >  1:2252270/1‑150 (MQ=255)
gTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACATGCTGCTGGTAAGTTTTGAACTCAGGCTCCATc                                                                                                                                                       >  2:3647723/1‑150 (MQ=255)
gTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACATGCTGCTGGTAAGTTTTGAACTCAGGCTCCATc                                                                                                                                                       >  2:4510580/1‑150 (MQ=255)
  cAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGc                                                                                                                                                     >  2:1879827/1‑150 (MQ=255)
   aGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGct                                                                                                                                                    >  1:928270/1‑150 (MQ=255)
    ggTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGctt                                                                                                                                                   >  1:2634101/1‑150 (MQ=255)
          tATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCAGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGAGACCTGATGCTGGTAAGTTTTGAACTCAGGCTCCATCTCTTCTTTCa                                                                                                                                             >  1:1175834/1‑150 (MQ=255)
            tcaaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCaga                                                                                                                                           <  2:3685143/150‑1 (MQ=255)
              aaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGc                                                                                                                                         >  2:2833403/1‑150 (MQ=255)
                   gatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCg                                                                                                                                    <  2:3454040/150‑1 (MQ=255)
                    atcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTAAGCGTTTTTAGCGACCTGCTTCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGc                                                                                                                                   >  1:346432/1‑150 (MQ=255)
                     tcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCt                                                                                                                                  <  2:2047/150‑1 (MQ=255)
                     tcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCt                                                                                                                                  <  1:764194/150‑1 (MQ=255)
                       aaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCttt                                                                                                                                >  2:2851052/1‑150 (MQ=255)
                         ccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTAc                                                                                                                              >  1:1420730/1‑150 (MQ=255)
                          cagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTAcc                                                                                                                             <  1:3246960/150‑1 (MQ=255)
                          cagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTAcc                                                                                                                             >  2:3829907/1‑150 (MQ=255)
                              aaCAGGTGGTTATCAGTGCCGCCGGAAACCTCTTTGTAGACGCGCTCGATGAACACTTCTTCCAGCGCTTTAGCGTTTTTATCGACCTGCTGCTTGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCTCTTTACCGGCg                                                                                                                         <  2:2728735/150‑1 (MQ=255)
                                cAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAt                                                                                                                       >  2:1803529/1‑150 (MQ=255)
                                 aGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAtt                                                                                                                      >  1:1827005/1‑150 (MQ=255)
                                 aGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGAtt                                                                                                                      <  2:4241715/150‑1 (MQ=255)
                                  ggtggtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTa                                                                                                                     <  1:994547/150‑1 (MQ=255)
                                       ttATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCACGAGGAACAATTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCATGCTACATCGCTTCTTTCAGAGCAACCGCTTTACCCGCGATTACGTGc                                                                                                                >  1:533839/1‑150 (MQ=255)
                                            aGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCaa                                                                                                           >  2:1932160/1‑150 (MQ=255)
                                            aGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCaa                                                                                                           >  2:1566325/1‑150 (MQ=255)
                                             gTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAAc                                                                                                          >  2:1521644/1‑150 (MQ=255)
                                               gccgccGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACgg                                                                                                        <  2:3083841/150‑1 (MQ=255)
                                                  gccgGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGAcc                                                                                                     >  1:2781548/1‑150 (MQ=255)
                                                   ccgGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCg                                                                                                    >  1:1301114/1‑150 (MQ=255)
                                                    cgGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGc                                                                                                   >  2:3514822/1‑150 (MQ=255)
                                                    cgGAAACCACTTTGTACACGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCAGCTGGTAATATTTGAACTCAGGCTCAATCGCTTCTTTAAGAGCAAACTCTATACCGGCGATTACGTGCATCAACGGACCGc                                                                                                   >  2:4323628/1‑150 (MQ=255)
                                                     ggAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGcc                                                                                                  >  2:4478631/1‑150 (MQ=255)
                                                      gAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGccc                                                                                                 >  2:2721872/1‑150 (MQ=255)
                                                      gAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGccc                                                                                                 <  2:2352259/150‑1 (MQ=255)
                                                       aaaCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCt                                                                                                <  2:4066034/150‑1 (MQ=255)
                                                       aaaCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCt                                                                                                >  1:4130989/1‑150 (MQ=255)
                                                       aaaCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCt                                                                                                >  1:3176312/1‑150 (MQ=255)
                                                        aaCCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTg                                                                                               <  1:1005648/150‑1 (MQ=255)
                                                            aCTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGacca                                                                                           >  2:124520/1‑150 (MQ=255)
                                                               ttGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGAccacca                                                                                        <  2:4624773/150‑1 (MQ=255)
                                                                  tAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAggg                                                                                     >  2:1465149/1‑150 (MQ=255)
                                                                   aGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGa                                                                                    >  1:2859441/1‑150 (MQ=255)
                                                                    gCCGCGCTCGAGGAACACTTCTACCATCGCTTTATCGTTTTTAGCGACCTGCTTCTGGTAAGTTTTGAACTCATGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCTCCCTGACCACCAGGGaa                                                                                   <  1:1879827/150‑1 (MQ=255)
                                                                    gCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGaa                                                                                   <  2:638086/150‑1 (MQ=255)
                                                                    gCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGaa                                                                                   >  1:4357402/1‑150 (MQ=255)
                                                                    gCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGaa                                                                                   <  2:1153496/150‑1 (MQ=255)
                                                                      cgcgCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGaaaa                                                                                 <  1:642841/150‑1 (MQ=255)
                                                                       gcgcTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAAc                                                                                >  2:2931423/1‑150 (MQ=255)
                                                                        cgcTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACg                                                                               <  2:1703726/150‑1 (MQ=255)
                                                                            cGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCag                                                                           <  1:3803655/150‑1 (MQ=255)
                                                                             gAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCaga                                                                          >  1:819563/1‑150 (MQ=255)
                                                                                     cTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtt                                                                  <  2:2966990/150‑1 (MQ=255)
                                                                                     cTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtt                                                                  <  1:1934822/150‑1 (MQ=255)
                                                                                      ttCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGttt                                                                 >  2:3364425/1‑150 (MQ=255)
                                                                                         tACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtttttt                                                              <  1:3352044/150‑1 (MQ=255)
                                                                                         tACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtttttt                                                              >  1:1370498/1‑150 (MQ=255)
                                                                                         tACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtttttt                                                              <  1:1760664/150‑1 (MQ=255)
                                                                                            cATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTa                                                           <  2:1480387/150‑1 (MQ=255)
                                                                                            cATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTa                                                           >  2:535121/1‑150 (MQ=255)
                                                                                             aTCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTAc                                                          >  1:4475676/1‑150 (MQ=255)
                                                                                             aTCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTAc                                                          <  1:2163142/150‑1 (MQ=255)
                                                                                               cGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGACCTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAg                                                        <  1:811729/150‑1 (MQ=255)
                                                                                               cGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAg                                                        <  2:1792997/150‑1 (MQ=255)
                                                                                                gCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGc                                                       <  2:4130827/150‑1 (MQ=255)
                                                                                                gCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGc                                                       <  2:4293797/150‑1 (MQ=255)
                                                                                                 cTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGct                                                      <  1:4191779/150‑1 (MQ=255)
                                                                                                   ttAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtttttt                                                              >  2:478268/1‑140 (MQ=255)
                                                                                                   ttAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGtttttt                                                              <  1:478268/140‑1 (MQ=255)
                                                                                                   ttAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGctct                                                    <  1:2601095/150‑1 (MQ=255)
                                                                                                    tAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCtt                                                   <  2:2245743/150‑1 (MQ=255)
                                                                                                     aGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTc                                                  <  1:1675674/150‑1 (MQ=255)
                                                                                                     aGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTc                                                  <  1:3946941/150‑1 (MQ=255)
                                                                                                     aGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATAGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTc                                                  <  2:2524930/150‑1 (MQ=255)
                                                                                                      gCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCg                                                 <  1:4266117/150‑1 (MQ=255)
                                                                                                        gTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCt                                               <  1:3062814/150‑1 (MQ=255)
                                                                                                           tttAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTacc                                            <  2:4537969/150‑1 (MQ=255)
                                                                                                             tAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTaccac                                          >  1:2323013/1‑150 (MQ=255)
                                                                                                             tAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGNTCAGTTTTTTGTACAGCTCTTCGCTaccac                                          <  1:4281677/150‑1 (MQ=255)
                                                                                                              aGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTaccacc                                         >  2:2159492/1‑150 (MQ=255)
                                                                                                               gCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCt                                        <  1:629346/150‑1 (MQ=255)
                                                                                                               gCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCt                                        >  2:113741/1‑150 (MQ=255)
                                                                                                               gCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCt                                        <  2:2634101/150‑1 (MQ=255)
                                                                                                                cGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCtt                                       <  1:3248955/150‑1 (MQ=255)
                                                                                                                    ctgctgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGc                                   <  1:4121159/150‑1 (MQ=255)
                                                                                                                     tgctgctgGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGcc                                  >  2:1078593/1‑150 (MQ=255)
                                                                                                                             gTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAg                          <  2:1208213/150‑1 (MQ=255)
                                                                                                                                aGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGgcc                       <  1:106817/150‑1 (MQ=255)
                                                                                                                                aGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGgcc                       >  1:4222481/1‑150 (MQ=255)
                                                                                                                                aGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGgcc                       <  1:1662942/150‑1 (MQ=255)
                                                                                                                                    ttGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGccgccg                   >  1:4601593/1‑150 (MQ=255)
                                                                                                                                        aCTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTTTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCcgcgcg               >  1:1999198/1‑150 (MQ=255)
                                                                                                                                         cTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCgg              <  1:3621037/150‑1 (MQ=255)
                                                                                                                                             ggCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGAccc          <  1:3381605/150‑1 (MQ=255)
                                                                                                                                                tCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGcc       <  1:2931423/150‑1 (MQ=255)
                                                                                                                                                 ccATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCa      >  2:1383276/1‑150 (MQ=255)
                                                                                                                                                  cATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAg     >  2:253219/1‑150 (MQ=255)
                                                                                                                                                  cATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAg     >  2:585817/1‑150 (MQ=255)
                                                                                                                                                  cATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAg     <  2:928270/150‑1 (MQ=255)
                                                                                                                                                  cATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGGTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAg     >  1:2068108/1‑150 (MQ=255)
                                                                                                                                                    tCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACCGCCCTGACCACCAGGGAAAACGGCAGAGTTCAGTTTTTTGTACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAggg   >  1:2430685/1‑150 (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: