Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,625 C→G 4.0% intergenic (+97/+97) 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,6250C→G 4.0% 13.2 ‑1.7 0.0259intergenic (+97/+97)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (10/1):  ref base (175/73):  total (185/74)
Fisher's exact test strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

AGACGAAAC.GGGGGCGCTGCCCCCCTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAG  >  NC_014034/2566528‑2566725
                                                                                                    |                                                                                                    
 gaCGAAACGGGGGGCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                                                     <  1:59639/101‑1 (MQ=21)
       gcagggggCGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg                                                                                         >  1:4762757/4‑101 (MQ=0)
        cggggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                                                                              <  1:3815916/99‑1 (MQ=25)
        cagggggCGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                                                                                        >  1:728053/3‑101 (MQ=0)
         agggggCGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                                                       >  1:4342168/2‑101 (MQ=0)
          gggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                                            <  1:2242135/101‑1 (MQ=25)
          gggggCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGCATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATtcc                                                                                            <  1:125582/101‑1 (MQ=25)
          gggggCGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggg                                                                                      >  2:3245842/1‑100 (MQ=0)
           ggggCGATGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg                                                                                           <  1:918241/101‑1 (MQ=25)
              gcgcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                        >  1:5236898/1‑101 (MQ=37)
                gcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                                                                      <  2:3056052/101‑1 (MQ=37)
                gcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                                                                      <  1:5498512/101‑1 (MQ=37)
                   gCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGAAATTTTGGGGAAGATAAAGAAGGGCTTTTTCTTTTATCTTTCCCCAAAAATCCCCGGGGAGcc                                                                                   >  1:1184556/1‑101 (MQ=16)
                     ccctcttggg..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                                               >  2:5275561/10‑101 (MQ=255)
                       ccctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                                               >  2:5454605/11‑101 (MQ=255)
                        cctcttggccTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCtt                                                                              <  1:3293446/92‑1 (MQ=38)
                        cctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCtt                                                                              <  2:4282080/92‑1 (MQ=255)
                         ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                                                                             <  1:856547/93‑1 (MQ=255)
                         ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                                                                             <  2:4070968/93‑1 (MQ=255)
                           cttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                                                                           >  2:4229130/7‑101 (MQ=255)
                            ttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                                          <  1:1478804/96‑1 (MQ=255)
                               gccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcg                                                                       <  2:3474278/99‑1 (MQ=255)
                               gccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcg                                                                       >  1:2428521/3‑101 (MQ=255)
                                 cTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGcgggcggg                                                                     >  2:3715136/1‑101 (MQ=255)
                                    cGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                                                                  <  2:1296680/101‑1 (MQ=255)
                                     ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGcg                                                                 >  1:3883801/1‑100 (MQ=255)
                                     ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa                                                                 >  1:4274821/1‑101 (MQ=255)
                                       ccAATTCAGCGCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTCACCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgc                                                               <  2:4643456/101‑1 (MQ=39)
                                       ccAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgc                                                               <  2:2187833/101‑1 (MQ=255)
                                        cAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcg                                                              >  2:2031281/1‑101 (MQ=255)
                                            tCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGACCGCTTGCGGGCGGGGGCAGCGcccc                                                          >  1:2170162/1‑101 (MQ=255)
                                             cACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGGCGCGccccc                                                         >  1:4183811/1‑101 (MQ=255)
                                              cccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg                                                        <  1:1778636/100‑1 (MQ=255)
                                              cccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg                                                        <  2:1220946/100‑1 (MQ=255)
                                               ccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTTATCTTTGCCGGAATATCCCCGGGGAGGCCGCTTGCGGGCGGGGGGAGGGCCCCCGa                                                       >  2:2073563/1‑101 (MQ=37)
                                                  cccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgc                                                    >  1:4069280/1‑101 (MQ=255)
                                                   ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                                                   >  1:3105209/1‑101 (MQ=255)
                                                   ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                                                   >  1:4391008/1‑101 (MQ=255)
                                                   ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                                                   >  1:1950834/1‑101 (MQ=255)
                                                     gggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                                                 >  2:5672985/1‑101 (MQ=255)
                                                      ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                                                >  1:2898764/1‑101 (MQ=255)
                                                       gATATTTCGGGCGAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCGGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                               <  1:5366197/101‑1 (MQ=255)
                                                        atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                                              >  1:5169160/1‑101 (MQ=255)
                                                        atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                                              <  2:3680679/101‑1 (MQ=255)
                                                          atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                                            >  1:2858485/1‑101 (MQ=255)
                                                            ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCTAAGCGCCCGAcgcg                                          <  2:1672394/101‑1 (MQ=255)
                                                            ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                                          >  2:2735754/1‑101 (MQ=255)
                                                             tCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgc                                         >  1:4123738/1‑101 (MQ=255)
                                                              cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg                                        >  2:1490716/1‑101 (MQ=255)
                                                               gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg                                       >  1:3945630/1‑101 (MQ=255)
                                                                ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                                      >  1:2011620/1‑101 (MQ=255)
                                                                  aaaGATAAAGAAGATCTTTTCCTTTCATCTATCACTAAGTATCGGGGGGGAGCCCGCTTGGGGGCGGGGGCAGCGCCCGCGAAGCGCCCGACGCGcggcgg                                    <  2:4162318/100‑1 (MQ=255)
                                                                   aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  2:86709/1‑101 (MQ=255)
                                                                   aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  2:1619183/1‑101 (MQ=255)
                                                                   aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  2:2866536/1‑101 (MQ=255)
                                                                   aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  2:3334500/1‑101 (MQ=255)
                                                                   aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  1:1223751/1‑101 (MQ=255)
                                                                   aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  1:3483321/1‑101 (MQ=255)
                                                                    aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc                                  >  1:3559271/1‑101 (MQ=255)
                                                                    aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc                                  >  1:1707415/1‑101 (MQ=255)
                                                                     gATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCt                                 <  2:258023/101‑1 (MQ=255)
                                                                     aaTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCt                                 <  2:1089292/100‑1 (MQ=255)
                                                                      aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                <  1:2626191/101‑1 (MQ=255)
                                                                      aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                >  2:4893439/1‑101 (MQ=255)
                                                                       tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGCCGCGCGGCGGGCttt                               >  2:533794/1‑100 (MQ=255)
                                                                       tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg                               >  1:956784/1‑101 (MQ=255)
                                                                       tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg                               >  2:1837615/1‑101 (MQ=255)
                                                                       tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg                               >  1:2661004/1‑101 (MQ=255)
                                                                        aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                              >  1:1602930/1‑101 (MQ=255)
                                                                         aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGACGCGCCCGACGCGGGGCGGGCTTccc                             >  1:4545012/1‑101 (MQ=255)
                                                                         aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc                             >  2:61536/1‑101 (MQ=255)
                                                                         aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCCCCCCCGAAGCGCCCGCCGCGCGGGGGGCTTgcc                             >  1:3058712/1‑101 (MQ=255)
                                                                          agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc                            <  2:5168874/101‑1 (MQ=255)
                                                                          agaagGTCTTTTCCTTTCATCTTGCCCGAAATACCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc                            <  1:4984098/101‑1 (MQ=255)
                                                                            aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc                          >  2:5650168/1‑101 (MQ=255)
                                                                            aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc                          >  1:211319/1‑101 (MQ=255)
                                                                            aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc                          >  2:3984266/1‑101 (MQ=255)
                                                                             agGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgcc                         >  1:4711779/1‑101 (MQ=255)
                                                                              ggTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgccc                        >  1:1352507/1‑101 (MQ=255)
                                                                                tCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa                      <  2:3883559/101‑1 (MQ=255)
                                                                                tCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa                      >  2:115668/1‑101 (MQ=255)
                                                                                  ttttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagg                    >  2:4691268/1‑101 (MQ=255)
                                                                                    ttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGTTGCCCGCCCAaggag                  >  2:5068179/1‑101 (MQ=255)
                                                                                    ttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag                  >  2:4422187/1‑101 (MQ=255)
                                                                                     tCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAATCGCCCGACGACCGGCGGCACTTCGCGCCCAaggcgg                 >  2:283892/1‑101 (MQ=25)
                                                                                     tCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCCCTTGCGGGCGGGGGCACCGCCCCCGCGGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagaggg                 >  1:1837167/1‑97 (MQ=37)
                                                                                      ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                >  2:679215/1‑101 (MQ=255)
                                                                                      ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                >  2:2658303/1‑101 (MQ=255)
                                                                                      ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCCCCCAAGGAggg                >  2:4593039/1‑101 (MQ=255)
                                                                                       cTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc               >  2:4965464/1‑101 (MQ=255)
                                                                                       cTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc               <  1:3715364/101‑1 (MQ=255)
                                                                                       cTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc               >  2:2969681/1‑101 (MQ=255)
                                                                                        tttCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGAGGCGGGCTTGCCCGCCCAAGGAGGgcg              >  1:5413910/1‑101 (MQ=255)
                                                                                           cATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCaa           >  2:3155685/1‑101 (MQ=255)
                                                                                           cATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCaa           >  1:2267701/1‑101 (MQ=255)
                                                                                              cTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCg        >  1:5667199/1‑101 (MQ=255)
                                                                                               ttGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGc       >  1:3673989/1‑101 (MQ=255)
                                                                                                tGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGcc      <  2:4729318/101‑1 (MQ=255)
                                                                                                  cccGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTc    >  1:1035467/1‑101 (MQ=255)
                                                                                                   ccGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTca   >  2:5390713/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
AGACGAAAC.GGGGGCGCTGCCCCCCTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAG  >  NC_014034/2566528‑2566725

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