Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,185,975 G→C 7.4% intergenic (+248/‑257) 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,9750G→C 7.4% 2.9 ‑0.5 ‑0.467intergenic (+248/‑257)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (11/51):  total (11/56)
Fisher's exact test strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.97e-01

AAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCG.CCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..CCC..G.CCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAG  >  NC_014034/2185875‑2186068
                                                                                                          |                                                                                                  
aaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..g                                                                                                    >  2:882043/1‑101 (MQ=255)
aaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCG.CCCTCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..g                                                                                                    <  2:6157796/101‑1 (MQ=255)
 agcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..c                                                                                                 <  1:3349982/101‑1 (MQ=255)
 agcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGGCCGGGGGCGG..A..G..c                                                                                                 >  2:4708284/1‑101 (MQ=38)
 agcagcGGCTTGCCCCGAAGCCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..c                                                                                                 <  1:7217022/101‑1 (MQ=37)
  tcaggtgCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cc                                                                                                <  1:2980912/95‑1 (MQ=37)
   cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccc                                                                                               >  2:8215895/1‑101 (MQ=255)
   cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccc                                                                                               >  1:8901137/1‑101 (MQ=255)
   cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccc                                                                                               >  1:8211260/1‑101 (MQ=255)
   cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccc                                                                                               >  1:767092/1‑101 (MQ=255)
    cgcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAGCGCCGAACCGCCCC.CTTTCACCTTGGCTCAGATATCCCGGGGGTCCGGGGGCGG..A..G..cccc                                                                                              <  2:7536884/100‑2 (MQ=37)
    agcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccc                                                                                              <  2:753010/101‑2 (MQ=37)
    agcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccc                                                                                              >  1:5106648/1‑100 (MQ=37)
    agcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGGCCGGGGGCGG..A..G..cccc                                                                                              >  2:2354475/1‑100 (MQ=37)
    agcGGCTTGCCCCGAAGGCCTTG.TATTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccc                                                                                              >  2:4374812/1‑100 (MQ=37)
     gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  2:2869756/101‑3 (MQ=37)
     gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  1:7532949/101‑3 (MQ=37)
     gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  1:587127/101‑3 (MQ=37)
     gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  2:5437737/101‑3 (MQ=37)
     gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  2:6178695/101‑3 (MQ=37)
     gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  1:1773115/101‑3 (MQ=37)
      cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccg                                                                                            <  1:8844289/101‑4 (MQ=37)
      cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccg                                                                                            <  2:1656953/101‑4 (MQ=37)
      cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccg                                                                                            <  2:472920/101‑4 (MQ=37)
      cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccg                                                                                            <  2:8389981/101‑4 (MQ=37)
       ggCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccgg                                                                                           <  1:8015490/101‑5 (MQ=25)
        gCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggc                                                                                          >  2:6016404/1‑96 (MQ=25)
          ttGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCGCCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggcgg                                                                                        >  2:1629877/1‑94 (MQ=25)
           tGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggcggg                                                                                       >  2:1883161/1‑93 (MQ=21)
            gCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggcgggg                                                                                      >  1:554025/1‑92 (MQ=21)
            gCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggcgggg                                                                                      >  2:2098289/1‑92 (MQ=21)
            gCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggcgggg                                                                                      >  2:1882937/1‑92 (MQ=21)
             ccccGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggcgggtc                                                                                     <  1:7356805/101‑11 (MQ=21)
             ccccGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggcgggtc                                                                                     <  2:6479185/101‑11 (MQ=21)
                                          gAACCGCCGACCCCCCCC.CTTTCCCCTGGGCTAAAATCTCCGGGGGGTCCGGGGGCGG..ACCC..CCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAacac                                                      <  1:5525199/101‑1 (MQ=0)
                                            aaCGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACAccc                                                    <  1:1118319/101‑1 (MQ=2)
                                             aCGCCGACCCGCCCC.CTTTCGCCTTGGCTCAAATATCCGGGGGGTCGGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACAcccc                                                   <  1:1116316/101‑1 (MQ=2)
                                             aCGCCGACCCGCCCC.CTTTCCCCTGGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACAcccc                                                   <  2:1511626/101‑1 (MQ=2)
                                             aCGCCGACCCCCCCC.CTTTCCCCTTGGCTCAAATATCCGGGGGGTCGGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACAcccc                                                   <  1:2773655/101‑1 (MQ=2)
                                              cGCCGAACCGCCCC.CTTTCACCTGGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCa                                                  <  2:4988315/101‑1 (MQ=2)
                                                  gAACCGCCCC.CTTCCACCTTGGCTCAAATATCCCGGGGCTCCGGGGCCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAgcgc                                              <  2:2687490/101‑1 (MQ=2)
                                                   acccgcccc.cTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGcc                                             <  1:6808976/95‑1 (MQ=0)
                                                   acccccccc.cTTTCCCCTTGGCTCAAATACCCGGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGcc                                             <  1:6858652/95‑1 (MQ=0)
                                                     ccgcccc.cTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCac                                           <  2:6501987/97‑1 (MQ=2)
                                                      cgcccc.cTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCacc                                          <  1:2895984/98‑1 (MQ=2)
                                                      cccccc.cTTCCACCTTGGCTCAAATACCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCacc                                          <  2:8665780/98‑1 (MQ=2)
                                                        cccc.cTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCaccga                                        <  2:1703336/100‑1 (MQ=9)
                                                         ccc.cTTCCCCCTTGGCTCAAATACCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCaccgac                                       <  1:8501997/101‑1 (MQ=2)
                                                          ccccgctccccTTGTCTCGAATCCCCCGGTGCTCCGGGGGCGC..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCaccgac                                       <  2:7349847/92‑1 (MQ=2)
                                                             cTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCAccgaccga                                    <  2:4852213/101‑1 (MQ=14)
                                                              ttgcaccggggCTCAAATCTCCGGGGGGTCGGGGGGCGG..A..G..CCCCGG.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACcgaccgac                                   <  1:2710784/92‑1 (MQ=2)
                                                                tccccTGGGCTCAAATACCCGGGGGGCCGGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACtt                                 <  1:871499/98‑1 (MQ=2)
                                                                tCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACtt                                 <  2:2158383/101‑1 (MQ=14)
                                                                tCACCTGGCCACAAATACCCGGGGGGCCGGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACtt                                 <  2:519103/101‑1 (MQ=2)
                                                                 cACCTGGGCTAAAATACCCCGGGGGTCCGGGGGCGG..C..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTg                                <  2:4622402/101‑1 (MQ=9)
                                                                  cccTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGc                               <  2:5545285/100‑1 (MQ=14)
                                                                  cccTTGGCTCAAATACCCGGGGGGTCGGGGGGCGG..AC.C..CCC..GTGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGc                               <  1:423312/100‑1 (MQ=9)
                                                                   ccTGGGCTCAAATCGCCGGGGGGTCGGGGGGCGG..A..C..CCCCGG.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCg                              <  1:8796984/101‑1 (MQ=9)
                                                                   ccTGGGCTCAAATACCCGGGGGGCCGGGGGGCGG..A..C..CCCCGG.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCg                              <  1:2354542/101‑1 (MQ=9)
                                                                     ttGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGAc                            <  1:8669472/101‑1 (MQ=14)
                                                                     ttGGCTAAAATACCCCGGGGGCCGGGGGGCGGCCC..C..CCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGAc                            <  1:780946/101‑1 (MQ=9)
                                                                      tGGCTCAAATACCCCGGGGGTCCGGGGGCGG..ACCC..CCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACg                           <  1:1615606/101‑1 (MQ=9)
                                                                         cactaaTACCCCGGGGGTCCGGTGTCGC..ACCC..CCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGcc                        <  1:2585898/97‑1 (MQ=9)
                                                                         cTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..CCC..C.CGGCGGGGGGCCCAGGGGCAAGCCCGCGGCGAAACACCCCCCCGCCCCCGCCCCACCCGCGCCGGcc                        >  2:2888435/1‑101 (MQ=2)
                                                                         cTCAAATACCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGcc                        <  2:5337728/101‑1 (MQ=9)
                                                                             aatcgcccGGGGGCCGGGGGGCGG..G..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTc                    <  1:1925722/96‑1 (MQ=9)
                                                                             aaaacccggggggTCCGGGGGGGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTc                    <  2:767079/93‑1 (MQ=9)
                                                                             aaTACCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTc                    <  2:5698457/101‑1 (MQ=12)
                                                                              atatCCCGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTcg                   <  2:5877156/101‑1 (MQ=17)
                                                                              atagCCCGGGGGTCGGGGGGCGG..G..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTcg                   <  1:9071844/101‑1 (MQ=12)
                                                                               taccccGGGGGTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTcgc                  <  2:8807057/98‑1 (MQ=12)
                                                                                  cccGGGGGTCCGGGGGCGG..ACCC..CCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCa               <  2:9070155/101‑1 (MQ=16)
                                                                                    ggggggCCGGGGGGCGG..C..G..CCCCGG.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCaaa             <  1:4788684/100‑1 (MQ=12)
                                                                                     gggggTCGGGGGGCGG..G..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCaaaa            <  1:7088895/101‑1 (MQ=12)
                                                                                     gggggTCGGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCaaaa            <  2:5106479/101‑1 (MQ=12)
                                                                                     gggggCCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCaaaa            <  2:8900731/101‑1 (MQ=12)
                                                                                     gggcgTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCaaaa            <  2:8759791/101‑1 (MQ=12)
                                                                                      ggggTCCGGGGGGGG..A..G..CCC..C.CGGCGGGTCGTCTCGGGCAAAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCAGTCGCGCAAAAg           <  1:6591085/101‑1 (MQ=12)
                                                                                      ggggTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAg           <  2:7728329/101‑1 (MQ=18)
                                                                                      ggggTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAg           <  2:2747861/101‑1 (MQ=18)
                                                                                      ggggTCCGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAg           <  2:9057280/101‑1 (MQ=18)
                                                                                       gggTCCGGGGGCGG..A..G..CCCCGG.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGc          <  2:3394380/101‑1 (MQ=18)
                                                                                        ggcccGGGGGCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGcc         <  2:379800/98‑1 (MQ=12)
                                                                                            cGGGGGCGG..C..G..CCCCGG.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGc     <  2:671107/101‑1 (MQ=16)
                                                                                             gggggCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGcc    <  2:7767353/101‑1 (MQ=18)
                                                                                             gggggCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGcc    <  2:8210911/101‑1 (MQ=18)
                                                                                             gggggCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGcc    <  2:554015/101‑1 (MQ=18)
                                                                                               gggCGG..A..GCCCCC..G.GCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCcag  <  2:7447771/101‑1 (MQ=18)
                                                                                                          |                                                                                                  
AAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCG.CCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..CCC..G.CCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAG  >  NC_014034/2185875‑2186068

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG