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% 883.6 233W217L (TGG→TTG) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  new base (143/90):  ref base (0/0):  total (143/90)

Only 100 of 240 total aligned reads displayed.

CGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTGGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTTCGGT  >  NC_014034/3042855‑3043052
                                                                                                  |                                                                                                   
cGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTtgg                                                                                                   >  2:4327246/1‑101 (MQ=255)
  aTCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCg                                                                                                 <  2:3774629/101‑1 (MQ=255)
    ctgctgCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGAc                                                                                               >  2:123939/1‑101 (MQ=255)
    ctgctgCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGAc                                                                                               >  2:231297/1‑101 (MQ=255)
      gctgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGc                                                                                             >  2:1728795/1‑101 (MQ=255)
       ctgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCa                                                                                            >  1:2487164/1‑101 (MQ=255)
          ctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCt                                                                                         <  1:3588315/101‑1 (MQ=255)
             gATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGTGACGCAGCTGGc                                                                                      >  2:3873885/1‑101 (MQ=255)
             gATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGc                                                                                      >  1:3224811/1‑101 (MQ=255)
              aTGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGcc                                                                                     >  2:119056/1‑101 (MQ=255)
                  cgcgGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCa                                                                                 >  2:911179/1‑101 (MQ=255)
                   gcgGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCAt                                                                                <  2:2156141/101‑1 (MQ=255)
                    cgGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATc                                                                               >  1:2040247/1‑101 (MQ=255)
                    cgGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATc                                                                               <  1:1649804/101‑1 (MQ=255)
                      ggcgggtcacggTCCGAGGGGGGGCGGGTCGAGGCGAAAAACACCCGGGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATccc                                                                             <  2:15595/91‑1 (MQ=255)
                        cGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCg                                                                           >  1:2129664/1‑101 (MQ=255)
                        cGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTCGGCGACGCAGCTGGCCCGCATccctg                                                                           >  2:4714229/1‑99 (MQ=255)
                          aGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGAt                                                                         >  1:703039/1‑101 (MQ=255)
                           ggTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATc                                                                        >  1:590528/1‑101 (MQ=255)
                            gTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCt                                                                       >  1:2790610/1‑101 (MQ=255)
                             tCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTg                                                                      <  2:1222862/101‑1 (MQ=255)
                              cTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTgc                                                                     >  2:1433792/1‑101 (MQ=255)
                               tGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTgcc                                                                    >  1:1217172/1‑100 (MQ=255)
                                ggTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTgcgc                                                                   <  2:3129750/101‑1 (MQ=255)
                                ggTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTgcgc                                                                   >  2:1291016/1‑101 (MQ=255)
                                 gTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGcc                                                                  >  1:4672717/1‑101 (MQ=255)
                                 gTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGcc                                                                  >  1:1876032/1‑101 (MQ=255)
                                 gTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGcc                                                                  >  1:964361/1‑101 (MQ=255)
                                      aGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgcc                                                             <  2:4210546/101‑1 (MQ=255)
                                      aGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgcc                                                             <  2:2376731/101‑1 (MQ=255)
                                          gccgcGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGTTGGCCCGCATCCCCGATCTGCGCCAGgccgccg                                                         >  1:4123274/1‑101 (MQ=255)
                                               gggTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgcc                                                    <  2:3644292/101‑1 (MQ=255)
                                                ggTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgccg                                                   >  1:2760023/1‑101 (MQ=255)
                                                  tCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgccgcc                                                 <  2:51493/101‑1 (MQ=255)
                                                  tCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgccgcc                                                 <  1:2406867/101‑1 (MQ=255)
                                                   cGAGACGAAAAACACCCATGGCACCGGCTGCCCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgccc                                                <  1:1852278/101‑1 (MQ=255)
                                                   cGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgccc                                                >  1:812779/1‑101 (MQ=255)
                                                   cGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgccc                                                >  1:4473203/1‑101 (MQ=255)
                                                   cGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgccc                                                >  1:2296122/1‑101 (MQ=255)
                                                   cGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgccc                                                >  1:2972214/1‑101 (MQ=255)
                                                    gagaCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgcccg                                               >  2:2288375/1‑101 (MQ=255)
                                                     agaCGAAAACCACCCATCGCACGGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCGTCCCCGATCTGCGCCAGGCCGCCGCCGccgcccgc                                              <  2:754182/101‑1 (MQ=255)
                                                     agaCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgcccgc                                              >  2:1028570/1‑101 (MQ=255)
                                                       aCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCgcgc                                            >  2:3244247/1‑101 (MQ=255)
                                                        cGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGcc                                           <  2:3541439/101‑1 (MQ=255)
                                                          aaaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCta                                         <  2:3449471/101‑1 (MQ=255)
                                                          aaaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCta                                         >  1:3510310/1‑101 (MQ=255)
                                                          aaaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCta                                         >  2:459102/1‑101 (MQ=255)
                                                          aaaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCta                                         >  1:4685584/1‑101 (MQ=255)
                                                           aaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCTCCtat                                        <  1:3488470/101‑1 (MQ=255)
                                                           aaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCtat                                        <  2:4465407/101‑1 (MQ=255)
                                                           aaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCtat                                        >  1:2252925/1‑101 (MQ=255)
                                                           aaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCtat                                        >  2:889165/1‑101 (MQ=255)
                                                             aaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATAc                                      <  2:1113717/101‑1 (MQ=255)
                                                             aaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATAc                                      <  1:4447300/101‑1 (MQ=255)
                                                             aaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATAc                                      <  1:4372623/101‑1 (MQ=255)
                                                             aaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATAc                                      <  1:1927169/101‑1 (MQ=255)
                                                              acacCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACg                                     >  2:2863262/1‑101 (MQ=255)
                                                                 cccATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTg                                  >  1:660297/1‑101 (MQ=255)
                                                                 cccATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTg                                  >  1:1398473/1‑101 (MQ=255)
                                                                    aTGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTgcgc                               >  2:2098375/1‑101 (MQ=255)
                                                                      ggCACCGGCTCCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTgcgcgc                             <  2:546639/101‑1 (MQ=255)
                                                                       gCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGcc                            >  1:4023231/1‑101 (MQ=255)
                                                                        cACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCa                           <  2:244658/101‑1 (MQ=255)
                                                                          ccGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATc                         <  2:1065480/101‑1 (MQ=255)
                                                                          ccGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATc                         <  1:980337/101‑1 (MQ=255)
                                                                          ccGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATc                         <  2:702975/101‑1 (MQ=255)
                                                                           cGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcg                        >  1:3116196/1‑101 (MQ=255)
                                                                            ggCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgc                       <  2:2487021/101‑1 (MQ=255)
                                                                            ggCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgc                       >  1:1005915/1‑101 (MQ=255)
                                                                              cTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                     <  1:2421457/101‑1 (MQ=255)
                                                                              cTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                     >  1:764417/1‑101 (MQ=255)
                                                                               tGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccgc                    >  2:1973083/1‑101 (MQ=255)
                                                                               tGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccgc                    >  1:1153728/1‑101 (MQ=255)
                                                                               tGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccgc                    <  1:3467232/101‑1 (MQ=255)
                                                                               tGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTACACGCTGCGCGCCATcgccgc                    >  2:2906112/1‑101 (MQ=255)
                                                                                gCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCcgcg                   >  2:2034382/1‑101 (MQ=255)
                                                                                gCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCcgcg                   <  1:504056/101‑1 (MQ=255)
                                                                                  aCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGc                 >  1:3700999/1‑101 (MQ=255)
                                                                                   cccTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGcc                >  2:481456/1‑101 (MQ=255)
                                                                                   cccTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGcc                <  1:3673126/101‑1 (MQ=255)
                                                                                   cccTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCCCTAGGCCGCCCCAGCCGCACGCGCCAATACGCTGCAGGCCATCGCCGCGGcc                >  2:603089/1‑101 (MQ=255)
                                                                                     cTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGa              >  1:3268966/1‑101 (MQ=255)
                                                                                     cTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGa              <  1:3441981/101‑1 (MQ=255)
                                                                                      tGTCTTCGGCCTTGGTGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCTGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGAc             <  1:3874127/101‑1 (MQ=255)
                                                                                       gTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGAcc            <  1:123955/101‑1 (MQ=255)
                                                                                       gTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGAcc            >  1:1551334/1‑101 (MQ=255)
                                                                                        tCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGGGCCATCGCCGCGGCCGACCg           >  2:346287/1‑101 (MQ=255)
                                                                                        tCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCg           <  1:1433886/101‑1 (MQ=255)
                                                                                         cTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCgg          >  1:3013621/1‑101 (MQ=255)
                                                                                            cGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCtt       <  1:3809422/101‑1 (MQ=255)
                                                                                             ggCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCttt      >  1:106550/1‑101 (MQ=255)
                                                                                             ggCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCttt      >  1:3835909/1‑101 (MQ=255)
                                                                                             ggCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCttt      <  2:2220933/101‑1 (MQ=255)
                                                                                              gCCTTGGCGCCGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTtc     <  2:1398379/101‑1 (MQ=255)
                                                                                              gCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTtc     >  2:4111499/1‑101 (MQ=255)
                                                                                              gCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTtc     <  1:103240/101‑1 (MQ=255)
                                                                                              gCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTtc     >  2:2802028/1‑101 (MQ=255)
                                                                                              gCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTtc     >  1:2007816/1‑101 (MQ=255)
                                                                                                 tGGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTtcggt  <  2:1670144/101‑1 (MQ=255)
                                                                                                  |                                                                                                   
CGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTGGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTTCGGT  >  NC_014034/3042855‑3043052

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