Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,557 +C 98.4% 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.→C98.4% 3.7 ‑2.6 ‑0.4255intergenic (+29/+165)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (71/180):  ref base (3/1):  total (74/181)
Fisher's exact test strand distribution p-value = 7.49e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.18e-01

Only 100 of 298 total aligned reads displayed.

ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GG.GG.GCGCT.GCCCCCCTC.CT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566458‑2566644
                                                                                                            |                                                                                        
aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..g                                                                                           <  2:2247535/101‑1 (MQ=255)
aTGCGATGGCCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..g                                                                                           <  1:6537275/101‑1 (MQ=255)
 tGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..gg                                                                                          <  2:2010294/101‑2 (MQ=255)
    gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCGGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..ggcct                                                                                       <  1:2529884/101‑5 (MQ=38)
     atcgTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..ggcctt                                                                                      <  2:1325268/98‑6 (MQ=38)
      tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GGGCT.GCCCCCCTC.TT..ggccttc                                                                                     >  1:2651264/1‑95 (MQ=255)
      tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..ggccttc                                                                                     <  2:5321246/101‑7 (MQ=38)
      tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..ggccttc                                                                                     <  2:3970458/101‑7 (MQ=38)
      tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..ggccttc                                                                                     <  1:2877183/101‑7 (MQ=38)
       ggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..ggccttcg                                                                                    <  2:3567438/101‑8 (MQ=38)
       ggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..ggcctccg                                                                                    <  1:1939216/101‑8 (MQ=38)
            aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCaa                                                                               <  2:6178110/101‑1 (MQ=37)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GGGCT.GCCCCCCTC.CTTGG..CCTCGGCCCATTc                                                                            >  2:5970852/1‑101 (MQ=255)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGGCAATTc                                                                            >  2:1172561/1‑101 (MQ=37)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTc                                                                            <  1:5837617/101‑1 (MQ=37)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTc                                                                            <  2:1038677/101‑1 (MQ=37)
                gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GGGCTGGCCCCCCCT.TT..GC.CTTCGGCCAATTcc                                                                           >  1:6027130/1‑100 (MQ=255)
                gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCa                                                                           <  2:2734722/101‑1 (MQ=37)
                 aTCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCCATTCAc                                                                          >  2:2997201/1‑101 (MQ=37)
                  tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCAcc                                                                         <  2:1401243/101‑1 (MQ=37)
                     aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCAccccc                                                                      <  2:429156/101‑1 (MQ=25)
                     aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCAccccc                                                                      <  1:3988735/101‑1 (MQ=25)
                        ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTCTTT..GC.CTTCGGCCACTTCACCCCCCgg                                                                   >  2:6068422/1‑101 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGatat                                                              <  2:2957829/101‑1 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.CTGGG..CCTCAGCCAATTCACCCCCCCGGacat                                                              >  2:7015512/1‑101 (MQ=37)
                              gggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATAtt                                                             <  2:4684788/100‑1 (MQ=21)
                               ggtgTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATAttt                                                            <  2:6633905/98‑1 (MQ=21)
                                gggTTGTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTc                                                           <  1:5731627/101‑1 (MQ=21)
                                gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TTTGG..CCTCGGCCAATTCACCCCCCCGGATCTTTc                                                           >  2:1133395/1‑101 (MQ=25)
                                gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCCC.TT..GGCCTTCGGCCAATTCACCCCCCCGGGTATTTc                                                           >  2:4847750/1‑101 (MQ=21)
                                 ggTTCTGCCGCCGCGACATCGGCTGGGGGGCGCTCTGAGACGAAAAG.GG.GG.GGGCT.GCCCCCCCC.TT..GGCCTTCGGCCCATTCCCCCCCCCGCGTATTTCg                                                          >  1:4713374/1‑101 (MQ=255)
                                 ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCg                                                          <  1:6697804/101‑1 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGc                                                       <  2:3431092/101‑1 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCCGGATATTTCGGGc                                                       >  2:3216138/1‑101 (MQ=21)
                                     cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCa                                                      <  1:463114/101‑1 (MQ=21)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTCTTT..GC.CTTCGGCCAATTCACCCCCCCGGATATTTCGGGCAAg                                                    >  1:1795964/1‑101 (MQ=21)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  1:5214826/101‑1 (MQ=21)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  1:2904986/101‑1 (MQ=21)
                                        ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGa                                                   <  1:1342072/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGAt                                                  <  2:5073261/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGACATTTCGGGCAAGaa                                                  >  2:173144/1‑100 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  1:5047668/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:3678331/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  1:3848319/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCCGGATAATTCGGGGCAGATa                                                 >  2:2562912/1‑101 (MQ=18)
                                           ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                                <  2:723611/101‑1 (MQ=21)
                                            cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGTGATATTTCGGGCAAGATaaa                                               <  1:6385428/101‑1 (MQ=21)
                                            cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaaa                                               <  1:1607458/101‑1 (MQ=21)
                                            cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaaa                                               <  2:6787163/101‑1 (MQ=21)
                                            caggACAGCGGCTGGCGCGCGCTCTGAGACG.AAAC.GG.GGTGCGCT.GCCCCCCTC.TT..GTCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaaa                                               <  2:2147433/98‑1 (MQ=21)
                                             gcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:5355032/101‑1 (MQ=21)
                                             gcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:303436/101‑1 (MQ=21)
                                                aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCGAGATAaagaag                                           <  2:317249/101‑1 (MQ=21)
                                                aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaag                                           <  1:5076772/101‑1 (MQ=21)
                                                aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCCGGATATTTCGGGGAAGATAaagaag                                           >  2:2384387/1‑101 (MQ=21)
                                                  aTCGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGt                                         <  1:4425903/101‑1 (MQ=21)
                                                    cGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                                       >  2:5643539/1‑101 (MQ=21)
                                                    cGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                                       <  1:3665733/101‑1 (MQ=21)
                                                    cGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCt                                       <  2:3860472/101‑1 (MQ=21)
                                                    cGGCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGAAATTTCGAGCAAGATAAAGAAGGACt                                       >  1:722852/1‑101 (MQ=18)
                                                     ggCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  1:4252317/101‑1 (MQ=21)
                                                      gCTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAATAAGGTCttt                                     <  2:6957642/101‑1 (MQ=21)
                                                       cttgcgcgctcTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GTCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  1:894438/98‑1 (MQ=21)
                                                       cTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  2:6361469/101‑1 (MQ=21)
                                                       cTGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  2:3372018/101‑1 (MQ=21)
                                                        tGGCGCGCGCTCTGAGACG.AAACGGG.GG.GCGGT.GCCCCCCTC.TT..GGCCTTCGGGCAATTCACCCCCCGGGAAATTTCGGCGAAGATAAAGAAGGCCTTTTc                                   >  1:6946152/1‑101 (MQ=18)
                                                          gcgcgcgcTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                 >  2:2787478/1‑101 (MQ=21)
                                                            gcgcgcTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:4426601/101‑1 (MQ=21)
                                                            gagcgcTCTGAGACG.TGAC.GGCCG.TCCCT.GCCGCCCGC.TT..GGCCTTCGGCCAAATCGCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTATGTACCttt                               <  2:7052738/99‑1 (MQ=12)
                                                              gcgcTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  1:2229172/101‑1 (MQ=21)
                                                               cgcTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                            >  1:2331370/1‑101 (MQ=21)
                                                                gcTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  2:1218978/101‑1 (MQ=21)
                                                                gcTCTGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           >  1:7070547/1‑101 (MQ=21)
                                                                 ctctGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                          <  1:1014956/101‑1 (MQ=21)
                                                                   ctGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  1:3082905/101‑1 (MQ=21)
                                                                   ctGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  1:3751797/101‑1 (MQ=21)
                                                                    tGAGACG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                       >  1:1102326/1‑101 (MQ=21)
                                                                     gagaCG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      <  2:2858374/101‑1 (MQ=21)
                                                                       gaCG.AAACGGG.GG.GCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                    <  2:4081874/101‑1 (MQ=21)
                                                                            aagcggg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                <  2:5768779/96‑1 (MQ=21)
                                                                            aaacggg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                <  2:2494965/96‑1 (MQ=21)
                                                                            aaacggg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                >  1:6475801/6‑101 (MQ=21)
                                                                            aaacggg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                >  1:5450484/6‑101 (MQ=21)
                                                                             ggcagg.gg.gCGCT.G.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgg         >  2:1872811/5‑101 (MQ=0)
                                                                             aacggg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata               >  1:4021229/5‑101 (MQ=25)
                                                                              acggg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat              <  2:650406/98‑1 (MQ=25)
                                                                               cagg.gg.gCGCT.G.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg       >  2:769997/3‑101 (MQ=0)
                                                                                ggg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            >  2:5529951/2‑101 (MQ=25)
                                                                                ggg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  1:5317831/100‑1 (MQ=25)
                                                                                ggg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            >  1:2138245/2‑101 (MQ=25)
                                                                                agg.gg.gCGCT.G.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:5424946/2‑101 (MQ=0)
                                                                                agg.gg.gCGCT.G.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:1957994/2‑101 (MQ=0)
                                                                                 gg.gg.gCGCT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc           >  1:796048/1‑101 (MQ=25)
                                                                                    gg.gCGCT.G.CCCCCTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag   >  1:1654507/1‑98 (MQ=2)
                                                                                        cgcT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggggg      >  1:4471253/1‑101 (MQ=25)
                                                                                        cgcT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggggg      >  2:2631240/1‑101 (MQ=37)
                                                                                        cgcT.GCCCCCCTC.TT..GGCCTCCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgcggg      <  1:5614949/101‑1 (MQ=25)
                                                                                         gcT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa     >  1:1211124/1‑101 (MQ=37)
                                                                                          cT.GCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGac    >  1:3493723/1‑100 (MQ=37)
                                                                                             gCCCCCCTC.TT..GGCCTTCGGCCAATTCACCCCTCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGcc  <  2:5257869/101‑1 (MQ=37)
                                                                                                            |                                                                                        
ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GG.GG.GCGCT.GCCCCCCTC.CT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566458‑2566644

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG