Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,572 Δ1 bp 2.5% intergenic (+44/+150) 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,5720C→. 2.5% 4.4 ‑0.5 ‑0.7159intergenic (+44/+150)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (49/106):  total (51/108)
Fisher's exact test strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.08e-01

Only 100 of 178 total aligned reads displayed.

ATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTC.CT..G..CTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCG  >  NC_014034/2566475‑2566671
                                                                                                         |                                                                                                    
aTCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.c                                                                                                      <  2:1593355/101‑1 (MQ=37)
    aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.ccccc                                                                                                  <  2:1832720/101‑1 (MQ=25)
     cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.cccccc                                                                                                 <  2:193483/101‑1 (MQ=25)
     cTGGCCAGAGGGGTTAAGCCGACACAACATCGGCTGGCGCGAGCTCTGAGAAGAAACGGGGGGCGCTGCCCCC.ATC.TT..GGCCTTCGGCCAATTCA.cccccc                                                                                                 <  2:4392946/101‑1 (MQ=12)
      tgtcccGAGGTGTTCTGCCGCCGCGAGAGCGGCTGGCCCGCGCTCTGAGACGAAGCGGGGGGCGCTACCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCg                                                                                                <  1:784225/98‑1 (MQ=16)
       ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCCATTCA.cccccccg                                                                                               >  2:3289298/1‑99 (MQ=25)
        gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGGCGCCCCC.CTC.TTGGG..CTTCGGGCAATTCA.CCCCCCcgg                                                                                              >  2:2311527/1‑101 (MQ=255)
          ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCCATTCA.CCCCCCCGGat                                                                                            >  2:3591090/1‑101 (MQ=25)
          ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGat                                                                                            >  2:747689/1‑101 (MQ=25)
          ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGat                                                                                            <  2:597069/101‑1 (MQ=25)
          ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGat                                                                                            <  1:1749655/101‑1 (MQ=25)
           cgagctgTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGata                                                                                           <  2:658737/95‑1 (MQ=21)
            gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGata.t                                                                                         <  2:8710/101‑1 (MQ=25)
            gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGata.t                                                                                         <  2:832863/101‑1 (MQ=25)
            gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGata.t                                                                                         <  1:2581725/101‑1 (MQ=25)
            gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGata.t                                                                                         <  1:1022552/101‑1 (MQ=25)
            gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCCGGata.t                                                                                         >  1:2267035/1‑101 (MQ=25)
            gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGCCCAATTCA.CCCCCCGGGata.t                                                                                         <  1:3286874/101‑1 (MQ=25)
               gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTc                                                                                      <  1:149209/101‑1 (MQ=21)
               gggTTCGGCCGCCGCGAGATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTc                                                                                      <  1:123515/101‑1 (MQ=18)
               gggTGCGGCCGCCGCGACATCGGCGGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTc                                                                                      <  1:1301414/101‑1 (MQ=18)
                  ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCggg                                                                                   <  2:3811/101‑1 (MQ=21)
                     tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCCCTTCCCCC.CTCTTT..GC.CTTCGGCCAATTCA.CCCCCCGCGGTA.TTTCGGGGaa                                                                                >  2:1266654/1‑101 (MQ=21)
                      gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTCTTT..GC.CTTCGGCCAATTCA.CCCCCCCGGATA.TTTCTGGGAAg                                                                               >  2:3968143/1‑101 (MQ=21)
                         gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATa                                                                            <  1:1595489/101‑1 (MQ=21)
                         gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATa                                                                            <  1:3625649/101‑1 (MQ=21)
                          ccgcGGCGTCGGCTGGCGCGCGCTGTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCCATTCC.CCCCCCGGGATA.TTTCGGGCAAGATaa                                                                           <  2:2853815/101‑1 (MQ=18)
                          ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATaa                                                                           <  2:4182062/101‑1 (MQ=21)
                          ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATaa                                                                           <  2:2601768/101‑1 (MQ=21)
                           cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATaaa                                                                          <  1:2582465/101‑1 (MQ=21)
                           cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATaaa                                                                          <  2:3837105/101‑1 (MQ=21)
                            tcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TCTCGGGCAAGATAaag                                                                         <  1:3750265/100‑1 (MQ=21)
                            gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGTCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAaag                                                                         <  1:1358392/101‑1 (MQ=21)
                              gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAaagaa                                                                       <  2:1061428/101‑1 (MQ=21)
                               aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCCGGAAA.ATTCGGGCAAGATAaagaag                                                                      >  1:2317880/1‑101 (MQ=21)
                                cATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAgg                                                                     <  1:215566/101‑1 (MQ=21)
                                cATCGGCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAgg                                                                     <  1:2985799/101‑1 (MQ=21)
                                    ggCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCtt                                                                 <  1:3569336/101‑1 (MQ=21)
                                    ggCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCtt                                                                 <  2:3464254/101‑1 (MQ=21)
                                    ggCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCtt                                                                 <  2:1394035/101‑1 (MQ=21)
                                    gacTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCtt                                                                 >  1:3966902/3‑101 (MQ=21)
                                     gCTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCttt                                                                >  1:4154432/1‑101 (MQ=21)
                                      cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCtttt                                                               <  2:2134379/101‑1 (MQ=21)
                                      cTGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCtttt                                                               <  1:885005/101‑1 (MQ=21)
                                       tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTc                                                              >  1:3459804/1‑101 (MQ=21)
                                       tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTc                                                              >  2:3361660/1‑101 (MQ=21)
                                       tGGCGCGCGCTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCCGGGTATTTTCGGGC.AGATAaaaaaagggttttt                                                              >  2:3027699/1‑92 (MQ=9)
                                         gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                            <  1:4349531/101‑1 (MQ=21)
                                         gcgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                            <  2:53519/101‑1 (MQ=21)
                                          cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCC.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                           <  1:2822833/101‑1 (MQ=21)
                                          cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                           <  2:3926459/101‑1 (MQ=21)
                                          cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                           <  2:1094175/101‑1 (MQ=21)
                                           gcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                          <  1:2359893/101‑1 (MQ=21)
                                            cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                         <  1:1807017/101‑1 (MQ=21)
                                            cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTc                                                         >  2:413218/1‑101 (MQ=21)
                                              cgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCCGGGAAA.TTTCGGGCAAAAAAAAAAAGGTTTTTTCTttttt                                                        >  1:3164572/1‑99 (MQ=12)
                                              cgcTCTGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                                       <  1:3596700/101‑1 (MQ=21)
                                                  ctGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                                   >  1:1274440/1‑101 (MQ=21)
                                                   tGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                                                  <  2:1574485/101‑1 (MQ=21)
                                                   tGAGACGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGAAA.TTTCCTGCAAGATAAAGAAGGTCTTTTTCTTTCATCTTGc                                                  >  1:495608/1‑101 (MQ=18)
                                                    gagaCGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                 <  2:271613/101‑1 (MQ=21)
                                                      gaCGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                               >  1:1777787/1‑101 (MQ=21)
                                                        cGAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                             <  1:740473/101‑1 (MQ=21)
                                                         gAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCCATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                            <  1:3107135/101‑1 (MQ=21)
                                                         gAAACGGGGGGCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                            <  2:4135017/101‑1 (MQ=21)
                                                          aaacggggggCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCT.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                                           <  1:4394270/96‑1 (MQ=21)
                                                           aacggggggCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                          >  1:2794101/5‑101 (MQ=25)
                                                           aacggggggCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                          <  1:3403555/97‑1 (MQ=25)
                                                               gggggCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                      <  1:2729715/101‑1 (MQ=25)
                                                               gggggCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                      <  1:3095488/101‑1 (MQ=25)
                                                               gggggCGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATGCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                      >  1:4105956/1‑101 (MQ=25)
                                                                ggggCGCTGCCCCCTCTT.GG..C..CTTTGGCCAATTCA..CCCCTGGGATA.TTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga                               <  2:2846502/101‑3 (MQ=2)
                                                                   gcgcTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgggg                                  <  1:550421/101‑1 (MQ=37)
                                                                     gcTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                <  1:2932405/101‑1 (MQ=37)
                                                                       tGCCCCC.CTC.TT..GGCCTTCGGCCAATTCA.CCCCCCGGGAAA.TTTCGGGGCAGATAAAGAGGGTTTTTTCCTTTCTTCTTTCCCGGAATATCCCCCGGGacc                              >  1:192838/1‑99 (MQ=16)
                                                                       tGCCCCC.CTC.TT..GGCCTTCGGCCAACTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAAGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGc                              <  2:2564462/101‑1 (MQ=37)
                                                                        gCCCCC.CTC.TT..GGCCTTCGGCCAATTCC.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGcc                             <  1:693658/101‑1 (MQ=38)
                                                                                ccctcttggccTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                         >  1:3981012/11‑101 (MQ=255)
                                                                                 cctcttggccTTCGTCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCtt                        <  2:3267489/92‑1 (MQ=37)
                                                                                  ctcttggccTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                       <  1:281610/93‑1 (MQ=255)
                                                                                    cttggccTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                     <  1:3958514/95‑1 (MQ=255)
                                                                                      tggccTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcggg                   <  1:747769/97‑1 (MQ=255)
                                                                                      tggccTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcggg                   >  1:2972992/5‑101 (MQ=255)
                                                                                          cTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGcgggcggg               >  2:2325259/1‑101 (MQ=255)
                                                                                             cGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc            <  1:329854/101‑1 (MQ=255)
                                                                                             cGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc            >  2:1158868/1‑101 (MQ=255)
                                                                                             cGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc            >  2:1456184/1‑101 (MQ=255)
                                                                                              ggCCAATTCG.CCCCCCGGGATA.TTTCGGGCGAGATAAAGAAGGTCTGTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa           <  2:3689381/101‑1 (MQ=255)
                                                                                                ccAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGGggggggggc         >  1:1153559/1‑97 (MQ=255)
                                                                                                 cAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcg        <  1:3361897/101‑1 (MQ=255)
                                                                                                    ttCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccc     >  1:3644170/1‑101 (MQ=255)
                                                                                                     tCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc    >  2:3359653/1‑101 (MQ=255)
                                                                                                     tCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc    <  2:2266890/101‑1 (MQ=255)
                                                                                                     tCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc    >  1:1776744/1‑101 (MQ=255)
                                                                                                      cA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTTCCCGAAATATCCCGGGGGACCCCGCCTGCGGGCGGGGGCCGGGccccc   >  1:209796/1‑101 (MQ=37)
                                                                                                      cA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc   >  2:1404180/1‑101 (MQ=255)
                                                                                                      cA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc   <  2:2392515/101‑1 (MQ=255)
                                                                                                      cA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGGCTTTTTCTTTCATCTTGCCCGAAATATCCCGGGGGGGCCCGCTTGCGGGGGGGGGGCGCGccccc   >  2:3072127/1‑101 (MQ=21)
                                                                                                        ccccctcCGGATA.TTTCGGGCAGGATAAAGCAGGTCTTTTCCCTCCACGTTGTCGCAAACACCCCGGGGGGGCCCGCCAGCGGGCGGGGGCAGCGCCCCCg  <  2:1427563/100‑1 (MQ=255)
                                                                                                        cccccccGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg  <  2:1374118/100‑1 (MQ=255)
                                                                                                         |                                                                                                    
ATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCC.CTC.CT..G..CTTCGGCCAATTCA.CCCCCCGGGATA.TTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCG  >  NC_014034/2566475‑2566671

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