Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,042,953 G→T 100% W217L (TGG→TTG)  thiD → phosphomethylpyrimidine kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,042,9530G→T100.0% 617.8 164W217L (TGG→TTG) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  new base (88/75):  ref base (0/0):  total (89/75)

Only 100 of 170 total aligned reads displayed.

CGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTGGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTTCG  >  NC_014034/3042855‑3043050
                                                                                                  |                                                                                                 
cGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTtgg                                                                                                 <  1:3503402/101‑1 (MQ=255)
cGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTtgg                                                                                                 >  2:2871762/1‑101 (MQ=255)
  aTCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCg                                                                                               >  2:929460/1‑101 (MQ=255)
  aTCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCg                                                                                               <  1:2966381/101‑1 (MQ=255)
  aTCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCg                                                                                               >  1:2827921/1‑101 (MQ=255)
   tCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGa                                                                                              >  2:4360843/1‑101 (MQ=255)
    ctgctgCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGAc                                                                                             <  2:2806907/101‑1 (MQ=255)
    ctgctgCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGAc                                                                                             <  2:3611328/101‑1 (MQ=255)
      gctgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGc                                                                                           >  1:3200862/1‑101 (MQ=255)
      gccgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGc                                                                                           <  2:1966649/98‑1 (MQ=255)
       ctgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCa                                                                                          >  1:534369/1‑101 (MQ=255)
        tgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAg                                                                                         >  2:356954/1‑101 (MQ=255)
        tgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAg                                                                                         <  2:4162978/101‑1 (MQ=255)
         gctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGc                                                                                        >  1:3301720/1‑101 (MQ=255)
          ctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCt                                                                                       <  2:2546645/101‑1 (MQ=255)
             gATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGc                                                                                    >  1:1178138/1‑101 (MQ=255)
              aTGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGcc                                                                                   <  2:1535428/101‑1 (MQ=255)
                    cgGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATc                                                                             <  2:688102/101‑1 (MQ=255)
                     gggCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATcc                                                                            >  2:148052/1‑101 (MQ=255)
                     gggCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATcc                                                                            <  1:1442842/101‑1 (MQ=255)
                      ggCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATccc                                                                           <  1:2619987/101‑1 (MQ=255)
                        ggAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCg                                                                         <  2:4183303/100‑1 (MQ=255)
                         gAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGa                                                                        >  2:2348665/1‑101 (MQ=255)
                         gAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGa                                                                        >  2:2629167/1‑101 (MQ=255)
                         gAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGa                                                                        >  2:4064212/1‑101 (MQ=255)
                          aGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGAt                                                                       >  2:1281767/1‑101 (MQ=255)
                           ggTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATc                                                                      <  1:1095785/101‑1 (MQ=255)
                            gTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCt                                                                     >  2:785598/1‑101 (MQ=255)
                               tGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTgcg                                                                  <  2:4209675/101‑1 (MQ=255)
                               tGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTgcg                                                                  <  1:1405886/101‑1 (MQ=255)
                                ggTCCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGCCACCGGCTGCTGCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTgcgc                                                                 <  2:4274759/101‑1 (MQ=255)
                                 gTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGcc                                                                <  1:252008/101‑1 (MQ=255)
                                 gTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGcc                                                                <  1:4246835/101‑1 (MQ=255)
                                  ttCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCa                                                               <  2:4224755/101‑1 (MQ=255)
                                  ttCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCa                                                               <  2:412974/101‑1 (MQ=255)
                                  ttCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCa                                                               >  1:2370442/1‑101 (MQ=255)
                                     gAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgc                                                            <  1:3029746/101‑1 (MQ=255)
                                     gAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgc                                                            >  2:1886178/1‑101 (MQ=255)
                                     gAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgc                                                            >  1:2561104/1‑101 (MQ=255)
                                      aGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTTGCCCGCATCCCCGATCTGCGCCAGgcc                                                           >  1:3104808/1‑101 (MQ=255)
                                         cgccgcGGGTCGAGACGAAAAACCCCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgcc                                                        <  1:1325446/101‑1 (MQ=255)
                                             gcgGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccg                                                    <  1:3239923/101‑1 (MQ=255)
                                             gcgGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccg                                                    <  2:2513150/101‑1 (MQ=255)
                                              cgGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgc                                                   >  1:3406946/1‑101 (MQ=255)
                                               gggTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgcc                                                  <  1:3851600/101‑1 (MQ=255)
                                                ggTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGggcgccgccgccg                                                 >  1:2700531/1‑101 (MQ=255)
                                                ggTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgccg                                                 >  1:4137954/1‑101 (MQ=255)
                                                  tCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgccgcc                                               >  2:397173/1‑101 (MQ=255)
                                                    gagaCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgcccg                                             >  1:1363823/1‑101 (MQ=255)
                                                     agaCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgcccgc                                            >  2:3190687/1‑101 (MQ=255)
                                                     agaCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgcccgc                                            >  1:3482116/1‑101 (MQ=255)
                                                      gaCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCcgcg                                           >  2:1829990/1‑101 (MQ=255)
                                                         gAAAAACCCCCAGGGCACCGGCTGCGCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCt                                        <  2:1505696/101‑1 (MQ=255)
                                                         gAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCt                                        >  1:2856989/1‑101 (MQ=255)
                                                         gAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCt                                        >  2:2788665/1‑101 (MQ=255)
                                                           aaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCtat                                      <  2:2264170/101‑1 (MQ=255)
                                                             aaCCCCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATAc                                    <  1:3131374/101‑1 (MQ=255)
                                                             aaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATAc                                    >  2:3488978/1‑101 (MQ=255)
                                                              acccccATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACg                                   <  1:2420871/98‑1 (MQ=255)
                                                              acacCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACg                                   >  1:210600/1‑101 (MQ=255)
                                                                acCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCt                                 >  1:3900228/1‑101 (MQ=255)
                                                                  ccATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTgc                               >  1:209515/1‑101 (MQ=255)
                                                                        cACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCa                         <  2:2375022/101‑1 (MQ=255)
                                                                         cccGGCTGCCCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCACGCAGCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCAt                        <  2:4417055/100‑1 (MQ=255)
                                                                         aCCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCAt                        >  2:187712/1‑101 (MQ=255)
                                                                         aCCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCAt                        >  1:1233985/1‑101 (MQ=255)
                                                                          cgggCGGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATc                       <  2:1014589/99‑1 (MQ=255)
                                                                           cGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcg                      <  1:4419556/101‑1 (MQ=255)
                                                                           cGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcg                      <  1:16281/101‑1 (MQ=255)
                                                                             gCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgcc                    >  2:1871611/1‑101 (MQ=255)
                                                                              cTGCCCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                   <  2:1401848/101‑1 (MQ=255)
                                                                              cTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                   >  1:876790/1‑101 (MQ=255)
                                                                              cTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                   >  2:1436442/1‑101 (MQ=255)
                                                                              cTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                   <  1:2377127/101‑1 (MQ=255)
                                                                              cTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                   >  2:3219661/1‑101 (MQ=255)
                                                                              cTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                   <  2:3289412/101‑1 (MQ=255)
                                                                               tGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccgc                  >  2:2197106/1‑101 (MQ=255)
                                                                               tGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccgc                  <  2:3601452/101‑1 (MQ=255)
                                                                                gCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCcgcg                 <  1:3577513/101‑1 (MQ=255)
                                                                                gCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCcgcg                 >  1:393465/1‑101 (MQ=255)
                                                                                gCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCcgcg                 <  1:3149282/101‑1 (MQ=255)
                                                                                gCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCcgcg                 <  1:2890756/101‑1 (MQ=255)
                                                                                 cACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCgg                >  2:696705/1‑101 (MQ=255)
                                                                                  aCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGc               >  2:627484/1‑101 (MQ=255)
                                                                                  aCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCACCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGACATCGCCGCGGc               >  2:1162949/1‑101 (MQ=255)
                                                                                   cccTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCAGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGcc              <  1:785680/101‑1 (MQ=255)
                                                                                     cTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGa            >  1:1853770/1‑101 (MQ=255)
                                                                                     cTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGCCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGa            <  1:683625/101‑1 (MQ=255)
                                                                                      tGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGAc           <  2:1693171/101‑1 (MQ=255)
                                                                                       gTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGAcc          >  1:430034/1‑101 (MQ=255)
                                                                                        tCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCg         >  1:3534187/1‑101 (MQ=255)
                                                                                          ttCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGc       <  2:3986587/101‑1 (MQ=255)
                                                                                          ttCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGc       >  1:1437090/1‑101 (MQ=255)
                                                                                           tCGGCCTTGGCGACGCGGCCGGCCCGCATCCCCGATCTGCGCCAGGCCGCCCCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCt      <  2:4254045/101‑1 (MQ=255)
                                                                                           tCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCt      >  1:2166779/1‑101 (MQ=255)
                                                                                            cGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCtt     >  1:3465894/1‑101 (MQ=255)
                                                                                            cGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCtt     >  1:592836/1‑101 (MQ=255)
                                                                                              gCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTtc   <  2:3817124/101‑1 (MQ=255)
                                                                                              gCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTtc   >  2:2932188/1‑101 (MQ=255)
                                                                                               ccTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTTCGTCCCACCGCCGCGCCCGACCCGCTTtcg  >  1:2982266/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTGGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTTCG  >  NC_014034/3042855‑3043050

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG