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% 1293.5 339W217L (TGG→TTG) thiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  new base (182/157):  ref base (0/0):  total (182/157)

Only 100 of 351 total aligned reads displayed.

CGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTGGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTTCGG  >  NC_014034/3042855‑3043051
                                                                                                  |                                                                                                  
cGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTtgg                                                                                                  <  2:542696/101‑1 (MQ=255)
cGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTtgg                                                                                                  <  1:4455002/101‑1 (MQ=255)
cGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTAGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTtgg                                                                                                  >  2:259868/1‑101 (MQ=255)
  ctCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCCTCGGGGGTCGAGACGAAAAACACCCATGGCACCGCTTGCAACCTGACTTCGGCCTTGGCg                                                                                                >  2:6401989/2‑101 (MQ=255)
  aTCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCg                                                                                                <  2:6598144/101‑1 (MQ=255)
    ctgctgCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCCCCGGCTGCACCCTGTCTTCGGCCTTGGCGAc                                                                                              <  2:1237476/101‑1 (MQ=255)
    ctgctgCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGAc                                                                                              >  1:1111306/1‑101 (MQ=255)
    ctgctgCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGAc                                                                                              >  1:5167844/1‑101 (MQ=255)
     tgctgcTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACg                                                                                             <  1:2608629/101‑1 (MQ=255)
      gctgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGc                                                                                            <  1:3247661/101‑1 (MQ=255)
       ctgctCTATGCGGCGGGCCAGGTCTGGTTCGAGGCGCTGCGAGTTGAGACGAAAGACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCa                                                                                           <  2:5801706/101‑1 (MQ=255)
       ctgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCa                                                                                           <  2:3882266/101‑1 (MQ=255)
        tgctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAg                                                                                          >  1:4169178/1‑101 (MQ=255)
         gctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGc                                                                                         >  2:686076/1‑101 (MQ=255)
         gctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGc                                                                                         <  2:2867494/101‑1 (MQ=255)
         gctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGc                                                                                         >  1:6222634/1‑101 (MQ=255)
         gctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGc                                                                                         >  2:3538202/1‑101 (MQ=255)
          ctCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCt                                                                                        >  1:1619898/1‑101 (MQ=255)
          cgcGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCt                                                                                        >  1:4253784/3‑101 (MQ=255)
           tCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTg                                                                                       <  1:4722089/101‑1 (MQ=255)
            cGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTgg                                                                                      <  1:4064208/101‑1 (MQ=255)
            cGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTgg                                                                                      <  1:1748722/101‑1 (MQ=255)
             gATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGc                                                                                     >  2:5932543/1‑101 (MQ=255)
             gATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGc                                                                                     >  2:6760845/1‑101 (MQ=255)
             gATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGc                                                                                     >  1:1868102/1‑101 (MQ=255)
                     gggCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATcc                                                                             >  1:2362370/1‑101 (MQ=255)
                     gggCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATcc                                                                             >  2:1165755/1‑101 (MQ=255)
                      ggCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATccc                                                                            <  1:5659354/101‑1 (MQ=255)
                      ggCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATccc                                                                            >  2:1924119/1‑101 (MQ=255)
                       gCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGTCCCggaaccac                                                                           >  2:4524367/1‑94 (MQ=255)
                        cGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCCTCCCCg                                                                          >  2:1676318/1‑101 (MQ=255)
                         gAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGa                                                                         >  2:5962842/1‑101 (MQ=255)
                          aGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGAt                                                                        >  1:1998675/1‑101 (MQ=255)
                            gTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCt                                                                      >  1:6738513/1‑101 (MQ=255)
                             tCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTg                                                                     >  1:2525536/1‑101 (MQ=255)
                              cTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTgc                                                                    <  1:11695/101‑1 (MQ=255)
                                 gTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGcc                                                                 <  2:1398989/101‑1 (MQ=255)
                                   tCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAg                                                               <  1:1501571/101‑1 (MQ=255)
                                   tCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAg                                                               >  1:6476820/1‑101 (MQ=255)
                                   tCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAg                                                               >  1:4112935/1‑101 (MQ=255)
                                      gggCGCCGCGGATCGAGACGAAAAACTCCCATGGCCCCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCTAGgcc                                                            <  2:2507795/100‑1 (MQ=255)
                                       ggCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccg                                                           >  1:5328489/1‑101 (MQ=255)
                                       ggCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccg                                                           >  1:1292938/1‑101 (MQ=255)
                                          gccgcGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccg                                                        <  2:4168983/101‑1 (MQ=255)
                                              cgGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgc                                                    >  1:4577286/1‑101 (MQ=255)
                                              cgGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgc                                                    >  2:3567458/1‑101 (MQ=255)
                                               gggTCGAGACGAAAAACACCCATGGCACCGNCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgcc                                                   <  1:1018894/101‑1 (MQ=255)
                                               gggTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgcc                                                   <  2:1075367/101‑1 (MQ=255)
                                                ggTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCGTCCCCGATCTGCGCCAGgccgccgccgccg                                                  >  2:775821/1‑101 (MQ=255)
                                                 gTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCATgccgcagccgccgc                                                 >  1:5831434/1‑101 (MQ=255)
                                                 gTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgccgc                                                 >  1:6167105/1‑101 (MQ=255)
                                                  tCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGgccgccgccgccgcc                                                >  1:5440111/1‑101 (MQ=255)
                                                   cGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgccc                                               >  2:2068358/1‑101 (MQ=255)
                                                   cGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgccc                                               >  2:2393903/1‑101 (MQ=255)
                                                   cGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgccc                                               <  2:6677474/101‑1 (MQ=255)
                                                     agaCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgcccgc                                             >  2:148999/1‑101 (MQ=255)
                                                     agaCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgcccgc                                             >  1:6535484/1‑101 (MQ=255)
                                                     agaCGAAAAACACCCATGGAACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGccgcccgc                                             >  2:3379621/1‑101 (MQ=255)
                                                       aCGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCgcgc                                           >  2:2163899/1‑101 (MQ=255)
                                                        cGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGcc                                          <  2:928644/101‑1 (MQ=255)
                                                         gAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCt                                         <  1:3485123/101‑1 (MQ=255)
                                                         gAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCt                                         <  2:3767223/101‑1 (MQ=255)
                                                          aaaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCTATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCta                                        >  1:7089458/1‑101 (MQ=255)
                                                          aaaaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCtc                                        >  1:3320825/1‑100 (MQ=255)
                                                            aaaCACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCtata                                      <  2:507843/101‑1 (MQ=255)
                                                              acacCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACg                                    >  2:612275/1‑101 (MQ=255)
                                                              acacCCATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACg                                    >  1:4951867/1‑101 (MQ=255)
                                                                 cccATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTg                                 <  2:3028832/101‑1 (MQ=255)
                                                                  ccATGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTgc                                >  1:92729/1‑101 (MQ=255)
                                                                    aTGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATCCGCTgcgc                              >  1:2066697/1‑101 (MQ=255)
                                                                    aTGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTgcgc                              >  1:6887592/1‑101 (MQ=255)
                                                                     tGGCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTgcgcg                             >  2:4370780/1‑101 (MQ=255)
                                                                       gCACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGcc                           >  1:1556615/1‑101 (MQ=255)
                                                                        cACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCa                          <  2:85106/101‑1 (MQ=255)
                                                                        cACCGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCa                          <  1:3766120/101‑1 (MQ=255)
                                                                         cccGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCAt                         <  2:5343535/100‑1 (MQ=255)
                                                                         cccGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCAt                         <  2:1593693/100‑1 (MQ=255)
                                                                         cccGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCAt                         <  1:3592160/100‑1 (MQ=255)
                                                                           cGGCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcg                       <  1:5972160/101‑1 (MQ=255)
                                                                            ggCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgc                      <  2:1795303/101‑1 (MQ=255)
                                                                            ggCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgc                      >  1:3941/1‑101 (MQ=255)
                                                                             gCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgcc                     >  1:3894983/1‑101 (MQ=255)
                                                                             gCTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgcc                     >  1:5872068/1‑101 (MQ=255)
                                                                              cTGCCCCCTGTCTTCGGCCTTGGCGACGCCGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                    <  1:2966172/101‑1 (MQ=255)
                                                                              cTGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccg                    <  1:229842/101‑1 (MQ=255)
                                                                               tGCGCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccgc                   <  1:143470/101‑1 (MQ=255)
                                                                               tGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccgc                   >  1:5037853/1‑101 (MQ=255)
                                                                               tGCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATcgccgc                   >  1:3093941/1‑101 (MQ=255)
                                                                                gCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCcgcg                  >  2:3088758/1‑101 (MQ=255)
                                                                                gCACCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCcgcg                  >  1:287469/1‑101 (MQ=255)
                                                                                  aCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGc                >  1:591432/1‑101 (MQ=255)
                                                                                  aCCCTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGc                >  1:5541993/1‑101 (MQ=255)
                                                                                   cccTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCTCGCGCCTATACGCTGCGCGCCATCGCCGCGGcc               <  1:6372096/101‑1 (MQ=255)
                                                                                     cTGTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGa             >  1:242108/1‑101 (MQ=255)
                                                                                       gTCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGAcc           >  2:4043596/1‑101 (MQ=255)
                                                                                        tCTTCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCg          >  1:756033/1‑101 (MQ=255)
                                                                                          ttCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCgtcgtca        <  2:5167584/101‑7 (MQ=255)
                                                                                           tCGGCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCt       <  2:4812119/101‑1 (MQ=255)
                                                                                             ggCCTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCttt     <  1:3416841/101‑1 (MQ=255)
                                                                                               ccTTGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTtcg   >  2:5651779/1‑101 (MQ=255)
                                                                                                  |                                                                                                  
CGATCTGCTGCTCGATGCCGCGGGCGAGGTCTGGTTCGAGGCGCCGCGGGTCGAGACGAAAAACACCCATGGCACCGGCTGCACCCTGTCTTCGGCCTGGGCGACGCAGCTGGCCCGCATCCCCGATCTGCGCCAGGCCGCCGCCGCCGCCCGCGCCTATACGCTGCGCGCCATCGCCGCGGCCGACCGGCTTTCGG  >  NC_014034/3042855‑3043051

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