Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 530,939 +T 98.6% intergenic (+17/+124) valS → / ← RCAP_rcc00497 valyl‑tRNA synthetase/UbiA prenyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034530,9391.→T98.6% 2.7 ‑1.1 ‑1.0289intergenic (+17/+124)valS/RCAP_rcc00497valyl‑tRNA synthetase/UbiA prenyltransferase
Reads supporting (aligned to +/- strand):  new base (121/164):  ref base (3/1):  total (124/165)
Fisher's exact test strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.20e-02

Only 100 of 342 total aligned reads displayed.

GAAGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTC.C.CGCGCTTG.TCGGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCAA  >  NC_014034/530839‑531022
                                                                                                     |                                                                                     
gAAGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTc                                                                                        <  2:4623167/101‑1 (MQ=255)
gAAGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTc                                                                                        <  2:8395112/101‑1 (MQ=255)
gAAGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTc                                                                                        <  1:4721731/101‑1 (MQ=255)
gAAGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTc                                                                                        >  1:4669333/1‑101 (MQ=255)
 aaccccGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTct                                                                                       <  1:8249890/98‑2 (MQ=255)
  aGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctg                                                                                      >  1:4021243/1‑99 (MQ=255)
    cccGCGAAAGCCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGCCTctgcc                                                                                    <  2:3406877/101‑5 (MQ=255)
      cgcgAAAGCCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgcccc                                                                                  <  1:7876249/101‑7 (MQ=255)
       gcgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccc                                                                                 <  2:5543408/101‑8 (MQ=255)
       gcgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccc                                                                                 >  2:1976268/1‑94 (MQ=255)
        cgAAACCCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgcccccc                                                                                <  2:2792505/101‑9 (MQ=255)
        cgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgcccccc                                                                                <  2:5152461/101‑9 (MQ=255)
        cgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgcccccc                                                                                >  1:7491471/1‑93 (MQ=255)
        cgAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgcccccc                                                                                >  1:691426/1‑93 (MQ=255)
         gAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccccg                                                                               >  1:4515172/1‑92 (MQ=255)
         gAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccccg                                                                               <  1:1616448/101‑10 (MQ=255)
         gAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccccg                                                                               >  1:1077537/1‑92 (MQ=255)
          aaaaCCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTctgccccccgc                                                                              <  1:3148398/101‑11 (MQ=255)
                         gcgcgGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAANTNGGATAAAAGAAGAAGGGNNNNNN.NTCCNNCCNC.NTNNCTCCCCCCggg                                                             >  1:8882499/1‑101 (MQ=255)
                             ggAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTC.T.GGCCCCCGCTTCGCTCCCCCCGGGatat                                                         >  1:2206821/1‑101 (MQ=25)
                              gAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATAtt                                                        >  1:3269975/1‑101 (MQ=25)
                              gAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATAtt                                                        >  2:1398932/1‑101 (MQ=25)
                               aaGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATAttt                                                       <  1:5577129/101‑1 (MQ=25)
                                agagGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGCAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATAtttt                                                      <  2:6267727/101‑1 (MQ=21)
                                agagGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATAtttt                                                      <  1:7655885/101‑1 (MQ=25)
                                 gagGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTc                                                     <  1:1398964/101‑1 (MQ=25)
                                    gcggcgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGc                                                  <  1:647149/101‑1 (MQ=25)
                                    gcggcgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGc                                                  <  1:5833702/101‑1 (MQ=25)
                                     cggcgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCa                                                 <  1:6673207/101‑1 (MQ=25)
                                     cggcgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCCGGGTATTTTCGGCa                                                 >  1:6023330/1‑101 (MQ=25)
                                       gcgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAg                                               <  1:3609556/101‑1 (MQ=25)
                                        cgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGa                                              <  2:7809475/101‑1 (MQ=21)
                                        cgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGa                                              <  2:5746347/101‑1 (MQ=21)
                                        cgAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGa                                              >  1:6782795/1‑101 (MQ=21)
                                         gAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGAt                                             >  2:5065474/1‑101 (MQ=21)
                                            gATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATaaa                                          <  2:7652059/101‑1 (MQ=21)
                                             aTCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAg                                         >  1:7027330/1‑101 (MQ=21)
                                               cggcgcggcgCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGc                                       <  1:1539834/101‑1 (MQ=21)
                                               cggcgcggcgCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGc                                       <  2:662194/101‑1 (MQ=21)
                                                 gcgcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGccc                                     <  1:2646731/101‑1 (MQ=21)
                                                 gcgcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGccc                                     <  2:6906819/101‑1 (MQ=21)
                                                 gcgcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGccc                                     <  1:2056866/101‑1 (MQ=21)
                                                  cgcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCg                                    <  2:2446119/101‑1 (MQ=21)
                                                   gcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCgg                                   <  2:4868728/101‑1 (MQ=21)
                                                   gcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCgg                                   >  1:294774/1‑101 (MQ=21)
                                                   gcggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTCGGGCAAGATgaaatcatgc                                   <  2:4798453/101‑8 (MQ=9)
                                                     ggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGaa                                 >  2:3263317/1‑101 (MQ=21)
                                                     ggcgcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGaa                                 <  2:4250330/101‑1 (MQ=21)
                                                        gcTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGcc                              <  1:3937712/101‑1 (MQ=21)
                                                          tGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTc                            >  2:2352215/1‑101 (MQ=21)
                                                          tGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTc                            >  2:7380892/1‑101 (MQ=21)
                                                          tGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTc                            <  2:4417193/101‑1 (MQ=21)
                                                           gAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTcc                           <  2:6134694/101‑1 (MQ=21)
                                                           gAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTcc                           >  1:1177749/1‑101 (MQ=21)
                                                           gAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTcc                           >  1:6342401/1‑101 (MQ=21)
                                                           gAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTcc                           >  1:864770/1‑101 (MQ=21)
                                                           gAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGATCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTcc                           <  1:5293209/101‑1 (MQ=21)
                                                             aGCTTCTGGCCGAACTCGGA.AAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCtt                        <  1:8040508/101‑1 (MQ=16)
                                                             aGCGGCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCt                         <  2:1563535/101‑1 (MQ=21)
                                                              gCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCtt                        <  1:2483536/101‑1 (MQ=21)
                                                               ggTCTGGCCGATCTCGGATAAAGGAAGAAGTGGGGCTCTG.CCCTCCGC.TTCGCTCCCCCCGGGATATTTTCGGTAAGATAAAGGCCCGGAAGCCTCCCttt                       <  2:3047519/100‑1 (MQ=16)
                                                               cGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCttt                       <  2:4021128/101‑1 (MQ=21)
                                                               cGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCttt                       <  2:5188448/101‑1 (MQ=21)
                                                                gTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCtttt                      >  2:1944297/1‑100 (MQ=21)
                                                                 tCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTAt                     <  2:4930445/101‑1 (MQ=21)
                                                                 tCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTAt                     <  2:8556057/101‑1 (MQ=21)
                                                                 tCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTAt                     <  1:4414863/101‑1 (MQ=21)
                                                                 tCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTAt                     <  2:3903593/101‑1 (MQ=21)
                                                                  cTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATc                    >  1:1882648/1‑101 (MQ=21)
                                                                    ggCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCtt                  <  1:2103404/101‑1 (MQ=21)
                                                                    ggCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCtt                  <  2:128916/101‑1 (MQ=21)
                                                                     gCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTg                 <  1:4143926/101‑1 (MQ=21)
                                                                     gCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTg                 >  1:8444539/1‑101 (MQ=21)
                                                                     gCCGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTg                 <  1:3700896/101‑1 (MQ=21)
                                                                      ccGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCCAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGc                >  2:8952976/1‑101 (MQ=21)
                                                                      ccGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGc                >  1:3468301/1‑101 (MQ=21)
                                                                      ccGAACTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGc                <  1:8335258/101‑1 (MQ=21)
                                                                         accTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCGTGCCgg             <  1:9026013/99‑1 (MQ=18)
                                                                         aaCTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCgg             >  2:8240703/1‑101 (MQ=21)
                                                                         aaCTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCgg             >  2:8031407/1‑101 (MQ=21)
                                                                          aCTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGa            <  2:8692262/101‑1 (MQ=21)
                                                                          aCTCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGa            <  1:7381112/101‑1 (MQ=21)
                                                                            tCGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGaaa          >  1:7134647/1‑101 (MQ=21)
                                                                             cGGATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAat         <  2:9173961/101‑1 (MQ=21)
                                                                               gATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGCAatat       >  2:217536/1‑101 (MQ=21)
                                                                               gATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAatat       <  2:5975767/101‑1 (MQ=21)
                                                                               gATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAatat       >  2:2233796/1‑101 (MQ=21)
                                                                               gATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAatat       <  2:7223023/101‑1 (MQ=21)
                                                                               gATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAatat       <  1:938413/101‑1 (MQ=21)
                                                                               gATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAatat       <  1:8716134/101‑1 (MQ=21)
                                                                               gATAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAatat       <  1:7977652/101‑1 (MQ=21)
                                                                                aTAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATc      <  2:1177731/101‑1 (MQ=21)
                                                                                aTAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATc      >  1:7772030/1‑101 (MQ=21)
                                                                                aTAAAGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATc      >  1:2261806/1‑101 (MQ=21)
                                                                                  aaaGGAAGAAGGGGGGCTCTG.CCCCCCTC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATccc    <  1:6185500/101‑1 (MQ=25)
                                                                                  aaaGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATccc    >  2:6346022/1‑101 (MQ=21)
                                                                                  aaaGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATccc    <  1:6541549/101‑1 (MQ=21)
                                                                                   aaGGAGGAAGGGGGGCTC.CGCCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCa   <  1:2788474/101‑1 (MQ=25)
                                                                                   aaGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCa   <  2:2093558/101‑1 (MQ=21)
                                                                                   aaGGAAGAAGGGGGGCTCTG.CCCCCCGC.TTCGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCa   >  2:8307616/1‑101 (MQ=21)
                                                                                                     |                                                                                     
GAAGCCCGCGAAAACCTCGCCGCGCGCGCGGAAGAGGCGGCGAAGATCGGCGCGGCGCTGAAGCGTCTGGCCGAACTCGGATAAAGGAAGAAGGGGGGCTC.C.CGCGCTTG.TCGGCTCCCCCCGGGATATTTTCGGCAAGATAAAGGCCCGGAAGCCTCCCTTTATCTTGCCGGAAATATCCCAA  >  NC_014034/530839‑531022

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG