Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,561 C→T 2.8% intergenic (+33/+161) 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,5610C→T 2.8% 2.9 ‑1.1 ‑0.8145intergenic (+33/+161)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (39/102):  total (41/104)
Fisher's exact test strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.77e-01

Only 100 of 160 total aligned reads displayed.

GCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGGCGCT.GCCCCC.CTC.CTG..CTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGC  >  NC_014034/2566472‑2566661
                                                                                                |                                                                                                      
   aTCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.c                                                                                            <  2:1593355/101‑1 (MQ=37)
     cGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.ccc                                                                                          <  2:2502891/101‑1 (MQ=37)
       aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGGGCGCGGTCTGAGACG.AAACGGGGGGCGCC.GGCCCC.CTC.TTGGCCTTGGGCCCATTCA.ccccc                                                                                        >  1:3689634/1‑101 (MQ=255)
       aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.ccccc                                                                                        <  2:1832720/101‑1 (MQ=25)
        cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.cccccc                                                                                       <  2:193483/101‑1 (MQ=25)
        cTGGCCAGAGGGGTTAAGCCGACACAACATCGGCTGGCGCGAGCTCTGAGAAG.AAACGGGGGGCGCT.GCCCCC.ATC.TTGGCCTTCGGCCAATTCA.cccccc                                                                                       <  2:4392946/101‑1 (MQ=12)
         tgtcccGAGGTGTTCTGCCGCCGCGAGAGCGGCTGGCCCGCGCTCTGAGACG.AAGCGGGGGGCGCT.ACCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCg                                                                                      <  1:784225/98‑1 (MQ=16)
         tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCg                                                                                      <  2:3257300/101‑1 (MQ=25)
          ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCCATTCA.cccccccg                                                                                     >  2:3289298/1‑99 (MQ=25)
           gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCggg                                                                                    <  2:821208/101‑1 (MQ=25)
             ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCCATTCA.CCCCCCCGGat                                                                                  >  2:3591090/1‑101 (MQ=25)
             ccGAGGGGTTATGCCGCCGCGACCTCGGCTGGCGCGCGCTCTGAGTCG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGat                                                                                  <  1:3406132/101‑1 (MQ=21)
              cgagctgTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGata                                                                                 <  2:658737/95‑1 (MQ=21)
               gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGGGGTGGCCCCC.CCT.TTGC.CTTCGGCCCATTCA.CCCCCCcggcaat                                                                                >  2:3764620/1‑97 (MQ=255)
               gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGatat                                                                                <  2:8710/101‑1 (MQ=25)
               gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGatat                                                                                <  2:832863/101‑1 (MQ=25)
               gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGatat                                                                                <  1:1022552/101‑1 (MQ=25)
               gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGatat                                                                                <  1:2581725/101‑1 (MQ=25)
               gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCCGGatat                                                                                >  1:2267035/1‑101 (MQ=25)
               gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGCCCAATTCA.CCCCCCGGGatat                                                                                <  1:3286874/101‑1 (MQ=25)
                  gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTc                                                                             <  2:716590/101‑1 (MQ=21)
                  gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTc                                                                             <  1:149209/101‑1 (MQ=21)
                   ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCg                                                                            <  1:3372459/101‑1 (MQ=21)
                     ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCggg                                                                          <  1:2045935/101‑1 (MQ=21)
                        tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCaa                                                                       <  1:1757897/101‑1 (MQ=21)
                         gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTCTTTGC.CTTCGGCCAATTCA.CCCCCCCGGATATTTCTGGGAAg                                                                      >  2:3968143/1‑101 (MQ=21)
                         gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAg                                                                      <  2:3544526/101‑1 (MQ=21)
                          ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGa                                                                     <  1:4402108/101‑1 (MQ=21)
                          ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGa                                                                     <  1:768504/101‑1 (MQ=21)
                           cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGAt                                                                    <  1:2301070/101‑1 (MQ=21)
                            gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                                   <  1:1595489/101‑1 (MQ=21)
                            gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATa                                                                   <  1:3625649/101‑1 (MQ=21)
                             ccgcGGCGTCGGCTGGCGCGCGCTGTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCCATTCC.CCCCCCGGGATATTTCGGGCAAGATaa                                                                  <  2:2853815/101‑1 (MQ=18)
                             ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaa                                                                  <  2:4182062/101‑1 (MQ=21)
                              cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATaaa                                                                 <  2:3837105/101‑1 (MQ=21)
                              cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAAATaaa                                                                 <  1:295598/101‑1 (MQ=21)
                               tcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATCTCGGGCAAGATAaag                                                                <  1:3750265/100‑1 (MQ=21)
                               gcgACATCGGCTGGCGCGCGCTCTGAGACGAAAAG.GGGGGCGCT.GCCCCC.CTC.TTGGCCTTCCGCCAATTCA.CCCCCCCGGGTATTTCGGGGAAGATaaaa                                                                >  2:3461737/1‑100 (MQ=18)
                               gcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGTCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaag                                                                <  1:1358392/101‑1 (MQ=21)
                                 gACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaagaa                                                              <  2:3920021/101‑1 (MQ=21)
                                 gACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaagaa                                                              <  2:1061428/101‑1 (MQ=21)
                                 gACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATAATTCGGGCAAGATAaagaa                                                              >  1:2318559/1‑101 (MQ=21)
                                  aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAaaaaag                                                             >  1:1271813/1‑101 (MQ=21)
                                   cATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                                            <  1:215566/101‑1 (MQ=21)
                                      cGGCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                                                         <  2:1414892/101‑1 (MQ=21)
                                       ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                        <  1:3569336/101‑1 (MQ=21)
                                       ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                        <  2:1394035/101‑1 (MQ=21)
                                       ggCTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                                        >  2:1299253/1‑101 (MQ=21)
                                         cTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                      <  2:2134379/101‑1 (MQ=21)
                                         cTGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                                      <  1:885005/101‑1 (MQ=21)
                                          tGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCGCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                     <  1:2647695/101‑1 (MQ=21)
                                          tGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGGTAATTTGGGCAAGATAAAAAAggggctttt                                                     >  1:2392665/1‑94 (MQ=12)
                                          tGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                     <  1:2981142/101‑1 (MQ=21)
                                          tGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                     >  2:3361660/1‑101 (MQ=21)
                                          tGGCGCGCGCTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                                     >  1:3459804/1‑101 (MQ=21)
                                             cgcgcgcTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                  <  1:2822833/101‑1 (MQ=21)
                                             cgcgcgcTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                  <  2:1094175/101‑1 (MQ=21)
                                              gcgcgcTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                 <  1:3001361/101‑1 (MQ=21)
                                               cgcgcTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                <  1:1807017/101‑1 (MQ=21)
                                               cgcgcTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTc                                                >  2:413218/1‑101 (MQ=21)
                                                gcgcTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                                               <  2:124700/101‑1 (MQ=21)
                                                 cgcTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCCGGGAAATTTCGGGCAAAAAAAAAAAGGTTTTTTCTttttt                                               >  1:3164572/1‑99 (MQ=12)
                                                 cgcTCTGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                              <  1:3596700/101‑1 (MQ=21)
                                                     ctGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                          >  2:2689914/1‑101 (MQ=21)
                                                     ctGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                          >  1:1274440/1‑101 (MQ=21)
                                                     ctGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCAAATTCG.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                                          <  1:4265530/101‑1 (MQ=21)
                                                      tGAGACG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGAAATTTCCTGCAAGATAAAGAAGGTCTTTTTCTTTCATCTTGc                                         >  1:495608/1‑101 (MQ=18)
                                                       gagaCG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                        <  2:271613/101‑1 (MQ=21)
                                                       gagaCG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTGCCTTTCATCTTGcc                                        <  1:519025/101‑1 (MQ=21)
                                                        agaCG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                       <  2:4443182/101‑1 (MQ=21)
                                                         gaCG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                      >  1:1777787/1‑101 (MQ=21)
                                                         gaCG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                      <  2:2243119/101‑1 (MQ=21)
                                                           cG.AAACGGGGGGCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                    <  1:740473/101‑1 (MQ=21)
                                                              aaacggggggCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCT.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                                  <  1:4394270/96‑1 (MQ=21)
                                                              aaacggggggCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                                  <  1:661025/96‑1 (MQ=21)
                                                               aacggggggCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                 >  1:2794101/5‑101 (MQ=25)
                                                                acggggggCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                <  1:1158219/98‑1 (MQ=25)
                                                                   gggggCGCT.GCCCCCTCTT.GGC..CTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.ggggg                      <  1:3971574/101‑2 (MQ=0)
                                                                   gggggCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                             <  1:2729715/101‑1 (MQ=25)
                                                                   gggggCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                             <  1:3095488/101‑1 (MQ=25)
                                                                   gggggCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATGCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                             >  1:4105956/1‑101 (MQ=25)
                                                                    ggggCGCT.GCCCCCTCTT.GGC..CTTTGGCCAATTCA..CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCg.ggggga                     <  2:2846502/101‑3 (MQ=2)
                                                                     gggCGCT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg.g                          <  1:3108485/101‑1 (MQ=37)
                                                                         gcT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGa                      <  1:2932405/101‑1 (MQ=37)
                                                                         gcT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGa                      <  1:4128436/101‑1 (MQ=37)
                                                                          cT.GCCCCC.CTC.TTGGCCTTCGGCCAATTCA.CCCCCCGGGATACTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAg                     >  2:3937016/1‑101 (MQ=37)
                                                                           t.GCCCCC.CTC.TTGGCCTTCGGCCAACTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAAGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGc                    <  2:2564462/101‑1 (MQ=37)
                                                                              ccccc.cTC.TGG..CTTC.GCC.ATTCA.CCCCCCGGGATATTTCGGGCACGATGA..AAGG.CATCCGGTTTCATCTTGCTGCAAATATCCCGAGGGGAGGCCGATTgc            >  2:2946351/1‑101 (MQ=18)
                                                                                   ccctcttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCt               >  1:3981012/11‑101 (MQ=255)
                                                                                    cctcttggccTTCGTCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCtt              <  2:3267489/92‑1 (MQ=37)
                                                                                    cctcttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGCAATATCCCG.GGGGAGCCCGCtt              >  2:4177887/10‑101 (MQ=37)
                                                                                     ctcttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTg             <  1:281610/93‑1 (MQ=255)
                                                                                      tcttggcctccGGCCAATTCC.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgc            <  2:1499170/91‑1 (MQ=38)
                                                                                       cttggccTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTgcg           <  1:3958514/95‑1 (MQ=255)
                                                                                             cTTCGGGCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGcgggcggg     <  1:725245/101‑1 (MQ=255)
                                                                                             cTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGcgggcggg     >  2:2325259/1‑101 (MQ=255)
                                                                                              ctCGGCCAATTCA.CCCTCCCGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCgggg    <  2:563537/100‑1 (MQ=255)
                                                                                                cGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGc  >  2:1158868/1‑101 (MQ=255)
                                                                                                cGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGc  >  2:1456184/1‑101 (MQ=255)
                                                                                                cGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGc  <  1:329854/101‑1 (MQ=255)
                                                                                                |                                                                                                      
GCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGGCGCT.GCCCCC.CTC.CTG..CTTCGGCCAATTCA.CCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCG.GGGGAGCCCGCTTGCGGGCGGGGGC  >  NC_014034/2566472‑2566661

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG