Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,836,072 C→A 31.7% A109E (GCG→GAG)  nudG → CTP pyrophosphohydrolase; also hydrolyzes 2‑hydroxy‑dATP, 8‑hydroxy‑dGTP, 5‑hydroxy‑CTP, dCTP and 5‑methyl‑dCTP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,836,0720CA31.7% 390.9 / 181.0 208A109E (GCG→GAG) nudGCTP pyrophosphohydrolase; also hydrolyzes 2‑hydroxy‑dATP, 8‑hydroxy‑dGTP, 5‑hydroxy‑CTP, dCTP and 5‑methyl‑dCTP
Reads supporting (aligned to +/- strand):  ref base C (53/89);  new base A (35/31);  total (88/120)
Fisher's exact test for biased strand distribution p-value = 3.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 214 total aligned reads displayed.

GAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTT  >  CP009273/1835924‑1836220
                                                                                                                                                    |                                                                                                                                                    
gAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGgcg                                                                                                                                                     <  2:3152241/150‑1 (MQ=255)
  aGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCt                                                                                                                                                   <  2:1253812/150‑1 (MQ=255)
    cAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGc                                                                                                                                                 >  2:565489/1‑150 (MQ=255)
     aaCTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCa                                                                                                                                                >  2:422410/1‑150 (MQ=255)
        tgtgGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAata                                                                                                                                             >  1:2540930/1‑150 (MQ=255)
          tgGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATc                                                                                                                                           >  2:929679/1‑150 (MQ=255)
           gggTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATcc                                                                                                                                          <  1:3821610/150‑1 (MQ=255)
             gTGAATATATTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGc                                                                                                                                        >  2:2691845/1‑150 (MQ=255)
               gAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTg                                                                                                                                      >  1:1901185/1‑150 (MQ=255)
               gAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTg                                                                                                                                      >  2:124940/1‑150 (MQ=255)
                aaTATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTgg                                                                                                                                     <  2:2625372/150‑1 (MQ=255)
                 atatGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGc                                                                                                                                    <  1:2902192/150‑1 (MQ=255)
                 atatGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGc                                                                                                                                    <  1:2755019/150‑1 (MQ=255)
                   atGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGccc                                                                                                                                  >  2:183170/1‑150 (MQ=255)
                    tGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGcccc                                                                                                                                 >  1:3537678/1‑150 (MQ=255)
                     gTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCct                                                                                                                                <  1:3149482/150‑1 (MQ=255)
                       tGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCctgc                                                                                                                              >  2:3430453/1‑150 (MQ=255)
                       tGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCctgc                                                                                                                              <  2:3614787/150‑1 (MQ=255)
                        gccagccaTCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCctgct                                                                                                                             <  2:3361077/150‑1 (MQ=255)
                        gccagccaTCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCctgct                                                                                                                             <  1:3508717/150‑1 (MQ=255)
                         ccagccaTCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCctgctg                                                                                                                            <  2:3883145/150‑1 (MQ=255)
                          cagccaTCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGa                                                                                                                           >  2:892513/1‑150 (MQ=255)
                          cagccaTCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGa                                                                                                                           >  2:1261759/1‑150 (MQ=255)
                            gccaTCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACa                                                                                                                         >  2:2668799/1‑150 (MQ=255)
                            gccaTCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACa                                                                                                                         <  1:2119404/150‑1 (MQ=255)
                                 cAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCa                                                                                                                    <  1:3469803/150‑1 (MQ=255)
                                      agaAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCattatt                                                                                                               <  1:178203/150‑1 (MQ=255)
                                       gaAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCAttatta                                                                                                              >  1:71887/1‑150 (MQ=255)
                                       gaAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCAttatta                                                                                                              >  1:696693/1‑150 (MQ=255)
                                         aGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTaga                                                                                                            <  1:1274365/150‑1 (MQ=255)
                                           tttCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAgg                                                                                                          >  1:31108/1‑150 (MQ=255)
                                            ttCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGc                                                                                                         >  1:1539447/1‑150 (MQ=255)
                                            ttCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGc                                                                                                         <  1:3748392/150‑1 (MQ=255)
                                            ttCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGc                                                                                                         >  1:3633388/1‑150 (MQ=255)
                                                 ggCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTa                                                                                                    >  2:1179785/1‑150 (MQ=255)
                                                 ggCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTa                                                                                                    >  1:249885/1‑150 (MQ=255)
                                                  gCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTAt                                                                                                   <  2:2774706/150‑1 (MQ=255)
                                                   cGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATg                                                                                                  <  1:561044/150‑1 (MQ=255)
                                                    ggATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATgg                                                                                                 <  1:1077948/150‑1 (MQ=255)
                                                    ggATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATgg                                                                                                 >  1:2756804/1‑150 (MQ=255)
                                                       ttATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCGGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATGATTAGAGGCGTTTATGGCtt                                                                                              >  2:3838344/1‑150 (MQ=255)
                                                        tATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCttt                                                                                             <  1:2947581/150‑1 (MQ=255)
                                                        tATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCttt                                                                                             <  1:3912987/150‑1 (MQ=255)
                                                        tATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCttt                                                                                             >  1:3219697/1‑150 (MQ=255)
                                                         aTCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTa                                                                                            >  1:3844706/1‑150 (MQ=255)
                                                         aTCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTa                                                                                            >  1:810800/1‑150 (MQ=255)
                                                          tCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTAc                                                                                           <  1:1735724/150‑1 (MQ=255)
                                                             aTCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTAcgcg                                                                                        <  1:1739669/150‑1 (MQ=255)
                                                                ttCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGcgccg                                                                                     <  1:962052/150‑1 (MQ=255)
                                                                ttCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGcgccg                                                                                     <  2:2948525/150‑1 (MQ=255)
                                                                ttCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGcgccg                                                                                     <  2:2312623/150‑1 (MQ=255)
                                                                   aTGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCa                                                                                  >  2:617478/1‑150 (MQ=255)
                                                                   aTGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCa                                                                                  <  2:3044219/150‑1 (MQ=255)
                                                                        tGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCa                                                                             >  2:3364246/1‑150 (MQ=255)
                                                                           cACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCg                                                                          <  2:3454915/150‑1 (MQ=255)
                                                                           cACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCg                                                                          >  2:666979/1‑150 (MQ=255)
                                                                           cACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCg                                                                          <  2:1626955/150‑1 (MQ=255)
                                                                            aCGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCCGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCgg                                                                         <  1:2671568/150‑1 (MQ=255)
                                                                            aCGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCgg                                                                         <  1:1827196/150‑1 (MQ=255)
                                                                                aCCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTa                                                                     <  1:1688621/150‑1 (MQ=255)
                                                                                 cccGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAg                                                                    <  2:2125554/150‑1 (MQ=255)
                                                                                 cccGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAg                                                                    >  2:3211484/1‑150 (MQ=255)
                                                                                  ccGACTTCCACGGGACGTTACAGGCACATGAACTTCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCATGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGt                                                                   <  2:1334972/150‑1 (MQ=255)
                                                                                  ccGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGt                                                                   <  2:1878694/150‑1 (MQ=255)
                                                                                    gACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGc                                                                 >  1:2545405/1‑150 (MQ=255)
                                                                                    gACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGc                                                                 >  1:1999964/1‑150 (MQ=255)
                                                                                     aCTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCt                                                                >  1:1722658/1‑150 (MQ=255)
                                                                                     aCTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCt                                                                >  1:4127702/1‑150 (MQ=255)
                                                                                       ttCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTaa                                                              <  1:1964616/150‑1 (MQ=255)
                                                                                        tCCACGGGACTTTACAGGCTCATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCTATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAg                                                             <  2:1756018/150‑1 (MQ=255)
                                                                                        tCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATGTCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAg                                                             >  1:3626019/1‑150 (MQ=255)
                                                                                          cACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAggg                                                           <  1:787838/150‑1 (MQ=255)
                                                                                          cACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAggg                                                           >  1:2172251/1‑150 (MQ=255)
                                                                                            cGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGtt                                                         >  1:2598348/1‑150 (MQ=255)
                                                                                            cGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGtt                                                         <  2:3410903/150‑1 (MQ=255)
                                                                                            cGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGtt                                                         <  1:2992695/150‑1 (MQ=255)
                                                                                             gggACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGttt                                                        >  1:1150260/1‑150 (MQ=255)
                                                                                                      cAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCAc                                               >  2:1053820/1‑150 (MQ=255)
                                                                                                         gCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCt                                            <  2:821514/150‑1 (MQ=255)
                                                                                                         gCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCt                                            <  2:4017647/150‑1 (MQ=255)
                                                                                                               gAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCAtt                                      <  1:1345517/150‑1 (MQ=255)
                                                                                                               gAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCAtt                                      >  1:2157481/1‑150 (MQ=255)
                                                                                                                  cATCAGGCGCTGGTCTGGTGCTCACCTGACGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTgca                                   <  2:3094991/150‑1 (MQ=255)
                                                                                                                  cATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTgca                                   <  1:1653939/150‑1 (MQ=255)
                                                                                                                    tCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTgcagc                                 >  1:1579753/1‑150 (MQ=255)
                                                                                                                          gcTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTAtt                           <  1:517497/150‑1 (MQ=255)
                                                                                                                                    tGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCtt                 <  1:1597916/150‑1 (MQ=255)
                                                                                                                                    tGCTCACCTGAAGAGGAGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCtt                 <  1:3585465/150‑1 (MQ=255)
                                                                                                                                      cTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTa               <  1:2438463/150‑1 (MQ=255)
                                                                                                                                      cTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTa               <  2:2510897/150‑1 (MQ=255)
                                                                                                                                       tCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTaa              <  2:1271398/150‑1 (MQ=255)
                                                                                                                                         aCCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAAtt            <  1:3270243/150‑1 (MQ=255)
                                                                                                                                          ccTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTa           <  2:3922193/150‑1 (MQ=255)
                                                                                                                                            tGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTAcc         <  1:124940/150‑1 (MQ=255)
                                                                                                                                            tGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTGTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTAATCCTTCGGCGTTAATTAcc         >  1:3230591/1‑150 (MQ=255)
                                                                                                                                             gAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTAATAGAGGCGTTTATGGCTTTAAGCGCCCCCACACCCGCCGCTTAGTGCGAAGGGGTTTGTCAACACCCATGCATTGCAGCAGAATACCTTCGGCTTTAATAACCg        >  1:4042132/1‑150 (MQ=255)
                                                                                                                                                gaggagCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGccc     >  2:796380/1‑150 (MQ=255)
                                                                                                                                                gaggagCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGccc     >  2:855877/1‑150 (MQ=255)
                                                                                                                                                  ggCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCt   <  1:452264/150‑1 (MQ=255)
                                                                                                                                                   gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCtt  <  2:2176149/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTT  >  CP009273/1835924‑1836220

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: