Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,557 +C 97.6% intergenic (+29/+165) purD → / ← sdaA phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,566,5572.→C97.6% 4.2 ‑2.7 ‑0.7168intergenic (+29/+165)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (44/120):  ref base (3/1):  total (47/121)
Fisher's exact test strand distribution p-value = 6.68e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.13e-01

Only 100 of 200 total aligned reads displayed.

ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGG.GCGCTGCCCCC.CTC.CT..G..C.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCT  >  NC_014034/2566458‑2566648
                                                                                                           |                                                                                               
aTGCGATGGTCGATGCGATCGCCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..g                                                                                                  <  1:1875668/101‑1 (MQ=255)
aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGACCCC.CTC.TT..g                                                                                                  <  1:4335753/101‑1 (MQ=255)
    gATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggcct                                                                                              <  2:5687034/101‑5 (MQ=39)
     aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggcctt                                                                                             <  1:3585313/101‑6 (MQ=38)
     aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggcctt                                                                                             <  2:3227032/101‑6 (MQ=38)
     aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggcctt                                                                                             <  1:4189479/101‑6 (MQ=38)
      tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGGTGCGCTGCCCCC.CTC.TT..ggccttc                                                                                            <  2:3301833/101‑7 (MQ=255)
       ggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..ggccttcg                                                                                           <  2:5796223/101‑8 (MQ=38)
          cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.cc                                                                                      <  1:2947866/101‑1 (MQ=37)
              gcgatgtccTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGTCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAAtt                                                                                  <  1:5501191/93‑1 (MQ=21)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTc                                                                                 <  2:2206998/101‑1 (MQ=37)
                gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCa                                                                                <  2:679278/101‑1 (MQ=37)
                gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCa                                                                                <  1:4639120/101‑1 (MQ=37)
                   cGACTGGCCCGGGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCGCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCAccc                                                                             <  1:1555670/101‑1 (MQ=25)
                   cGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCAccc                                                                             <  2:2854771/101‑1 (MQ=37)
                    gACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCAcccc                                                                            <  2:1303849/101‑1 (MQ=25)
                      cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCAcccccc                                                                          <  1:2375308/101‑1 (MQ=25)
                         gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCggg                                                                       <  2:558639/101‑1 (MQ=25)
                          cccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.a                                                                     <  1:1160426/101‑1 (MQ=25)
                           ccaaGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGT.at                                                                    <  2:3684768/98‑1 (MQ=21)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGGGCGCGCTCTGAGACGAAAAG.GGGG.GCGCGGCCCCC.CTCTTT..G..CGTTCGG.GCACTTCACCCCCCgggcaat                                                                   >  1:464977/1‑97 (MQ=255)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.atat                                                                  <  1:739255/101‑1 (MQ=25)
                             gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.atat                                                                  <  1:4816648/101‑1 (MQ=25)
                               ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATAttt                                                                <  1:612127/101‑1 (MQ=21)
                                gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTc                                                               <  2:5321542/101‑1 (MQ=21)
                                gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTc                                                               <  2:1382023/101‑1 (MQ=21)
                                 ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCg                                                              <  2:3138014/101‑1 (MQ=21)
                                  gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGAAGAAAAG.GGGG.GCGCTGCCCCC.CCC.TTTGG..C.TTCGGCCCAATCCCCCCCCCgggtcctttgg                                                               >  2:2061801/1‑93 (MQ=16)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GGGCTGCCCCC.CCC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGc                                                           >  1:2831878/1‑101 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCgggg                                                           >  1:1413396/1‑100 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGc                                                           <  2:3323608/101‑1 (MQ=21)
                                     cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCa                                                          <  2:1914657/101‑1 (MQ=21)
                                      tGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCaa                                                         <  1:3805170/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCCATTCACCCCCCCGG.AGATTTCGGGCAAGaa                                                      >  2:1564995/1‑100 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGAt                                                      <  1:1521637/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGAt                                                      <  2:552918/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGAt                                                      <  1:2123216/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATa                                                     <  2:1150900/101‑1 (MQ=21)
                                            cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATaaa                                                   <  1:3981860/101‑1 (MQ=21)
                                               gACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaagaa                                                <  1:5218538/101‑1 (MQ=21)
                                                aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAaagaag                                               <  1:3420105/101‑1 (MQ=21)
                                                 cATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAgg                                              <  2:1941308/101‑1 (MQ=21)
                                                 cATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAgg                                              <  1:1479045/101‑1 (MQ=21)
                                                  aTCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGt                                             <  2:3351824/101‑1 (MQ=21)
                                                   tCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTc                                            <  1:2857171/101‑1 (MQ=21)
                                                    cGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCt                                           <  1:3000140/101‑1 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCtt                                          <  1:485086/101‑1 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCtt                                          <  1:5350399/101‑1 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGGCtt                                          >  2:3293246/1‑101 (MQ=21)
                                                       cTGGCGCGTGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.GCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCtttt                                        <  2:1437643/101‑1 (MQ=21)
                                                       cTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCtttt                                        <  2:155874/101‑1 (MQ=21)
                                                       cTGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCtttt                                        <  2:4869272/101‑1 (MQ=21)
                                                        tGGCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAAGTCATTTc                                       >  2:5260801/1‑101 (MQ=18)
                                                         ggCGCGCGCTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTcc                                      <  1:1745297/101‑1 (MQ=21)
                                                          gcgcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCCATTCACCCCCCGGGAAAATTTGGGCCAAGAAAAAAAAGGGTTTTTct                                      >  2:4171903/1‑100 (MQ=9)
                                                          gcgcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                     <  1:4587556/101‑1 (MQ=21)
                                                          gcgcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                     <  1:2668849/101‑1 (MQ=21)
                                                           cgcgcggTCTGAGACG.AAACGGGGG.GCGCAGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                    <  1:4756132/101‑1 (MQ=21)
                                                           cgcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                    <  1:4215930/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                   <  2:214335/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                   <  2:2039042/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                  <  1:993367/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                  <  2:4004732/101‑1 (MQ=21)
                                                               cgcTCTGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                <  1:1003710/101‑1 (MQ=21)
                                                                  tctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                             >  2:1505667/1‑101 (MQ=21)
                                                                  tctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                             <  1:3873153/101‑1 (MQ=21)
                                                                  tctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                             <  1:2501970/101‑1 (MQ=21)
                                                                  tctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                             <  1:5314409/101‑1 (MQ=21)
                                                                  tctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                             <  1:2074366/101‑1 (MQ=21)
                                                                  tctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                             >  2:2496858/1‑101 (MQ=21)
                                                                  tcggagaCG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATGAAGAAGGTCTTTTGCTTACATCtt                             <  2:173324/98‑1 (MQ=17)
                                                                   ctGAGACG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                            <  1:5036225/101‑1 (MQ=21)
                                                                     gagaCG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTTGGGCAAGATAAAGAAGGGCTTTTTCTTTTTTCTTTcc                          >  2:4983817/1‑101 (MQ=17)
                                                                      agaCG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                         >  1:4487235/1‑101 (MQ=21)
                                                                       gaCG.AAACGGGGG.GCGCTGCCTCT.CTA.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                        <  2:5761941/101‑1 (MQ=21)
                                                                       gaCG.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CAAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                        <  2:3555961/101‑1 (MQ=21)
                                                                          g.AAACGGGGG.GCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                     <  1:1183627/101‑1 (MQ=21)
                                                                            aaacggggg.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                    >  1:802128/6‑101 (MQ=21)
                                                                            aaacggggg.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                    <  2:1746015/96‑1 (MQ=21)
                                                                              gcagggg.gCGCTG.CCCC.CTC.TT..GGCC.TTTGG.CCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg            >  2:5474240/4‑101 (MQ=0)
                                                                              gcagggg.gCGCTG.CCCC.CTC.TT..GGCC.TTTGG.CCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg            >  1:476340/4‑101 (MQ=0)
                                                                              acggggg.gCGCTGCCCCC.CTT.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTT.GGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                 <  1:3639029/98‑1 (MQ=21)
                                                                               cggggg.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                 <  1:3815916/99‑1 (MQ=25)
                                                                               cggggg.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                 <  2:3847742/99‑1 (MQ=25)
                                                                               cagggg.gCGCTG.CCCC.CTC.TT..GGCC.TTTGG.CCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg           >  1:728053/3‑101 (MQ=0)
                                                                                ggggg.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                >  2:335983/2‑101 (MQ=25)
                                                                                agggg.gCGCTGCCCCCTCTT.GG..C..C.TTTGG.CCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg          >  1:1969377/2‑101 (MQ=0)
                                                                                agggg.gCGCTG.CCCC.CTC.TT..GGCC.TTTGG.CCAATTCA.CCCCTGGG.ATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg          >  1:3552899/2‑101 (MQ=0)
                                                                                 gggg.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc               <  1:2242135/101‑1 (MQ=25)
                                                                                 gggg.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc               <  1:125581/101‑1 (MQ=25)
                                                                                 ggcg.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc               <  2:950081/98‑1 (MQ=25)
                                                                                  ggg.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAATAACCCg              >  1:1221029/1‑101 (MQ=25)
                                                                                  ggg.gCGATGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg              <  1:918241/101‑1 (MQ=25)
                                                                                    g.gCGCTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggg            <  1:5368097/101‑1 (MQ=37)
                                                                                       cgcTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggggg          <  2:1271088/101‑1 (MQ=37)
                                                                                        gcTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa         >  2:4478594/1‑101 (MQ=37)
                                                                                        gcTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa         <  2:3056052/101‑1 (MQ=37)
                                                                                        gcTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa         >  1:5471362/1‑101 (MQ=37)
                                                                                         cTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg        >  1:4357674/1‑101 (MQ=37)
                                                                                         cTGCCCCC.CTC.TT..GGCC.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg        >  1:553744/1‑101 (MQ=37)
                                                                                                           |                                                                                               
ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGG.GCGCTGCCCCC.CTC.CT..G..C.TTCGG.CCAATTCACCCCCCGGG.ATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCT  >  NC_014034/2566458‑2566648

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