Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,557 +C 97.2% 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.2% 5.1 ‑2.7 ‑0.2143intergenic (+29/+165)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (37/102):  ref base (3/1):  total (40/103)
Fisher's exact test strand distribution p-value = 6.64e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.71e-01

Only 100 of 168 total aligned reads displayed.

ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGG.GCGCTGCCCCC.CTC.CT..G..CTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566458‑2566644
                                                                                                           |                                                                                           
aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..g                                                                                              <  1:628652/101‑1 (MQ=255)
aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..g                                                                                              <  2:4014949/101‑1 (MQ=255)
 tGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..gg                                                                                             <  1:1286644/101‑2 (MQ=255)
  gCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggc                                                                                            <  1:1394923/101‑3 (MQ=255)
    gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggcct                                                                                          <  1:4143399/101‑5 (MQ=39)
    gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggcct                                                                                          <  1:3895620/101‑5 (MQ=39)
      tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggccttc                                                                                        <  1:966156/101‑7 (MQ=38)
        gTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggccttcgg                                                                                      <  1:2903263/101‑9 (MQ=38)
         tCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggccttcggc                                                                                     >  2:1480413/1‑92 (MQ=37)
          cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGcc                                                                                   <  2:2103579/101‑1 (MQ=37)
          cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGcc                                                                                   <  1:4189647/101‑1 (MQ=37)
           gATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCa                                                                                  <  2:373762/101‑1 (MQ=37)
            aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTAGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCaa                                                                                 <  1:3924085/101‑1 (MQ=25)
              gCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAAtt                                                                               <  2:2292218/101‑1 (MQ=37)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTc                                                                              <  1:3216613/101‑1 (MQ=37)
                gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCa                                                                             <  2:2878599/101‑1 (MQ=37)
                gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCa                                                                             <  1:2629475/101‑1 (MQ=37)
                  tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.cc                                                                          <  1:2179857/101‑1 (MQ=37)
                    gACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAAG.GGGG.GGGCTGCCCCC.CTC.TT..GGCCTTC.GGCCCATTCA.cccc                                                                        >  2:1825213/1‑101 (MQ=25)
                    gACTGGCCCGAGAGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.CGCCCATTCC.cccc                                                                        >  1:2529402/1‑101 (MQ=25)
                     aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.ccccc                                                                       <  2:898937/101‑1 (MQ=25)
                      cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.cccccc                                                                      <  1:1662172/101‑1 (MQ=25)
                       tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCg                                                                     <  1:730590/101‑1 (MQ=25)
                        ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCgg                                                                    <  2:3842005/101‑1 (MQ=25)
                        ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCgg                                                                    <  1:52255/101‑1 (MQ=25)
                        ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCgg                                                                    <  1:4585449/101‑1 (MQ=25)
                          gccGGGGGTTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.a                                                                 <  2:2989507/100‑1 (MQ=21)
                              aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATAtt                                                             <  1:2930665/101‑1 (MQ=21)
                               ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATAttt                                                            <  1:634884/101‑1 (MQ=21)
                                gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATACTTc                                                           <  1:3591918/101‑1 (MQ=21)
                                 ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGATACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCg                                                          <  1:2312179/101‑1 (MQ=21)
                                  gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCgg                                                         <  1:1572687/101‑1 (MQ=21)
                                  gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCCTC.TGCCCATTCA.CCCCCCCGG.AGATTTCgg                                                         >  2:501699/1‑101 (MQ=21)
                                   ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCggg                                                        <  2:555777/101‑1 (MQ=21)
                                   ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCggg                                                        >  1:2827380/1‑101 (MQ=21)
                                   ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCCGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCCGG.AAATTTTggg                                                        >  2:4079914/1‑101 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAAG.GGGG.GCGCTGCCCCC.CCC.TT..GGCCTTT.GGCCAATTCA.CCCCCCCGG.AAATGTCGGGc                                                       >  2:2268209/1‑101 (MQ=18)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGc                                                       <  2:2554543/101‑1 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGc                                                       <  1:1433966/101‑1 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGc                                                       <  1:3531817/101‑1 (MQ=21)
                                     cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TTTGG..CTTC.GGCCCATTCACCCCCCCGGGCAAATTTCggg                                                        >  2:3330359/1‑101 (MQ=21)
                                      tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAAG.GGGG.GGGCTGCCCCC.CTCTTT..G..CCTCTGGCCAATTCA.CCCCCCCGG.GAATTTTCGGCaa                                                     >  1:3995048/1‑101 (MQ=16)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAg                                                    <  1:1027741/101‑1 (MQ=21)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAg                                                    <  1:4538410/101‑1 (MQ=21)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAg                                                    <  2:4245141/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGAt                                                  <  2:4310610/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATa                                                 <  2:1969619/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATa                                                 <  2:4475018/101‑1 (MQ=21)
                                           ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATaa                                                <  1:1149516/101‑1 (MQ=21)
                                              cgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAaaga                                             >  1:2809323/1‑101 (MQ=21)
                                                aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAaagaag                                           <  2:3109085/101‑1 (MQ=21)
                                                    cGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCt                                       <  1:3690360/101‑1 (MQ=21)
                                                    cGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCt                                       <  1:2317065/101‑1 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  1:823683/101‑1 (MQ=21)
                                                      gCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCC.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCttt                                     <  1:3311393/101‑1 (MQ=21)
                                                      gCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCttt                                     <  1:399936/101‑1 (MQ=21)
                                                      gCTGGCGCGCGCGCTGATACG.AAAC.GGGGTGCGCTCCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.TTATTTCGGGCAAGATAAAGAAGGTCttt                                     <  2:776877/101‑1 (MQ=21)
                                                       cTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  2:269921/101‑1 (MQ=21)
                                                          gcgcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  2:2512811/101‑1 (MQ=21)
                                                          gcgcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 <  1:3271464/101‑1 (MQ=21)
                                                           cgcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCC.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:1621234/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  1:2863293/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGGTAAAGAAGGTCTTTTCCTTTc                              <  1:1909872/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  2:927318/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  1:4442676/101‑1 (MQ=21)
                                                              gcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  2:355323/101‑1 (MQ=21)
                                                              gcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  2:1107032/101‑1 (MQ=21)
                                                              gcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  1:3657419/101‑1 (MQ=21)
                                                              gatctctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  1:949411/98‑1 (MQ=21)
                                                               cgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                            >  1:2541633/1‑101 (MQ=21)
                                                                gcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  2:358972/101‑1 (MQ=21)
                                                                gcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  2:3829338/101‑1 (MQ=21)
                                                                   ctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  1:2426478/101‑1 (MQ=21)
                                                                   ctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  1:3848187/101‑1 (MQ=21)
                                                                       gaCG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                    <  1:206365/101‑1 (MQ=21)
                                                                        aCG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                   >  1:3476249/1‑101 (MQ=21)
                                                                          g.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  1:433325/101‑1 (MQ=21)
                                                                          g.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  1:4492940/101‑1 (MQ=21)
                                                                             ggcagggg.gCGCTG.CCCC.CTC.TT..GGCCTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgg         >  1:4310023/5‑101 (MQ=0)
                                                                             aacggggg.gCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata               <  2:2855141/97‑1 (MQ=25)
                                                                              gcagggg.gCGCTG.CCCC.CTC.TT..GGCCTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg        >  1:1479095/4‑101 (MQ=0)
                                                                              gcagggg.gCGCTG.CCCC.CTC.TT..GGCCTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg        >  2:1564774/4‑101 (MQ=0)
                                                                              acggggg.gCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGACTTTTTCTTTCATCTTTCCCGAAaaat              >  1:266360/4‑101 (MQ=18)
                                                                                ggggg.gCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  2:3092070/100‑1 (MQ=25)
                                                                                agggg.gCGCTGCCCCCTCTT.GG..C..CTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:3136529/2‑101 (MQ=0)
                                                                                agggg.gCGCTG.CCCC.CTC.TT..GGCCTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:3636398/2‑101 (MQ=0)
                                                                                agggg.gCGCTG.CCCC.CTC.TT..GGCCTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:4097678/2‑101 (MQ=0)
                                                                                agggg.gCGCTG.CCCC.CTC.TT..GGCCTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  2:2229699/2‑101 (MQ=0)
                                                                                agggg.gCGCTG.CCCC.CTC.TT..GGCCTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCCTGCCGCAAATATCCCggggg      >  1:4292733/2‑101 (MQ=2)
                                                                                 gggg.gCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCC.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc           <  2:1496807/101‑1 (MQ=25)
                                                                                  ggg.gCGCTGCCCCCTCTT.GG..C..CTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga    <  1:1280857/101‑3 (MQ=0)
                                                                                  ggg.gCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg          <  1:4310501/101‑1 (MQ=25)
                                                                                   gg.gCGCTGCCCCCTCTT.GG..C..CTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  2:453797/1‑98 (MQ=2)
                                                                                   gg.gCGCTG.CCCC.CTC.TT..GGCCTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  2:2606470/1‑98 (MQ=2)
                                                                                   gg.gCGCTG.CCCC.CTC.TT..GGCCTTT.GGCCAATTCA..CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  2:2471942/1‑98 (MQ=2)
                                                                                    g.gCGCTGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggg        <  1:1041361/101‑1 (MQ=37)
                                                                                          tGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTTGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTTCCCGAAATATCCCGGGGGAGc   >  2:1791967/1‑101 (MQ=25)
                                                                                          tGCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGc   >  1:1858119/1‑101 (MQ=37)
                                                                                           gCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.GGATTTTGGGCCAGATAAAGAAGGGCTTTTTCTTTTATCTTTCCCGAAATATTCCGGGGGGGcc  >  1:4652424/1‑101 (MQ=16)
                                                                                           gCCCCC.CTC.TT..GGCCTTC.GGCCAATTCA.CCCCCCGGG.AGATTTCGGGCCAGATAAAAAAGGGCTTTTTCTTTTATCTTTTCCGGAAAATCCCCGGGGGGcc  >  2:1755715/1‑101 (MQ=12)
                                                                                                           |                                                                                           
ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGG.GCGCTGCCCCC.CTC.CT..G..CTTC.GGCCAATTCA.CCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566458‑2566644

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG