Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,555 C→T 100% intergenic (+27/+167) 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,5550C→T100.0% 789.2 278intergenic (+27/+167)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (75/201):  ref base (0/0):  total (76/202)

Only 100 of 300 total aligned reads displayed.

CTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GG.GG.GCGCTGCCCCCCTC.CT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566456‑2566644
                                                                                                       |                                                                                             
  tggCGGTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGGGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..g                                                                                           <  1:3691558/99‑1 (MQ=255)
  aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..g                                                                                           <  2:2247535/101‑1 (MQ=255)
     cGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..ggcc                                                                                        <  1:1414556/101‑4 (MQ=39)
      gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCGGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..ggcct                                                                                       <  1:2529884/101‑5 (MQ=38)
       atcgTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..ggcctt                                                                                      <  2:1325268/98‑6 (MQ=38)
        tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GGGCTGCCCCCCTC.TT..ggccttc                                                                                     >  1:2651264/1‑95 (MQ=255)
        tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..ggccttc                                                                                     <  2:5321246/101‑7 (MQ=38)
        tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..ggccttc                                                                                     <  2:3970458/101‑7 (MQ=38)
        tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..ggccttc                                                                                     <  1:2877183/101‑7 (MQ=38)
         ggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..ggccttcg                                                                                    <  2:3567438/101‑8 (MQ=38)
         ggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..ggccttcg                                                                                    <  1:1939217/101‑8 (MQ=38)
            cGATGCGATCGAGTGGCCCGAGGAGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGcc                                                                                 <  1:7075642/101‑1 (MQ=25)
              aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCaa                                                                               <  2:6178110/101‑1 (MQ=37)
              aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCaa                                                                               <  1:4447130/101‑1 (MQ=37)
                 cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GGGCTGCCCCCCTC.CTTGG..CCTCGGCCCATTc                                                                            >  2:5970852/1‑101 (MQ=255)
                 cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGGCAATTc                                                                            >  2:1172561/1‑101 (MQ=37)
                 cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTc                                                                            <  1:6691494/101‑1 (MQ=37)
                 cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTc                                                                            <  1:638906/101‑1 (MQ=37)
                  gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCa                                                                           <  2:2734722/101‑1 (MQ=37)
                   aTCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCCATTCAc                                                                          >  2:2997201/1‑101 (MQ=37)
                    tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCAcc                                                                         <  2:1401243/101‑1 (MQ=37)
                       aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCAccccc                                                                      <  1:1752097/101‑1 (MQ=25)
                       aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCAccccc                                                                      <  2:429156/101‑1 (MQ=25)
                         tGGCCCGAGGGGTTCTGCCGTCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCGCCg                                                                    <  1:5601447/101‑1 (MQ=25)
                          ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTCTTT..GC.CTTCGGCCACTTCACCCCCCgg                                                                   >  2:6068422/1‑101 (MQ=25)
                           gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCggg                                                                  <  2:197939/101‑1 (MQ=25)
                             ccGAGGGGTTCTGCCGCCGCGACATCGGGTGGGGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCACCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCcgggt                                                                >  1:5058475/1‑99 (MQ=21)
                             ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGat                                                                <  1:5714711/101‑1 (MQ=25)
                               gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGatat                                                              <  2:2957829/101‑1 (MQ=25)
                               gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.CTGGG..CCTCAGCCAATTCACCCCCCCGGacat                                                              >  2:7015512/1‑101 (MQ=37)
                                gggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATAtt                                                             <  2:4684788/100‑1 (MQ=21)
                                aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCCATTCACCCCCCCGGCTAtt                                                             >  1:6354303/1‑101 (MQ=21)
                                aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATAtt                                                             <  1:1885179/101‑1 (MQ=21)
                                 ggtgTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATAttt                                                            <  2:6633905/98‑1 (MQ=21)
                                  gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCCC.TT..GGCCTTCGGCCAATTCACCCCCCCGGGTATTTc                                                           >  2:4847750/1‑101 (MQ=21)
                                      tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGc                                                       <  2:3431092/101‑1 (MQ=21)
                                      tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCCGGATATTTCGGGc                                                       >  2:3216138/1‑101 (MQ=21)
                                         gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  1:2044285/101‑1 (MQ=21)
                                         gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  1:2904986/101‑1 (MQ=21)
                                           cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:3051239/101‑1 (MQ=21)
                                           cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGACATTTCGGGCAAGaa                                                  >  2:173144/1‑100 (MQ=21)
                                            gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:3678331/101‑1 (MQ=21)
                                            gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCCGGATAATTCGGGGCAGATa                                                 >  2:2562912/1‑101 (MQ=18)
                                             ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                                <  1:887977/101‑1 (MQ=21)
                                             ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                                <  2:723611/101‑1 (MQ=21)
                                              cgcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaaa                                               <  2:6787163/101‑1 (MQ=21)
                                              caggACAGCGGCTGGCGCGCGCTCTGAGACGAAAC.GG.GGTGCGCTGCCCCCCTC.TT..GTCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaaa                                               <  2:2147433/98‑1 (MQ=21)
                                               gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:6537118/101‑1 (MQ=21)
                                               gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:303436/101‑1 (MQ=21)
                                                cgACATCTTCTGGCGCGCGCTCTGAGACGAAAC.GGTGT.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGTATATTTCGGGCAATATAaaga                                             <  1:4544027/101‑1 (MQ=18)
                                                cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                             <  1:1616270/101‑1 (MQ=21)
                                                  aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaag                                           >  1:3928096/1‑101 (MQ=21)
                                                  aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaag                                           <  2:1003257/101‑1 (MQ=21)
                                                  aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCCGGATATTTCGGGGAAGATAaagaag                                           >  2:2384387/1‑101 (MQ=21)
                                                   cATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                          >  1:5277468/1‑101 (MQ=21)
                                                      cGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                                       >  2:5643539/1‑101 (MQ=21)
                                                      cGGCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                                       <  2:3860472/101‑1 (MQ=21)
                                                       ggCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAAATAAAGAAGGTCtt                                      >  1:3659943/1‑101 (MQ=21)
                                                        gCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAATAAGGTCttt                                     <  2:6957642/101‑1 (MQ=21)
                                                        gCTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                     <  1:4306564/101‑1 (MQ=21)
                                                         cTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  2:6361469/101‑1 (MQ=21)
                                                         cTGGCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  2:3372018/101‑1 (MQ=21)
                                                           ggCGCGCGCTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTcc                                  <  1:4011077/101‑1 (MQ=21)
                                                            gcgcgcgcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 >  2:2787478/1‑101 (MQ=21)
                                                              gcgcgcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:5011686/101‑1 (MQ=21)
                                                              gcgcgcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:4426601/101‑1 (MQ=21)
                                                              gagcgcTCTGAGACGTGAC.GGCCG.TCCCTGCCGCCCGC.TT..GGCCTTCGGCCAAATCGCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTATGTACCttt                               <  2:7052738/99‑1 (MQ=12)
                                                               cgcgcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  1:4705344/101‑1 (MQ=21)
                                                                gcgcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  1:2229172/101‑1 (MQ=21)
                                                                 cgcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                            <  1:533247/101‑1 (MQ=21)
                                                                 cgcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                            >  1:2331370/1‑101 (MQ=21)
                                                                 cgcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                            <  1:4144111/101‑1 (MQ=21)
                                                                  gcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:5448908/101‑1 (MQ=21)
                                                                  gcTCTGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  2:1218978/101‑1 (MQ=21)
                                                                   ctctGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                          <  1:1014956/101‑1 (MQ=21)
                                                                     ctGAGACGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  1:3082905/101‑1 (MQ=21)
                                                                       gagaCGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      <  2:2858374/101‑1 (MQ=21)
                                                                         gaCGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                    <  2:4081874/101‑1 (MQ=21)
                                                                          gcgagccggg.gg.gCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                   <  1:5422155/93‑1 (MQ=21)
                                                                           cGAAACGGG.GG.GCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                  >  1:1230074/1‑101 (MQ=21)
                                                                             aagcggg.gg.gCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                <  2:5768779/96‑1 (MQ=21)
                                                                             aaacggg.gg.gCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                <  2:2494965/96‑1 (MQ=21)
                                                                              ggcagg.gg.gCGCTG.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgg         >  2:1872811/5‑101 (MQ=0)
                                                                               gcagg.gg.gCGCTG.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg        >  1:7009280/4‑101 (MQ=0)
                                                                               acggg.gg.gCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat              <  2:650406/98‑1 (MQ=25)
                                                                                cagg.gg.gCGCTG.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGTCCCTTCATCTTGCCGCAAATATCCCgggg       >  1:4996583/3‑101 (MQ=2)
                                                                                cagg.gg.gCGCTG.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg       >  2:769997/3‑101 (MQ=0)
                                                                                 ggg.gg.gCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            >  2:5529951/2‑101 (MQ=25)
                                                                                 ggg.gg.gCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  1:5920528/100‑1 (MQ=25)
                                                                                 ggg.gg.gCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            >  1:2138245/2‑101 (MQ=25)
                                                                                 agg.gg.gCGCTG.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:3838067/2‑101 (MQ=0)
                                                                                  gggcg.gCGCTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTGTCATCTTGCCCGAAATATcc            >  1:6847284/5‑101 (MQ=25)
                                                                                  gg.gg.gCGCTG.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggg     >  1:1206135/1‑100 (MQ=0)
                                                                                         cgcTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggggg      >  2:2631240/1‑101 (MQ=37)
                                                                                         cgcTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggggg      <  1:5190003/101‑1 (MQ=37)
                                                                                           cTGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGac    >  1:3493723/1‑100 (MQ=37)
                                                                                            tGCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGc   >  2:1128118/1‑101 (MQ=37)
                                                                                             gCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCTCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGcc  <  2:5257869/101‑1 (MQ=37)
                                                                                             gCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGcc  >  1:781716/1‑101 (MQ=37)
                                                                                             gCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGcc  >  1:5758333/1‑101 (MQ=25)
                                                                                                       |                                                                                             
CTATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GG.GG.GCGCTGCCCCCCTC.CT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566456‑2566644

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG