Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,394,542 2 bp→CC 100% coding (603‑604/708 nt) modC2 → molybdate ABC transporter ATP‑binding protein ModC

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,394,5420G→C100.0% 287.7 76E201D (GAG→GACmodC2molybdate ABC transporter ATP‑binding protein ModC
Reads supporting (aligned to +/- strand):  new base (19/57):  ref base (0/0):  total (19/57)
*NC_0140342,394,5430A→C100.0% 289.5 78T202P (ACC→CCC) modC2molybdate ABC transporter ATP‑binding protein ModC
Reads supporting (aligned to +/- strand):  new base (19/59):  ref base (0/0):  total (19/59)

TCGGCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGAGACCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATCCCTGCGCCTC  >  NC_014034/2394447‑2394635
                                                                                                ||                                                                                            
tCGGCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCccgc                                                                                          <  2:1243281/101‑1 (MQ=255)
 cGGCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCccgcc                                                                                         <  1:4635036/101‑1 (MQ=255)
 cGGCCCTTGATCCCGACCCCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCccgcc                                                                                         <  1:1850135/101‑1 (MQ=255)
  ggCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCccgccg                                                                                        <  1:436904/101‑1 (MQ=255)
    cgcTTGATCCCGACACCCGGCCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAg                                                                                      <  2:510243/101‑1 (MQ=255)
     gcTTGATCCCGCCCCCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAgg                                                                                     <  1:4878675/101‑1 (MQ=255)
        tGATCCCGACCCCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgc                                                                                  <  1:3903305/101‑1 (MQ=255)
        tGATCCCGACACCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgc                                                                                  <  2:3309728/101‑1 (MQ=255)
        tGACCCCGACCCCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgc                                                                                  <  2:963434/101‑1 (MQ=255)
         gATCCCGACACCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgcg                                                                                 >  1:1557750/1‑101 (MQ=255)
            cccGACACCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgcggcg                                                                              >  1:3543858/1‑101 (MQ=255)
             ccGACACCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAggcgcggcgc                                                                             >  1:5032957/1‑101 (MQ=255)
                  ccccGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGc                                                                        <  1:518530/100‑1 (MQ=255)
                     cGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGa                                                                     <  2:4285688/101‑1 (MQ=255)
                      gggCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGAc                                                                    >  2:1502228/1‑101 (MQ=255)
                       ggCCCGGCTGCAAGATCGC.CTGCCGCGCCGGCACGATCTCTTTGGCGCCACGCCGCGGAGGACAAGCGACGACCCCGCCGAGGCGCGGCGCCTGGCCGAcc                                                                   <  2:1815822/101‑1 (MQ=255)
                       ggCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGAcc                                                                   <  1:3503699/101‑1 (MQ=255)
                       ggCCCGGATGCAAGA.CAGCCGGCGCCGCCTGCACGATCGGTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGAcc                                                                   <  2:2118839/101‑1 (MQ=255)
                        gCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACcg                                                                  <  2:4679051/101‑1 (MQ=255)
                        gCCCGGCTGCAAGAGCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACcg                                                                  <  2:2041890/101‑1 (MQ=255)
                         cccGGCGGCAAGATCAG.CGGCGCCGCCTGCCCGATCTCTGGGGGGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACcgc                                                                 <  1:2014529/101‑1 (MQ=255)
                            ggCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtg                                                              <  1:2394768/101‑1 (MQ=255)
                            ggCCGCAAGGTCAG.CGGCCCCGCCTGCCCGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtg                                                              <  2:1229102/101‑1 (MQ=255)
                             gCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgc                                                             >  1:720806/1‑101 (MQ=255)
                             gCTGCAAGATCAG.CGGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgc                                                             <  1:3490009/101‑1 (MQ=255)
                              ctgaaggATCCG.CTGCGCCGCCTCCACGACCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgct                                                            <  1:1021463/95‑1 (MQ=255)
                              ctccAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgct                                                            <  2:1261514/98‑1 (MQ=255)
                                gCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgctgc                                                          <  2:3278639/101‑1 (MQ=255)
                                gCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgctgc                                                          >  2:3976461/1‑101 (MQ=255)
                                gCAAGAGCAG.CTGCGCCGCCTGCACGATCTTTTTGGCGCCACGACGCTGATAACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGtgctgc                                                          <  1:2464215/101‑1 (MQ=255)
                                 cAAGATCAG.CTGCGCCGCCTGCCCGAGCTCTTTGGCGCCAGGACGCTGATGAGAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCg                                                         <  2:1640172/101‑1 (MQ=255)
                                 cAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCg                                                         >  2:4169047/1‑101 (MQ=255)
                                 cAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCg                                                         <  2:3543588/101‑1 (MQ=255)
                                 cAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCg                                                         >  2:2312017/1‑101 (MQ=255)
                                  ttccaatctcTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCg                                                         <  2:4918637/92‑1 (MQ=255)
                                     aTCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTg                                                     >  1:1695778/1‑101 (MQ=255)
                                     aTCAG.CTGCCCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTg                                                     <  1:3265262/101‑1 (MQ=255)
                                      tCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTgg                                                    <  2:4055960/101‑1 (MQ=255)
                                      tCAG.CTGCGCCCCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTgg                                                    <  2:3685573/101‑1 (MQ=255)
                                      accg.CTCCGCCGCCTGCACGTTCTCTTTGGCGCCACGACTTTGTTGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTgg                                                    <  2:578872/98‑1 (MQ=255)
                                       cAG.CTGCCCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGa                                                   <  1:126878/101‑1 (MQ=255)
                                        aG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAg                                                  <  1:2418609/101‑1 (MQ=255)
                                         accTGCCCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAg                                                  <  2:1557620/99‑1 (MQ=255)
                                           cTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCgagatcggaa                                                >  1:4918964/1‑92 (MQ=255)
                                            tgcccgaccTGCACGGTTTTGTTGGCGCGACGAGGCGGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGcc                                               <  2:3751415/94‑1 (MQ=255)
                                              cgccgcCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCgg                                             >  1:583061/1‑101 (MQ=255)
                                               gccgccTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCggg                                            >  1:369630/1‑101 (MQ=255)
                                                ccgccTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGggc                                           <  1:3645452/101‑1 (MQ=255)
                                                  gccTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGggcgg                                         <  1:238757/101‑1 (MQ=255)
                                                  gccTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGggcgg                                         >  2:1344209/1‑101 (MQ=255)
                                                  gccCGCACTATTTCTTTGGCGCCGCGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGggcgg                                         <  2:3180395/101‑1 (MQ=255)
                                                   ccTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGgggggc                                        >  2:1488356/1‑101 (MQ=255)
                                                   ccTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGggcggc                                        <  1:158853/101‑1 (MQ=255)
                                                    cggcccGACCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCt                                       <  2:4867999/96‑1 (MQ=255)
                                                    cTGCCCGATCTCTTGGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCt                                       <  2:1894979/101‑1 (MQ=255)
                                                    cTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCt                                       <  1:4336533/101‑1 (MQ=255)
                                                     tGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGgggggccg                                      >  1:4891036/1‑99 (MQ=255)
                                                      gCACGCTCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTgg                                     <  1:2000396/101‑1 (MQ=255)
                                                      gCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTgg                                     >  2:1719590/1‑101 (MQ=255)
                                                       cACGGTCTCTTGGGCGCCAGCCCCCGGACGACAACCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTggg                                    <  2:2612000/101‑1 (MQ=255)
                                                       cACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTggg                                    <  2:4444837/101‑1 (MQ=255)
                                                              tctTTGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgggccgcg                             >  1:2279753/1‑101 (MQ=255)
                                                                tttGGCGCCACGACGCTGATGACAAGCCACGCCCCCGTCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgccgcgcc                           <  1:4312964/101‑1 (MQ=255)
                                                                tttGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgccgcgcc                           <  1:4169364/101‑1 (MQ=255)
                                                                tttGGCGCCACGACGCTGAGGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGgcgccgcgcc                           <  2:720717/101‑1 (MQ=255)
                                                                 ttGGCGCCACGAGGCGTATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGcgccgcgccg                          <  1:4071179/101‑1 (MQ=255)
                                                                 ttGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGcgccgcgccg                          <  1:1864692/101‑1 (MQ=255)
                                                                 ttGGCGCCACGACGCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGcgccgcgccg                          <  2:3251521/101‑1 (MQ=255)
                                                                   ggCGCCAGGCCGCAGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgcc                        <  2:3991188/101‑1 (MQ=255)
                                                                   ggCGCCACGACGCGGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgcc                        <  2:2760340/101‑1 (MQ=255)
                                                                   ggCCCCACGCCGCGGATGAAAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGACGcgccgcc                        <  1:1976785/101‑1 (MQ=255)
                                                                      gcCACGCCCCTGATGACAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgccgcc                     <  2:5032627/101‑1 (MQ=255)
                                                                      cccACGCCGTGGAGGACAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgccgcc                     <  2:369577/100‑1 (MQ=255)
                                                                         ccgccgCTGATGACAAGCCCCGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGAGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgccgccgcc                  <  1:1056159/97‑1 (MQ=255)
                                                                         ccgccgCTGAGAACAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGcgccgccgccgcc                  <  1:497471/97‑1 (MQ=255)
                                                                            ccgCGGATGAAAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCgccgccgccgccgcc               <  2:1703665/100‑1 (MQ=255)
                                                                            acgTTGATACCAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCgccgccgccgccgcc               <  2:1695649/101‑1 (MQ=255)
                                                                                 gatgaCAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATccc          <  1:2436358/101‑1 (MQ=255)
                                                                                 gatgaCAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGGGGCTGGGCGCCGCGCCGCCGCCGCCGCCCTccc          >  1:1062413/1‑101 (MQ=255)
                                                                                 gatgaAAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATccc          <  1:1488486/101‑1 (MQ=255)
                                                                                   tgaCAAGCCACGCCCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATCCCTg        <  2:4890709/101‑1 (MQ=255)
                                                                                    gaCAAGCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGGCGCCGCCCTCCCTgc       >  1:4845002/1‑101 (MQ=255)
                                                                                         gCCACGACCCCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATCCCTGCGCCTc  <  2:582982/101‑1 (MQ=255)
                                                                                                ||                                                                                            
TCGGCGCTTGATCCCGACACCCGGGCCCGGCTGCAAGATCAG.CTGCGCCGCCTGCACGATCTCTTTGGCGCCACGACGCTGATGACAAGCCACGAGACCGCCGAGGCGCGGCGCCTGGCCGACCGCGTGCTGCGGCTGGAGGCCGGGCGGCTGGGCGCCGCGCCGCCGCCGCCGCCATCCCTGCGCCTC  >  NC_014034/2394447‑2394635

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG