Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,185,849 +A 2.6% intergenic (+122/‑383) mtlK → / → RCAP_rcc02028 mannitol 2‑dehydrogenase/FeoA family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,185,8490.→A 2.6% 2.2 ‑3.2 0.0113intergenic (+122/‑383)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (21/89):  total (24/89)
Fisher's exact test strand distribution p-value = 4.07e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

ACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCG.C.CCCGCAC.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCG.CCCCTTTCACCTTGG  >  NC_014034/2185751‑2185945
                                                                                                      |                                                                                                 
aCAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCAc.c                                                                                                  <  1:5048947/101‑1 (MQ=255)
 cAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCAc.cc                                                                                                 <  1:4573547/101‑1 (MQ=255)
  aGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCAc.ccc                                                                                                <  2:4078286/101‑1 (MQ=255)
  aGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCAc.ccc                                                                                                <  2:1221750/101‑1 (MQ=255)
   ggcggcCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCAc.cccc                                                                                               <  1:2621333/101‑1 (MQ=255)
   ggcgccCCCGCCCCCTCCGGGGCGGGCCCGCCTGATGTGCAAAGCGGTTTCGCCTGGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCAc.cccc                                                                                               <  2:3432327/101‑1 (MQ=25)
    gcggcCCCGCCCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCAc.ccccg                                                                                              <  2:4983117/101‑2 (MQ=255)
     cggcCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCAc.ccccgg                                                                                             <  1:3498628/101‑3 (MQ=255)
      ggcCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCAc.ccccggc                                                                                            <  2:3206208/101‑4 (MQ=255)
      ggcCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCCGGGGG.GCCG.C.CCCCCCc.cccccccc                                                                                           >  2:2093032/1‑98 (MQ=255)
          ccGGTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTAt                                                                                       <  2:4036619/101‑1 (MQ=37)
          ccGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTAt                                                                                       <  1:1841854/101‑1 (MQ=255)
          ccGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTAt                                                                                       <  2:3679457/101‑1 (MQ=255)
           cGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTTCGGGTGCGTCGC.C.CCCGCA..CCCCCGCCGTAtt                                                                                      >  2:2462156/1‑101 (MQ=37)
            gCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTAttt                                                                                     <  1:228923/101‑1 (MQ=255)
            gCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTAttt                                                                                     <  1:3517443/101‑1 (MQ=255)
            gCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTAttt                                                                                     <  1:95148/101‑1 (MQ=255)
            gCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTAttt                                                                                     <  1:730499/101‑1 (MQ=255)
               ccctcccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATc                                                                                  <  1:1868562/101‑1 (MQ=255)
               ccctcccGGGCGGGGCCCCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATc                                                                                  <  2:4537091/101‑1 (MQ=37)
                 ctcccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCG.CACCCGCA..CCCCCGGCGTATTTATccc                                                                                >  2:228241/1‑101 (MQ=255)
                  tcccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCa                                                                               <  1:4749860/101‑1 (MQ=255)
                  tcccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCa                                                                               <  1:1223768/101‑1 (MQ=255)
                    ccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAg                                                                             <  2:3021215/101‑1 (MQ=255)
                    ccGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAg                                                                             <  1:46104/101‑1 (MQ=255)
                     cgggcgggGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGa                                                                            <  1:459671/101‑1 (MQ=255)
                      gggcgggGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGaa                                                                           <  1:3739960/101‑1 (MQ=255)
                       ggcgggGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGaaa                                                                          <  1:3249757/101‑1 (MQ=255)
                         cgggGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGaaaaa                                                                        <  1:2119488/101‑1 (MQ=255)
                         cgggCCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGaaaaa                                                                        <  1:2126203/101‑1 (MQ=37)
                          ggggCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAag                                                                       <  2:2777606/101‑1 (MQ=255)
                           gggCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAagc                                                                      <  1:3400001/101‑1 (MQ=255)
                           gggCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAagc                                                                      >  1:2484105/1‑101 (MQ=255)
                           gggCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAagc                                                                      <  1:1898738/101‑1 (MQ=255)
                            ggCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAagca                                                                     >  2:4759680/1‑101 (MQ=255)
                            ggCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAagca                                                                     >  2:3292261/1‑101 (MQ=255)
                            ggCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAagca                                                                     <  1:3198208/101‑1 (MQ=255)
                                gccTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGGGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCCTATTTATCCCAAGAAAAAGCAGCgc                                                                 >  1:1682943/1‑100 (MQ=25)
                                  cTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGCCCCTGCGGGGGT.TCCG.C.CCCCGGA.CCCCCGGCGCATTTATCCCAATAAAAAGCAGCGGCt                                                               >  1:2176644/1‑101 (MQ=25)
                                  cTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCt                                                               <  2:1640060/101‑1 (MQ=255)
                                  cTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCt                                                               >  1:1795061/1‑101 (MQ=255)
                                  cTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.ACCG...CTCACACACCCCCGGCCTATTTATTCCAAGAAAAAGGGGCGGCt                                                               >  1:3301980/1‑101 (MQ=255)
                                   tGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCtt                                                              <  2:1642858/101‑1 (MQ=255)
                                   tGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCtt                                                              <  1:4489633/101‑1 (MQ=255)
                                   tGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCtt                                                              <  2:4332183/101‑1 (MQ=255)
                                    gATGTGGAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTg                                                             <  1:809972/101‑1 (MQ=37)
                                       gtgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGccc                                                          <  2:1936531/101‑1 (MQ=255)
                                       gtgAAAGGCGGTTTCGCCCCGGCGGTGCCAAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGccc                                                          <  1:4048547/101‑1 (MQ=25)
                                        tgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGcccc                                                         <  2:66189/101‑1 (MQ=255)
                                        tgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGcccc                                                         <  1:490190/101‑1 (MQ=255)
                                        tgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGcccc                                                         <  1:4059741/101‑1 (MQ=255)
                                        tgCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGcccc                                                         <  1:4463196/101‑1 (MQ=255)
                                         gCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCg                                                        <  2:280206/101‑1 (MQ=255)
                                          cAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCG.C.CCCCCCA.CCCCCGGCGGATTTATCCCAAGAAAAAGCAGCGGCTTGccccca                                                       >  1:419643/1‑99 (MQ=37)
                                            aaGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTTGCCCCAAg                                                     >  1:31765/1‑101 (MQ=25)
                                             aGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTTCCCCCaaag                                                    >  2:2196765/1‑99 (MQ=25)
                                                 gTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCtt                                                <  2:2980386/101‑1 (MQ=255)
                                                  tttCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTg                                               <  1:4773584/101‑1 (MQ=255)
                                                  tttCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTg                                               <  1:4643577/101‑1 (MQ=255)
                                                   ttCGCCTCGGCGGTGCCGAGGCGCCCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGt                                              <  1:2196651/101‑1 (MQ=37)
                                                     cGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCTAAGGCCTTGttt                                            <  2:2440291/101‑1 (MQ=37)
                                                     cGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGttt                                            <  1:1283345/101‑1 (MQ=255)
                                                      gCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGtttt                                           <  1:2871970/101‑1 (MQ=255)
                                                       ccTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGttttg                                          <  2:3075115/101‑2 (MQ=37)
                                                       ccTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGttttg                                          <  1:3419891/101‑2 (MQ=37)
                                                        cTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCCCAAGCCCTGttttgt                                         >  2:4533735/1‑99 (MQ=21)
                                                         tCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCG.C.CCCCCCA.CCCCCCCCGAATTTACCCCAAGAAAAATCAGCGGGTTGCCCCGGAGGCCTTGtttttat                                        >  2:3850357/1‑99 (MQ=25)
                                                            gcggTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGc                                    <  1:3777973/101‑1 (MQ=37)
                                                            gcggTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGc                                    <  2:581207/101‑1 (MQ=37)
                                                              ggTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGCcac                                  <  2:4277961/101‑1 (MQ=25)
                                                               gTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGcaca                                 <  2:2665562/101‑1 (MQ=37)
                                                               gTGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGcaca                                 <  1:900188/101‑1 (MQ=37)
                                                                tGCCGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACaa                                <  2:3187750/101‑1 (MQ=37)
                                                                 gCCGAGGCGCCCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGTAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAg                               <  1:2084202/101‑1 (MQ=25)
                                                                   cGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGaa                             >  2:2302790/1‑101 (MQ=37)
                                                                   cGAGGCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGaa                             <  2:1144203/101‑1 (MQ=37)
                                                                      ggCGACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACg                          <  1:1311180/101‑1 (MQ=37)
                                                                         gACCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCg                       <  1:1649405/101‑1 (MQ=37)
                                                                         gACCCTGCGGGGGA.TACGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGGAAAAGCAGCGGCTTGCCCCGCAGGCCATG.TTTTAATGGGACCAGAAACGccc                       >  1:847232/1‑100 (MQ=16)
                                                                          aCCCTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGa                      <  1:3169649/101‑1 (MQ=37)
                                                                            ccTGCGGGGGC.TACGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAAc                    <  2:3038815/101‑1 (MQ=37)
                                                                             cTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAAcc                   <  1:2620239/101‑1 (MQ=37)
                                                                             cTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAAcc                   <  2:1368301/101‑1 (MQ=37)
                                                                             cTGCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAAcc                   <  1:4683533/101‑1 (MQ=37)
                                                                               gCGGGGGC.TCCGCC.CCCGCAC.CCCCC.GCGTATTTATCCCAAGAAAAAGCAGCAGCTTGACCCCCAAGCCCTG.TTTTGCTGGCACAAAGAAAGCGGAAACG.c                >  2:3049471/1‑101 (MQ=9)
                                                                               gCGGGGGC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.c                <  2:618519/101‑1 (MQ=37)
                                                                                  ggggC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTTCCCCCAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.cccc             >  2:1962068/1‑101 (MQ=25)
                                                                                  ggggC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAAAGCCGAACCG.cccc             >  2:1354543/1‑101 (MQ=37)
                                                                                   gggC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.ccccc            <  2:393602/101‑2 (MQ=37)
                                                                                   gggC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.ccccc            <  2:1677436/101‑2 (MQ=37)
                                                                                    ggC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAAAGCCGAACCG.ccccct           >  2:2714577/1‑99 (MQ=37)
                                                                                     gC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCctt          <  2:4530431/101‑1 (MQ=37)
                                                                                     gC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCctt          <  1:1364705/101‑1 (MQ=37)
                                                                                     gC.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCctt          <  1:4525240/101‑1 (MQ=37)
                                                                                      c.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCcttt         >  2:3924689/1‑100 (MQ=37)
                                                                                      c.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCcttt         <  2:4879063/101‑2 (MQ=37)
                                                                                      c.TCCGCC.CCCGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCcttt         <  1:3259981/101‑2 (MQ=37)
                                                                                        tccgcc.cccGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCCctttc        >  1:4713053/6‑99 (MQ=37)
                                                                                          cgcc.cccGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCAc       <  2:1283000/98‑1 (MQ=38)
                                                                                               cccGCA..CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCCCTTTCACCTTg   <  1:3678633/101‑1 (MQ=38)
                                                                                                      |                                                                                                 
ACAGGCGGCCCCGCTCCCTCCCGGGCGGGGCCGCCTGATGTGCAAAGCGGTTTCGCCTCGGCGGTGCCGAGGCGACCCTGCGGGGGC.TCCG.C.CCCGCAC.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCG.CCCCTTTCACCTTGG  >  NC_014034/2185751‑2185945

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG