Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,089,887 C→G 6.3% intergenic (+42/‑63) 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,8870CG 6.3% 221.3 / 7.2 64intergenic (+42/‑63)kdgD/gudP5‑dehydro‑4‑deoxyglucarate dehydratase/glucarate transporter
Reads supporting (aligned to +/- strand):  ref base C (21/39);  new base G (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 = 7.79e-02

ACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATG  >  AE015451/4089661‑4090122
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
acgaTTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACt                                                                                                                                                                                                                                    >  1:412738/1‑236 (MQ=255)
                                                 cGCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACaaa                                                                                                                                                                                                        <  1:842864/215‑1 (MQ=255)
                                                  gCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGc                                                                                                                                                                                   <  1:628887/235‑1 (MQ=255)
                                                     gTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGcaca                                                                                                                                                                                                          <  1:933709/209‑1 (MQ=255)
                                                         gCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGca                                                                                                                                                                                                            >  2:452399/1‑203 (MQ=255)
                                                              gTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGa                                                                                                                                                                                          <  1:73294/216‑1 (MQ=255)
                                                              gTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCg                                                                                                                                                                      <  1:464615/236‑1 (MQ=255)
                                                                   ggCCGGCGCCAGACTTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGCCCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGccc                                                                                                                                                                                                               <  1:350990/190‑1 (MQ=255)
                                                                   ggCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACaa                                                                                                                                                                                                         <  1:581703/196‑1 (MQ=255)
                                                                     ccGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGCCCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTa                                                                                                                                                                  <  1:123732/233‑1 (MQ=255)
                                                                     ccGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTa                                                                                                                                                                  >  2:650044/1‑233 (MQ=255)
                                                                      cGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCa                                                                                                                                                                                 >  2:691729/1‑217 (MQ=255)
                                                                           ccAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGGGGGCCTA‑C‑CCCACTCCGGCTa                                                                                                                                                                                                                             >  2:73294/1‑168 (MQ=255)
                                                                                  tGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGCTGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTgc                                                                                                                                                   <  1:920033/235‑1 (MQ=255)
                                                                                   cccGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGCCCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTACACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACaaaaa                                                                                                                                                                                                      <  2:977108/182‑1 (MQ=255)
                                                                                      ggTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTCCCAg                                                                                                                                                                                >  2:1274028/1‑202 (MQ=255)
                                                                                           cgacgCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGa                                                                                                                                          <  1:452402/235‑1 (MQ=255)
                                                                                              cgCCGGCCCGGGGCGGACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGccc                                                                                                                                                                                                               <  2:26572/163‑1 (MQ=255)
                                                                                                                cccGCTGACCGCCCTCCCCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGc                                                                                                                              <  1:295891/226‑1 (MQ=255)
                                                                                                                 ccGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTg                                                                                                                   <  1:962229/236‑1 (MQ=255)
                                                                                                                      gaccgaccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCAc                                                                                                              <  1:409458/236‑1 (MQ=255)
                                                                                                                       accgaccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCAc                                                                                                              <  1:944090/235‑1 (MQ=255)
                                                                                                                         cgaccTCCCCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACaa                                                                                                                                                                                                         <  1:246338/142‑1 (MQ=255)
                                                                                                                         cgaccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGa                                                                                                            <  1:650047/235‑1 (MQ=255)
                                                                                                                         cgaccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCAGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGa                                                                                                            <  1:1236955/235‑1 (MQ=255)
                                                                                                                            ccTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAAc                                                                                                                                                           <  2:554778/185‑1 (MQ=255)
                                                                                                                                 ccgccgAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATaa                                                                                                                                                                                                   >  2:1062701/1‑140 (MQ=255)
                                                                                                                                   gccgAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGa                                                                                                            <  1:946695/225‑1 (MQ=255)
                                                                                                                                    ccgAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGCGCGGGCCTA‑C‑CCCACTCCGGCTACAATTCGt                                                                                                                                                                                                                     >  2:197953/1‑119 (MQ=255)
                                                                                                                                       aaGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGAc                                                                                              >  1:4481/1‑235 (MQ=255)
                                                                                                                                       aaGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGAc                                                                                              <  1:444061/235‑1 (MQ=255)
                                                                                                                                            tACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACcgcgcg                                                                                        <  1:278408/236‑1 (MQ=255)
                                                                                                                                               gAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACcgcgcg                                                                                        >  2:511256/1‑233 (MQ=255)
                                                                                                                                                     cTCGCCGCGCTGATGGACCAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCACGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACGAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGAcc                                                                                                                <  1:229255/203‑1 (MQ=255)
                                                                                                                                                                gATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCa                                                                                    >  2:499159/1‑220 (MQ=255)
                                                                                                                                                                  tGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCa                                                                  <  1:691732/236‑1 (MQ=255)
                                                                                                                                                                      cAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACg                                                                                                                                                          <  2:1169770/144‑1 (MQ=255)
                                                                                                                                                                       aaGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCCGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGt                                                                                                                    <  1:296005/181‑1 (MQ=255)
                                                                                                                                                                          aTGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCa                                                           <  1:840045/235‑1 (MQ=255)
                                                                                                                                                                              gCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGaa                                                        >  2:1241632/1‑234 (MQ=255)
                                                                                                                                                                               cccGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCa                                                                                    <  1:28834/205‑1 (MQ=255)
                                                                                                                                                                                ccGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGAc                                                    <  1:80105/236‑1 (MQ=255)
                                                                                                                                                                                    aGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGGGAGCCCGCACAAAAATAACTAGCGGGa                                                                                                                                                                                          >  2:296004/1‑98 (MQ=255)
                                                                                                                                                                                     gTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGAcc                                                                                             <  1:565916/190‑1 (MQ=255)
                                                                                                                                                                                      tAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGG‑GTAGG‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAAc                                                                                                                                                           <  2:1104648/127‑1 (MQ=255)
                                                                                                                                                                                      tAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGG‑GTA‑GGCCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAAc                                                                                                                                                           >  1:1104651/1‑127 (MQ=255)
                                                                                                                                                                                        aaCGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAgg                                                                       <  1:20244/209‑1 (MQ=255)
                                                                                                                                                                                                ccGCGGCCTGCCAGACGCCGCGCCCGGAGTGGG‑GTAGG‑CCCACTCCGGCTACAATTCGTGAGCCCGCACaa                                                                                                                                                                                                         >  2:1189239/1‑71 (MQ=25)
                                                                                                                                                                                                ccGCGGCCTGCCAGACGCCGCGCCCGGAGTGGG‑GTAGG‑CCCACTCCGGCTACAATTCGTGAGCCCGCACaa                                                                                                                                                                                                         <  1:1189242/71‑1 (MQ=255)
                                                                                                                                                                                                  gcgGCCTGCCAGACGCCGGGCCCGGAGTGGGCCTA‑C‑CCCACTCCgg                                                                                                                                                                                                                                <  1:235600/46‑1 (MQ=255)
                                                                                                                                                                                                  gcgGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACCAAAATAACTAGCGGGa                                                                                                                                                                                          >  2:161194/1‑84 (MQ=255)
                                                                                                                                                                                                    ggCCTGCCAGACGCCGCGCACGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATGGCCAGCATGCAAGCGACTAAGAAAACGCAtgt                                                                                                                                                     <  2:603439/119‑1 (MQ=255)
                                                                                                                                                                                                                 ccgcGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTtata                     >  2:753497/1‑234 (MQ=255)
                                                                                                                                                                                                                     gcCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATg                                                                                                                               <  1:344042/124‑1 (MQ=255)
                                                                                                                                                                                                                     gcCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATg                                                                                                                               >  2:344041/1‑124 (MQ=255)
                                                                                                                                                                                                                     agccGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCAAGCCTGCAGAAAGACCTCGGCCTTGACGCCGTAACCCTCGGTta                       >  2:1186082/3‑228 (MQ=255)
                                                                                                                                                                                                                       ccGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAgg                                                                       >  2:796052/1‑178 (MQ=255)
                                                                                                                                                                                                                         ggAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGc           >  1:1040867/1‑236 (MQ=255)
                                                                                                                                                                                                                             tGGGCCTA‑C‑CCCCCTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGc                                                            <  2:276242/183‑1 (MQ=255)
                                                                                                                                                                                                                             tGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGTCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACTCTCgtt                         >  2:1088614/1‑216 (MQ=255)
                                                                                                                                                                                                                             tGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCtt        >  1:750364/1‑235 (MQ=255)
                                                                                                                                                                                                                             tGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTc       <  1:1241635/236‑1 (MQ=255)
                                                                                                                                                                                                                               ggCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCg      <  1:125970/235‑1 (MQ=255)
                                                                                                                                                                                                                                gCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTc       <  2:1040864/233‑1 (MQ=255)
                                                                                                                                                                                                                                  cTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATg  >  1:161462/1‑236 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
ACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCGGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATGGACAAGATGGGCCCGCAGTAAACGCTGCCCGCGGCCTGCCAGACGCCGCGCCCGGAGTGGGCCTA‑C‑CCCACTCCGGCTACAATTCGTGAGCCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCATGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATG  >  AE015451/4089661‑4090122

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: ‑