Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,557 +C 97.7% intergenic (+29/+165) purD → / ← sdaA phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,566,5572.→C97.7% 2.4 0.0 ‑0.9129intergenic (+29/+165)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (34/92):  ref base (1/2):  total (35/94)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.19e-01

Only 100 of 151 total aligned reads displayed.

ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGC  >  NC_014034/2566458‑2566651
                                                                                                       |                                                                                                 
gtgcgagaggcGATGCGCTCGACTGGCCCGAGGGGTTTTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTg                                                                                                    <  2:1442056/91‑1 (MQ=255)
aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTg                                                                                                    <  2:1717287/101‑1 (MQ=255)
 tGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTgg                                                                                                   <  2:3120958/101‑2 (MQ=255)
    gagggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTggcct                                                                                                <  2:1539579/98‑5 (MQ=38)
     aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCGGCCCCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTggcctt                                                                                               <  1:3812637/101‑6 (MQ=38)
      tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTggccttc                                                                                              <  1:1302888/101‑7 (MQ=38)
      tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTggccttc                                                                                              <  2:2142979/101‑7 (MQ=38)
      tGGTCGATGCGAGCGATTGCCCCGAGGGCTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCGGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTggccttc                                                                                              <  2:2970871/101‑7 (MQ=25)
       ggTGGATGCGATCGACTGGCCCGGGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGGCG.AAACGGGGGGCGCTGCCCCC.CTCTTggccttcg                                                                                             <  2:1433428/101‑8 (MQ=37)
              gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAAtt                                                                                      <  1:441387/101‑1 (MQ=37)
                 aTCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.c                                                                                  <  2:1593355/101‑1 (MQ=37)
                   cGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.ccc                                                                                <  2:2502891/101‑1 (MQ=37)
                     aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGGGCGCGGTCTGAGACG.AAACGGGGGGCGCCGGCCCC.CTCTTGGCCTTGGGCCCATTCA.ccccc                                                                              >  1:3689634/1‑101 (MQ=255)
                      cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.cccccc                                                                             <  2:193483/101‑1 (MQ=25)
                       tgtcccGAGGTGTTCTGCCGCCGCGAGAGCGGCTGGCCCGCGCTCTGAGACG.AAGCGGGGGGCGCTACCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCg                                                                            <  1:784225/98‑1 (MQ=16)
                       tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCg                                                                            <  2:3257300/101‑1 (MQ=25)
                        ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCCATTCA.cccccccg                                                                           >  2:3289298/1‑99 (MQ=25)
                         gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCggg                                                                          <  2:821208/101‑1 (MQ=25)
                           ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCCATTCA.CCCCCCCGGat                                                                        >  2:3591090/1‑101 (MQ=25)
                           ccGAGGGGTTATGCCGCCGCGACCTCGGCTGGCGCGCGCTCTGAGTCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGat                                                                        <  1:3406132/101‑1 (MQ=21)
                            cgagctgTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGata                                                                       <  2:658737/95‑1 (MQ=21)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGatat                                                                      <  2:832863/101‑1 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGatat                                                                      <  2:8710/101‑1 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGatat                                                                      <  1:1022552/101‑1 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGatat                                                                      <  1:2581725/101‑1 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCCGGatat                                                                      >  1:2267035/1‑101 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGCCCAATTCA.CCCCCCGGGatat                                                                      <  1:3286874/101‑1 (MQ=25)
                                gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTc                                                                   <  1:149209/101‑1 (MQ=21)
                                gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTc                                                                   <  2:716590/101‑1 (MQ=21)
                                 ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCg                                                                  <  1:3372459/101‑1 (MQ=21)
                                   ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCggg                                                                <  2:3811/101‑1 (MQ=21)
                                   ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCggg                                                                <  1:2045935/101‑1 (MQ=21)
                                      tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCaa                                                             <  1:1757897/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                          <  1:2301070/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                         <  1:1595489/101‑1 (MQ=21)
                                          cccgcGACATCGGCTGGCGCGCCCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                         <  2:3598144/100‑1 (MQ=21)
                                           ccgcGGCGTCGGCTGGCGCGCGCTGTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCCATTCC.CCCCCCGGGATATTTCGGGCAAGATaa                                                        <  2:2853815/101‑1 (MQ=18)
                                           ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaa                                                        <  2:4182062/101‑1 (MQ=21)
                                            cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaaa                                                       <  2:1980484/101‑1 (MQ=21)
                                            cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaaa                                                       <  2:3837105/101‑1 (MQ=21)
                                            cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaaa                                                       <  1:2582465/101‑1 (MQ=21)
                                            cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAAATaaa                                                       <  1:295598/101‑1 (MQ=21)
                                             gcgACATCGGCTGGCGCGCGCTCTGAGACGAAAAG.GGGGGCGCTGCCCCC.CTCTTGGCCTTCCGCCAATTCA.CCCCCCCGGGTATTTCGGGGAAGATaaaa                                                      >  2:3461737/1‑100 (MQ=18)
                                             gcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGTCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaag                                                      <  1:1358392/101‑1 (MQ=21)
                                               gACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaagaa                                                    <  2:1061428/101‑1 (MQ=21)
                                               gACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATAATTCGGGCAAGATAaagaa                                                    >  1:2318559/1‑101 (MQ=21)
                                                aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaaaag                                                   >  1:1271813/1‑101 (MQ=21)
                                                 cATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                                  <  1:215566/101‑1 (MQ=21)
                                                    cGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                                               <  2:1414892/101‑1 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                              <  1:3569336/101‑1 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                              <  2:3464254/101‑1 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                              >  2:1299253/1‑101 (MQ=21)
                                                      gCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                             >  1:4154432/1‑101 (MQ=21)
                                                       cTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                            <  1:4081669/101‑1 (MQ=21)
                                                        tGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                           >  1:3459804/1‑101 (MQ=21)
                                                        tGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                           <  1:2981142/101‑1 (MQ=21)
                                                          gcgcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                         <  1:4349531/101‑1 (MQ=21)
                                                          gcgcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                         <  2:53519/101‑1 (MQ=21)
                                                           cgcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                        <  1:2822833/101‑1 (MQ=21)
                                                           cgcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                        <  2:3926459/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                       <  1:2359893/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                       <  1:3001361/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                      <  1:1807017/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTc                                      >  2:413218/1‑101 (MQ=21)
                                                              gcgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                                     <  2:124700/101‑1 (MQ=21)
                                                               cgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCCGGGAAATTTCGGGCAAAAAAAAAAAGGTTTTTTCTttttt                                     >  1:3164572/1‑99 (MQ=12)
                                                               cgcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                    <  1:3596700/101‑1 (MQ=21)
                                                                gcTCTGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                                   <  1:817442/101‑1 (MQ=21)
                                                                   ctGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                >  1:1274440/1‑101 (MQ=21)
                                                                   ctGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                >  2:2689914/1‑101 (MQ=21)
                                                                    tGAGACG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGAAATTTCCTGCAAGATAAAGAAGGTCTTTTTCTTTCATCTTGc                               >  1:495608/1‑101 (MQ=18)
                                                                     gagaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                              <  2:271613/101‑1 (MQ=21)
                                                                      agaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                             <  2:4443182/101‑1 (MQ=21)
                                                                       gaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                            >  1:1777787/1‑101 (MQ=21)
                                                                       gaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                            <  1:1785610/101‑1 (MQ=21)
                                                                       gaCG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                            <  2:2243119/101‑1 (MQ=21)
                                                                         cG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                          <  1:740473/101‑1 (MQ=21)
                                                                         cG.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                          >  1:3643487/1‑101 (MQ=21)
                                                                          g.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTTGGGGAAGATAAAGAAGGGCTTTTCCTTTTATCTTTCCCCaaa                         >  1:3883112/1‑101 (MQ=17)
                                                                          g.AAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGAAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                         >  2:1377123/1‑101 (MQ=21)
                                                                            aaacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCT.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                        <  1:4394270/96‑1 (MQ=21)
                                                                            aaacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                        <  1:661025/96‑1 (MQ=21)
                                                                              acggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                      <  1:1158219/98‑1 (MQ=25)
                                                                              acggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                      >  2:4193579/4‑101 (MQ=25)
                                                                                agggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.gggg             >  1:2760966/2‑101 (MQ=0)
                                                                                 gggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.ggggg            <  1:3971574/101‑2 (MQ=0)
                                                                                 gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                   <  1:2729715/101‑1 (MQ=25)
                                                                                 gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                   <  1:3095488/101‑1 (MQ=25)
                                                                                 gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                   <  1:2802997/101‑1 (MQ=25)
                                                                                  ggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.ggggga           <  2:2846502/101‑3 (MQ=2)
                                                                                   gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg.g                <  1:3108485/101‑1 (MQ=37)
                                                                                     gcgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg.ggg              <  1:550421/101‑1 (MQ=37)
                                                                                      cgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg.gggg             <  2:4078109/101‑1 (MQ=37)
                                                                                      cgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCCAAATATCCCg.gggg             >  2:3331108/1‑101 (MQ=25)
                                                                                       gcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCGTTTCATCTTGCCCGAAATATCCCG.GGGGa            <  2:1102050/101‑1 (MQ=37)
                                                                                       gcTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGa            <  1:4128436/101‑1 (MQ=37)
                                                                                        cTGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATACTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAg           >  2:3937016/1‑101 (MQ=37)
                                                                                         tGCCCCC.CTCTTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGc          >  2:1008137/1‑101 (MQ=37)
                                                                                         tGCCCCC.CTCTTGGCCTTCGGCCAACTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAAGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGc          <  2:2564462/101‑1 (MQ=37)
                                                                                           ccccc.cTCTGG..CTTC.GCC.ATTCA.CCCCCCGGGATATTTCGGGCACGATGA..AAGG.CATCCGGTTTCATCTTGCTGCAAATATCCCGAGGGGAGGCCGATTgc  >  2:2946351/1‑101 (MQ=18)
                                                                                                       |                                                                                                 
ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGC  >  NC_014034/2566458‑2566651

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG