Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,089,886 Δ1 bp 6.2% intergenic (+41/‑64) kdgD → / → gudP 5‑dehydro‑4‑deoxyglucarate dehydratase/glucarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,089,8860C. 6.2% 222.3 / 7.9 64intergenic (+41/‑64)kdgD/gudP5‑dehydro‑4‑deoxyglucarate dehydratase/glucarate transporter
Reads supporting (aligned to +/- strand):  ref base C (21/39);  new base . (2/2);  total (23/41)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-02

ACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCG  >  AE015451/4089661‑4090118
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          
acgaTTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACt                                                                                                                                                                                                                                >  1:412738/1‑236 (MQ=255)
                                                 cGCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACaaa                                                                                                                                                                                                    <  1:842864/215‑1 (MQ=255)
                                                  gCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGc                                                                                                                                                                               <  1:628887/235‑1 (MQ=255)
                                                     gTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGcaca                                                                                                                                                                                                      <  1:933709/209‑1 (MQ=255)
                                                         gCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGca                                                                                                                                                                                                        >  2:452399/1‑203 (MQ=255)
                                                              gTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGa                                                                                                                                                                                      <  1:73294/216‑1 (MQ=255)
                                                              gTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCg                                                                                                                                                                  <  1:464615/236‑1 (MQ=255)
                                                                   ggCCGGCGCCAGACTTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGCCCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGccc                                                                                                                                                                                                           <  1:350990/190‑1 (MQ=255)
                                                                   ggCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACaa                                                                                                                                                                                                     <  1:581703/196‑1 (MQ=255)
                                                                     ccGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGCCCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTa                                                                                                                                                              <  1:123732/233‑1 (MQ=255)
                                                                     ccGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTa                                                                                                                                                              >  2:650044/1‑233 (MQ=255)
                                                                      cGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCa                                                                                                                                                                             >  2:691729/1‑217 (MQ=255)
                                                                           ccAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGGGGGCCTACCCCACTCCGGCTa                                                                                                                                                                                                                         >  2:73294/1‑168 (MQ=255)
                                                                                  tGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGCTGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTgc                                                                                                                                               <  1:920033/235‑1 (MQ=255)
                                                                                   cccGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGCCCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTACACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACaaaaa                                                                                                                                                                                                  <  2:977108/182‑1 (MQ=255)
                                                                                      ggTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTCCCAg                                                                                                                                                                            >  2:1274028/1‑202 (MQ=255)
                                                                                           cgacgCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGa                                                                                                                                      <  1:452402/235‑1 (MQ=255)
                                                                                              cgCCGGCCCGGGGCGGACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGccc                                                                                                                                                                                                           <  2:26572/163‑1 (MQ=255)
                                                                                                                cccGCTGACCGCCCTCCCCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGc                                                                                                                          <  1:295891/226‑1 (MQ=255)
                                                                                                                 ccGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTg                                                                                                               <  1:962229/236‑1 (MQ=255)
                                                                                                                      gaccgaccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCAc                                                                                                          <  1:409458/236‑1 (MQ=255)
                                                                                                                       accgaccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCAc                                                                                                          <  1:944090/235‑1 (MQ=255)
                                                                                                                         cgaccTCCCCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACaa                                                                                                                                                                                                     <  1:246338/142‑1 (MQ=255)
                                                                                                                         cgaccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGa                                                                                                        <  1:650047/235‑1 (MQ=255)
                                                                                                                         cgaccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCAGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGa                                                                                                        <  1:1236955/235‑1 (MQ=255)
                                                                                                                            ccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAAc                                                                                                                                                       <  2:554778/185‑1 (MQ=255)
                                                                                                                                 ccgccgAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATaa                                                                                                                                                                                               >  2:1062701/1‑140 (MQ=255)
                                                                                                                                   gccgAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGa                                                                                                        <  1:946695/225‑1 (MQ=255)
                                                                                                                                    ccgAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGCGCGGGCCTACCCCACTCCGGCTACAATTCGt                                                                                                                                                                                                                 >  2:197953/1‑119 (MQ=255)
                                                                                                                                       aaGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGAc                                                                                          >  1:4481/1‑235 (MQ=255)
                                                                                                                                       aaGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGAc                                                                                          <  1:444061/235‑1 (MQ=255)
                                                                                                                                            tACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACcgcgcg                                                                                    <  1:278408/236‑1 (MQ=255)
                                                                                                                                               gAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACcgcgcg                                                                                    >  2:511256/1‑233 (MQ=255)
                                                                                                                                                     cTCGCCGCGCTGATGGACCAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCACGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACGAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGAcc                                                                                                            <  1:229255/203‑1 (MQ=255)
                                                                                                                                                                gATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCa                                                                                >  2:499159/1‑220 (MQ=255)
                                                                                                                                                                  tGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCa                                                              <  1:691732/236‑1 (MQ=255)
                                                                                                                                                                      cAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACg                                                                                                                                                      <  2:1169770/144‑1 (MQ=255)
                                                                                                                                                                       aaGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCCGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGt                                                                                                                <  1:296005/181‑1 (MQ=255)
                                                                                                                                                                          aTGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCa                                                       <  1:840045/235‑1 (MQ=255)
                                                                                                                                                                              gCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGaa                                                    >  2:1241632/1‑234 (MQ=255)
                                                                                                                                                                               cccGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCa                                                                                <  1:28834/205‑1 (MQ=255)
                                                                                                                                                                                ccGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGAc                                                <  1:80105/236‑1 (MQ=255)
                                                                                                                                                                                    aGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGGGAGCCCGCACAAAAATAACTAGCGGGa                                                                                                                                                                                      >  2:296004/1‑98 (MQ=255)
                                                                                                                                                                                     gTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGAcc                                                                                         <  1:565916/190‑1 (MQ=255)
                                                                                                                                                                                      tAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGGTAGGCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAAc                                                                                                                                                       <  2:1104648/127‑1 (MQ=255)
                                                                                                                                                                                      tAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGGTAGGCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAAc                                                                                                                                                       >  1:1104651/1‑127 (MQ=255)
                                                                                                                                                                                        aaCGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAgg                                                                   <  1:20244/209‑1 (MQ=255)
                                                                                                                                                                                                ccGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGGTAGGCCCACTCCGGCTACAATTCGTGAGCCCGCACaa                                                                                                                                                                                                     >  2:1189239/1‑71 (MQ=25)
                                                                                                                                                                                                ccGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGGTAGGCCCACTCCGGCTACAATTCGTGAGCCCGCACaa                                                                                                                                                                                                     <  1:1189242/71‑1 (MQ=255)
                                                                                                                                                                                                  gcgGCCTGCCAGACGCCGGGCCCGGAGTGGGCCTACCCCACTCCgg                                                                                                                                                                                                                            <  1:235600/46‑1 (MQ=255)
                                                                                                                                                                                                  gcgGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACCAAAATAACTAGCGGGa                                                                                                                                                                                      >  2:161194/1‑84 (MQ=255)
                                                                                                                                                                                                    ggCCTGCCAGACGCCGCGCACGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATGGCCAGCATGCAAGCGACTAAGAAAACGCAtgt                                                                                                                                                 <  2:603439/119‑1 (MQ=255)
                                                                                                                                                                                                                 ccgcGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTtata                 >  2:753497/1‑234 (MQ=255)
                                                                                                                                                                                                                     gcCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATg                                                                                                                           >  2:344041/1‑124 (MQ=255)
                                                                                                                                                                                                                     gcCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATg                                                                                                                           <  1:344042/124‑1 (MQ=255)
                                                                                                                                                                                                                     agccGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCAAGCCTGCAGAAAGACCTCGGCCTTGACGCCGTAACCCTCGGTta                   >  2:1186082/3‑228 (MQ=255)
                                                                                                                                                                                                                       ccGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAgg                                                                   >  2:796052/1‑178 (MQ=255)
                                                                                                                                                                                                                         ggAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGc       >  1:1040867/1‑236 (MQ=255)
                                                                                                                                                                                                                             tGGGCCTACCCCCCTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGc                                                        <  2:276242/183‑1 (MQ=255)
                                                                                                                                                                                                                             tGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGTCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACTCTCgtt                     >  2:1088614/1‑216 (MQ=255)
                                                                                                                                                                                                                             tGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCtt    >  1:750364/1‑235 (MQ=255)
                                                                                                                                                                                                                             tGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTc   <  1:1241635/236‑1 (MQ=255)
                                                                                                                                                                                                                               ggCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCg  <  1:125970/235‑1 (MQ=255)
                                                                                                                                                                                                                                gCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTc   <  2:1040864/233‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          
ACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTACCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCG  >  AE015451/4089661‑4090118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: