Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,185,974 +C 4.5% intergenic (+247/‑258) 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,9742.→C 4.5% 2.6 0.0 ‑0.567intergenic (+247/‑258)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (11/53):  total (11/56)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.08e-01

AAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCG.CCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..CCC..G.CCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAG  >  NC_014034/2185874‑2186068
                                                                                                         |                                                                                                    
aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..a                                                                                                       >  1:226891/1‑101 (MQ=255)
aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..a                                                                                                       >  2:6590825/1‑101 (MQ=255)
aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..a                                                                                                       >  1:6215721/1‑101 (MQ=255)
aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGATCAAATATCCCGGGGGTCCGGGGGCGG..a                                                                                                       >  1:6403777/1‑101 (MQ=38)
 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                                                                                               >  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)
    cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccc                                                                                               >  2:8215895/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                                                                                             <  1:1773115/101‑3 (MQ=37)
      gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  2:6178695/101‑3 (MQ=37)
      gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  1:7532949/101‑3 (MQ=37)
      gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  2:5437737/101‑3 (MQ=37)
      gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  2:2869756/101‑3 (MQ=37)
      gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..ccccc                                                                                             <  1:587127/101‑3 (MQ=37)
       cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccg                                                                                            <  2:472920/101‑4 (MQ=37)
       cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccg                                                                                            <  2:1656953/101‑4 (MQ=37)
       cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccg                                                                                            <  2:8389981/101‑4 (MQ=37)
       cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccg                                                                                            <  1:8844289/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                                                                                      >  2:1882937/1‑92 (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)
              ccccGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggcgggtc                                                                                     <  2:6479185/101‑11 (MQ=21)
              ccccGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCGCCCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..cccccggcgggtc                                                                                     <  1:7356805/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)
                                                                                                         |                                                                                                    
AAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAACCG.CCC.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGG..A..G..CCC..G.CCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCACCGACCGACTTGCGACGGCCTGTCGCGCAAAAGCCCGGCCAG  >  NC_014034/2185874‑2186068

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