Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,566,557 +G 100% 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,5571.→G100.0% 506.1 184intergenic (+29/+165)purD/sdaAphosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand):  new base (33/148):  ref base (0/1):  total (33/151)

Only 100 of 222 total aligned reads displayed.

ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GG.GGGCGCTGCCCCCCTCCT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566458‑2566644
                                                                                                        |                                                                                        
aTGTGATGGTCGATGCGATCGACTGGCGCGAGGGGTTCTGCCTCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..g                                                                                           <  2:4231738/101‑1 (MQ=255)
aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..g                                                                                           >  2:3418250/1‑101 (MQ=255)
aTGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGACGCGACATCGGCTGGCGCGAGCGCTGAGACGAAACGGG.GGGCGCCGCCCCCCCCCT..g                                                                                           >  1:4319592/1‑101 (MQ=255)
   cGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..ggcc                                                                                        <  1:146102/101‑4 (MQ=39)
     ggggTTGATGCGATCGACTGGCCCGGGGTGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GGAGGGCGCTGCCCCCCTCTT..gttctt                                                                                      <  2:162798/99‑6 (MQ=255)
     aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGGGGCCCCCCCCCT..gggcct                                                                                      >  1:3057443/1‑96 (MQ=255)
     aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..ggcctt                                                                                      <  2:2803819/101‑6 (MQ=38)
      tGGTCGAGCCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..ggccttc                                                                                     <  2:4084918/101‑7 (MQ=37)
          cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGcc                                                                                 <  1:3815222/101‑1 (MQ=37)
          cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGcc                                                                                 <  1:936316/101‑1 (MQ=37)
            aTGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCaa                                                                               <  1:2214065/101‑1 (MQ=37)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTc                                                                            <  1:408834/101‑1 (MQ=37)
               cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTc                                                                            <  1:1152670/101‑1 (MQ=37)
               cGAGCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTc                                                                            <  1:1697672/101‑1 (MQ=25)
                gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCa                                                                           <  2:1860831/101‑1 (MQ=37)
                  tCGACTGGCCCGGGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAcc                                                                         <  1:2288985/101‑1 (MQ=25)
                    gACTGTCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAcccc                                                                       <  2:1103212/101‑1 (MQ=25)
                    gACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGGTGCCCCCCTCCTTGG..CCTCGGCCCATTTAcccc                                                                       >  2:195032/1‑101 (MQ=255)
                    gACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTTGGG..CCTCGGGCACTTCAcccc                                                                       >  2:546213/1‑101 (MQ=255)
                      cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAcccccc                                                                     <  2:3157648/101‑1 (MQ=25)
                      cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAcccccc                                                                     <  1:2895076/101‑1 (MQ=25)
                       tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCg                                                                    <  2:2580723/101‑1 (MQ=25)
                         gCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGTCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCggg                                                                  <  1:4309561/101‑1 (MQ=25)
                              aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATAtt                                                             <  1:15296/101‑1 (MQ=21)
                               ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATAttt                                                            <  1:4823843/101‑1 (MQ=21)
                               ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATAttt                                                            <  1:3399827/101‑1 (MQ=21)
                                 ggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCCCTT..GGCCTTCGGCCCATTCACCCCCCCGGAGATTTct                                                          >  2:201523/1‑100 (MQ=21)
                                  gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGGGCTGCCCCCCTCTT..GGCCTTCGGCCCATTTCCCCCCCGGGGGATTTCgg                                                         >  1:2398914/1‑101 (MQ=21)
                                  gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCCAATCACCCCCCGGGATATTTTgg                                                         >  2:4745212/1‑101 (MQ=21)
                                   ttCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCggg                                                        >  2:2644884/1‑101 (MQ=21)
                                    tCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGc                                                       <  2:2196847/101‑1 (MQ=21)
                                     cTGCCGCCGCGACTTCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCa                                                      <  2:213462/101‑1 (MQ=21)
                                     cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCa                                                      <  1:1928170/101‑1 (MQ=21)
                                       gcctccGCGACATCGGCTGGCGCGCGTTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  2:558382/101‑1 (MQ=18)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  2:404170/101‑1 (MQ=21)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  1:4004963/101‑1 (MQ=21)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  1:1811994/101‑1 (MQ=21)
                                       gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  2:612616/101‑1 (MQ=21)
                                        ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGa                                                   <  1:3738925/101‑1 (MQ=21)
                                         cgccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGAt                                                  <  1:3134645/101‑1 (MQ=21)
                                          gccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  1:758907/101‑1 (MQ=21)
                                           ccgcGACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGGGGCGCCCCCCTCCT..GGCCTTTCGCCCATTCACCCCCCGGGATACTTCGGGCGAGACaa                                                >  2:812208/1‑101 (MQ=16)
                                             gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  2:641627/101‑1 (MQ=21)
                                             gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:3206906/101‑1 (MQ=21)
                                             gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:2556790/101‑1 (MQ=21)
                                             gcgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:3001494/101‑1 (MQ=21)
                                              cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGCTATATAGGGGAAGATaaaaa                                             >  2:3690992/1‑99 (MQ=18)
                                              cgACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                             <  2:2991831/101‑1 (MQ=21)
                                               gACATCGGCTGGCGCGCGCTCTGGGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaa                                            <  2:1840456/101‑1 (MQ=21)
                                               gACATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCAACCCCCGGGATATTTCGGGCAAGATAaagaa                                            <  1:1375187/101‑1 (MQ=21)
                                                aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTTGGGCAAGATAaaaaag                                           >  1:1738075/1‑101 (MQ=18)
                                                aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaag                                           <  1:3942432/101‑1 (MQ=21)
                                                aCATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTACGGGCAAGATAaagaag                                           <  2:2525570/101‑1 (MQ=21)
                                                 cATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                          <  2:3278843/101‑1 (MQ=21)
                                                 cATCGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                          <  1:48548/101‑1 (MQ=21)
                                                    cGGCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCt                                       >  1:598450/1‑101 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  2:1275649/101‑1 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      >  2:4666595/1‑101 (MQ=21)
                                                     ggCTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  2:337048/101‑1 (MQ=21)
                                                     ggCTGGCGCGCCCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  2:4794865/101‑1 (MQ=21)
                                                       cTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  2:174030/101‑1 (MQ=21)
                                                       cTGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  1:1041996/101‑1 (MQ=21)
                                                        tGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   <  2:3994153/101‑1 (MQ=21)
                                                        tGGCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   <  2:5020865/101‑1 (MQ=21)
                                                         ggCGCGCGCTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTcc                                  >  1:4956338/1‑101 (MQ=21)
                                                          gcgcgcgcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCt                                 >  1:4505293/1‑101 (MQ=21)
                                                           cgcgcgcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:2093234/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  1:3285597/101‑1 (MQ=21)
                                                            gcgcgcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  1:3186209/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  2:3302022/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  2:2870611/101‑1 (MQ=21)
                                                             cgcgcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  2:192013/101‑1 (MQ=21)
                                                              gcgcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  1:2952744/101‑1 (MQ=21)
                                                              gcgcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCa                             <  2:426189/101‑1 (MQ=21)
                                                                gcTCTGAGACGAACCGGG.GGGCGCTGCCCCCTTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:2652893/101‑1 (MQ=21)
                                                                gcTCTGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATc                           <  1:2154159/101‑1 (MQ=21)
                                                                 ctctGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                          >  2:4678419/1‑101 (MQ=21)
                                                                  tctGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                         <  1:88760/101‑1 (MQ=21)
                                                                   ctGAGGCGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  1:2566418/101‑1 (MQ=21)
                                                                   ctGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  2:4898629/101‑1 (MQ=21)
                                                                    tGAGACGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGc                       <  1:4249700/101‑1 (MQ=21)
                                                                     gagaCGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCCGGAGATTTTGGGGCAGATAAAAAAGGGGTTTTTCTTTTATTttttc                      >  1:3572030/1‑98 (MQ=2)
                                                                       gaCGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                    <  2:3407080/101‑1 (MQ=21)
                                                                       gaCGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                    <  1:514012/101‑1 (MQ=21)
                                                                       gaCGAAACGGG.GGGCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                    <  1:1894315/101‑1 (MQ=21)
                                                                           aaacggg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                <  2:107585/96‑1 (MQ=21)
                                                                           aaacggg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGAAATTCCGGCCAAAAAAAGAAAGGACTTTTCTTTTCTTCTTCCCCGaaaa                >  2:4359411/6‑100 (MQ=2)
                                                                            aacggg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTATTTTCCTTTCATCTTGCCCGAAata               >  1:557794/5‑101 (MQ=21)
                                                                             acggg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat              <  2:3198269/98‑1 (MQ=25)
                                                                             acggg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat              <  2:908820/98‑1 (MQ=25)
                                                                             acggg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat              <  1:2129400/98‑1 (MQ=25)
                                                                              cggg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             <  2:1023372/99‑1 (MQ=25)
                                                                              cggg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             <  2:237318/99‑1 (MQ=25)
                                                                              cggg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc             >  2:3038455/3‑101 (MQ=25)
                                                                                gg.gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc           <  2:2453882/101‑1 (MQ=25)
                                                                                gg.gggCGCTG.CCCCCTCTT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggg     >  2:1480212/1‑100 (MQ=0)
                                                                                 g.gggCCCTGCCCCCTTTTG..GC.CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggga    <  2:359618/101‑3 (MQ=2)
                                                                                   gggCGCTGCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgg         <  2:2353969/101‑1 (MQ=37)
                                                                                          gCCCCCCTCTT..GGCCTTCGGCCAATTCTCCCCCCGGGATATTTCGGTCAAGATAAAGAAGGTCTTTTCCTTTCATCTTTCCCGAAATATCCCGGGGGAGcc  <  1:709533/101‑1 (MQ=25)
                                                                                          gCCCCCCTCTT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGGCCGAAATATCCCGGGGGAGcc  >  2:4973558/1‑101 (MQ=25)
                                                                                                        |                                                                                        
ATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACGAAAC.GG.GGGCGCTGCCCCCCTCCT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566458‑2566644

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