Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,626 G→C 3.9% intergenic (+98/+96) 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,6260G→C 3.9% 15.2 ‑2.7 ‑0.3259intergenic (+98/+96)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (10/0):  ref base (178/71):  total (188/71)
Fisher's exact test strand distribution p-value = 6.62e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.04e-01

Only 100 of 286 total aligned reads displayed.

GACGAAACGGGGGCGCTGCCCCCCTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAGC  >  NC_014034/2566529‑2566726
                                                                                                   |                                                                                                      
       ggggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                                                >  1:701147/2‑101 (MQ=25)
        gggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGCATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATtcc                                                                                               <  1:125582/101‑1 (MQ=25)
        gggggCGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCC.gggggg                                                                                        >  2:3245842/1‑100 (MQ=0)
              gcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGa                                                                                        <  1:728443/101‑1 (MQ=37)
              gcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGa                                                                                        <  2:3056052/101‑1 (MQ=37)
               cTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAg                                                                                       >  1:4357674/1‑101 (MQ=37)
                 gCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGAAATTTTGGGGAAGATAAAGAAGGGCTTTTTCTTTTATCTTTCCCCAAAAATCCC.CGGGGAGcc                                                                                     >  1:1184556/1‑101 (MQ=16)
                   ccctcttggg..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCt                                                                                 >  2:5275561/10‑101 (MQ=255)
                     ccctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCt                                                                                 >  2:5454605/11‑101 (MQ=255)
                     ccctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCt                                                                                 >  1:5486878/11‑101 (MQ=255)
                      cctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCtt                                                                                <  2:4282080/92‑1 (MQ=255)
                       ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.g                                                                              <  2:4070968/93‑1 (MQ=255)
                       ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.g                                                                              <  1:3358734/93‑1 (MQ=255)
                       ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.g                                                                              >  1:5383904/9‑101 (MQ=255)
                         cttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.gcg                                                                            >  2:4229130/7‑101 (MQ=255)
                          ttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.gcgg                                                                           <  1:1478804/96‑1 (MQ=255)
                           tggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.gcggg                                                                          >  1:3680907/5‑101 (MQ=255)
                             gccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.gcgggcg                                                                        <  2:3474278/99‑1 (MQ=255)
                               cTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.Gcgggcggg                                                                      >  2:3715136/1‑101 (MQ=255)
                                  cGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGc                                                                   <  2:1296680/101‑1 (MQ=255)
                                   ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGcc                                                                  >  1:5553447/1‑100 (MQ=255)
                                   ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCa                                                                  >  1:2246156/1‑101 (MQ=255)
                                     ccAATTCAGCGCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTCACCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAgc                                                                <  2:4643456/101‑1 (MQ=39)
                                     ccAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAgc                                                                <  2:1114626/101‑1 (MQ=255)
                                     ccAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAgc                                                                <  2:2187833/101‑1 (MQ=255)
                                      cAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAgcg                                                               >  2:2031281/1‑101 (MQ=255)
                                          tCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGcccc                                                           >  1:3896590/1‑101 (MQ=255)
                                          tCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGACCGCTT.GCGGGCGGGGGCAGCGcccc                                                           >  1:2170162/1‑101 (MQ=255)
                                           cACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGccccc                                                          <  2:1181547/101‑1 (MQ=255)
                                            cccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCg                                                         <  1:1778636/100‑1 (MQ=255)
                                             ccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTTATCTTTGCCGGAATATCCC.CGGGGAGGCCGCTT.GCGGGCGGGGGGAGGGCCCCCGa                                                        >  2:2073563/1‑101 (MQ=37)
                                                 ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAgcg                                                    >  1:1950834/1‑101 (MQ=255)
                                                  cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAgcgc                                                   >  1:5701705/1‑101 (MQ=255)
                                                   gggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGcc                                                  >  2:5672985/1‑101 (MQ=255)
                                                     gATATTTCGGGCGAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCGGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                                <  1:5366197/101‑1 (MQ=255)
                                                     gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                                >  1:4202269/1‑101 (MQ=255)
                                                      atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                                               <  2:3680679/101‑1 (MQ=255)
                                                        ttttTGGGCAAAAAAAAGGAGGTTTTTTCCTTTTCTTTTTCCCGAAATTCCCCGGGGGGACCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAc                                              <  2:1074223/100‑1 (MQ=19)
                                                        atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                                             >  1:5176522/1‑101 (MQ=255)
                                                        atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                                             >  1:923982/1‑101 (MQ=255)
                                                        atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                                             >  1:3745703/1‑101 (MQ=255)
                                                          ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCTAAGCGCCCGAcgcg                                           <  2:1672394/101‑1 (MQ=255)
                                                          ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                                           >  2:2735754/1‑101 (MQ=255)
                                                            cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg                                         >  2:1490716/1‑101 (MQ=255)
                                                              ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                                       >  1:4405431/1‑101 (MQ=255)
                                                              ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                                       >  1:3870149/1‑101 (MQ=255)
                                                              ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                                       >  1:2011620/1‑101 (MQ=255)
                                                              ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGC.GGGGCAGGGGCCCCCGAGCCCCCGACGCccggcg                                      >  1:3620813/1‑101 (MQ=37)
                                                                cAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGcggcgg                                     >  1:3238479/1‑101 (MQ=255)
                                                                aaaGATAAAGAAGATCTTTTCCTTTCATCTATCACTAAGTATCGG.GGGGGAGCCCGCTT.GGGGGCGGGGGCAGCGCCCGCGAAGCGCCCGACGCGcggcgg                                     <  2:4162318/100‑1 (MQ=255)
                                                                 aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                    >  2:1619183/1‑101 (MQ=255)
                                                                 aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                    >  2:2866536/1‑101 (MQ=255)
                                                                 aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                    >  2:3334500/1‑101 (MQ=255)
                                                                 aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                    >  2:86709/1‑101 (MQ=255)
                                                                 aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                    >  1:1223751/1‑101 (MQ=255)
                                                                  aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc                                   >  1:1707415/1‑101 (MQ=255)
                                                                   gATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCt                                  <  2:258023/101‑1 (MQ=255)
                                                                    aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                 >  2:4893439/1‑101 (MQ=255)
                                                                     tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGCCGCGCGGCGGGCttt                                >  2:533794/1‑100 (MQ=255)
                                                                     tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg                                >  2:1837615/1‑101 (MQ=255)
                                                                      aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCTCCCGACGCGCGGCGGGCTTgc                               >  1:3541890/1‑101 (MQ=255)
                                                                      aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                               >  1:870992/1‑101 (MQ=255)
                                                                      aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                               >  1:1602930/1‑101 (MQ=255)
                                                                       aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc                              >  2:61536/1‑101 (MQ=255)
                                                                        agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc                             <  2:5168874/101‑1 (MQ=255)
                                                                         aaatgttttttTCTTTTCTCCTGCCCGAAATTCCCG.GGGGGCGCCCGGTTGGGGGGGGGGGGCGCCCCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc                             <  1:4893720/94‑1 (MQ=12)
                                                                          aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGTTgcccgc                           >  1:2628142/1‑101 (MQ=255)
                                                                          aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc                           >  1:211319/1‑101 (MQ=255)
                                                                          aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc                           >  2:3984266/1‑101 (MQ=255)
                                                                          aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc                           >  2:5650168/1‑101 (MQ=255)
                                                                            ggTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgccc                         >  1:4299947/1‑101 (MQ=255)
                                                                            ggTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgccc                         >  1:1352507/1‑101 (MQ=255)
                                                                              tCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa                       <  2:3883559/101‑1 (MQ=255)
                                                                               cTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaac                      >  1:2882303/1‑100 (MQ=255)
                                                                                ttttCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagg                     >  1:4083803/1‑101 (MQ=255)
                                                                                ttttCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagg                     >  2:4691268/1‑101 (MQ=255)
                                                                                  ttCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGTTGCCCGCCCAaggag                   >  2:5068179/1‑101 (MQ=255)
                                                                                  ttCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag                   >  1:4159021/1‑101 (MQ=255)
                                                                                  ttCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag                   >  2:4422187/1‑101 (MQ=255)
                                                                                   tCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAATCGCCCGACGACCGGCGGCACTTCGCGCCCAaggcgg                  >  2:283892/1‑101 (MQ=25)
                                                                                   tCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCCCTT.GCGGGCGGGGGCACCGCCCCCGCGGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagaggg                  >  1:1837167/1‑97 (MQ=37)
                                                                                    ccTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                 >  2:2658303/1‑101 (MQ=255)
                                                                                    ccTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                 >  2:679215/1‑101 (MQ=255)
                                                                                    ccTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCCCCCAAGGAggg                 >  2:4593039/1‑101 (MQ=255)
                                                                                     cTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc                >  1:2347856/1‑101 (MQ=255)
                                                                                     cTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc                >  2:2969681/1‑101 (MQ=255)
                                                                                     cTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc                >  2:4965464/1‑101 (MQ=255)
                                                                                     cTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc                <  1:3097904/101‑1 (MQ=255)
                                                                                       ttCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcgc              >  1:2981227/1‑101 (MQ=255)
                                                                                         cATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCaa            >  2:3155685/1‑101 (MQ=255)
                                                                                          aTCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAc           >  1:2689956/1‑101 (MQ=255)
                                                                                          aTCTGGCCCGAAATACCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAc           <  1:2451054/101‑1 (MQ=255)
                                                                                            cTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCg         >  1:5127850/1‑101 (MQ=255)
                                                                                            cTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCg         >  1:403816/1‑101 (MQ=255)
                                                                                              tGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGcc       <  2:4729318/101‑1 (MQ=255)
                                                                                               gCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCt      <  1:4730459/101‑1 (MQ=255)
                                                                                                cccGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTc     >  1:1035467/1‑101 (MQ=255)
                                                                                                 ccGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTca    >  2:5390713/1‑101 (MQ=255)
                                                                                                  cGAAATATCCC.GGGGGAGCCCGCGT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTcag   >  1:4656463/1‑101 (MQ=255)
                                                                                                   gAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTcagc  >  2:1239035/1‑101 (MQ=255)
                                                                                                   |                                                                                                      
GACGAAACGGGGGCGCTGCCCCCCTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCC.GGGGGAGCCCGCTT.GCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAGC  >  NC_014034/2566529‑2566726

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