Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,185,930 +C 84.5% intergenic (+203/‑302) 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,9301.→C84.5% 30.1 ‑0.3 ‑1.571intergenic (+203/‑302)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (24/36):  ref base (5/6):  total (29/42)
Fisher's exact test strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.28e-02

GCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..CCCGCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGC  >  NC_014034/2185830‑2186024
                                                                                                       |                                                                                                  
gCGGGGGCTCCGCCCCCGCACCCCCC.GCGTATTTATCCCAAGAAAAAGCAGCAGCTTGACCCCCAAGCCCTG.TTTTGCTGGCACAAAGAAAGCGGAA.ACG.c                                                                                                   >  2:3049471/1‑101 (MQ=9)
gCGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.c                                                                                                   <  2:618519/101‑1 (MQ=37)
 cGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cc                                                                                                  >  2:3062775/1‑101 (MQ=37)
 cGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCCAA.CCG.cc                                                                                                  >  2:2986263/1‑101 (MQ=25)
  gggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cc.c                                                                                                <  2:4044468/101‑1 (MQ=37)
   ggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTTCCCCCAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cc.c.c                                                                                              >  2:1962068/1‑101 (MQ=25)
   ggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cc.c.c                                                                                              <  1:626703/101‑1 (MQ=37)
   ggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAAAGCCGAA.CCG.cc.c.c                                                                                              >  2:1354543/1‑101 (MQ=37)
    gggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cc.c.cc                                                                                             <  2:393602/101‑2 (MQ=37)
    gggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cc.c.cc                                                                                             <  2:1677436/101‑2 (MQ=37)
     ggcgccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cc.c.cct                                                                                            <  1:183810/93‑3 (MQ=37)
     ggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAAAGCCGAA.CCG.cc.c.cct                                                                                            >  2:2714577/1‑99 (MQ=37)
      gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.Cctt                                                                                           >  2:1619715/1‑101 (MQ=37)
      gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.Cctt                                                                                           <  1:4525240/101‑1 (MQ=37)
      gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.Cctt                                                                                           <  1:1364705/101‑1 (MQ=37)
      gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.Cctt                                                                                           <  2:4530431/101‑1 (MQ=37)
       cTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.Ccttt                                                                                          >  2:3924689/1‑100 (MQ=37)
       cTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.Ccttt                                                                                          <  2:4879063/101‑2 (MQ=37)
       cTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.Ccttt                                                                                          <  1:3259981/101‑2 (MQ=37)
        tccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.Cctttc                                                                                         >  1:4713053/6‑99 (MQ=37)
          cgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCAc                                                                                        <  2:1283000/98‑1 (MQ=38)
           cgcccccgcgcccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCAc                                                                                        <  2:511729/91‑1 (MQ=38)
             ccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCtt                                                                                     <  2:1166255/101‑1 (MQ=38)
              cccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTg                                                                                    <  1:3678633/101‑1 (MQ=38)
               ccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTgg                                                                                   >  2:3748391/1‑101 (MQ=255)
                 cgcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGc                                                                                  >  2:139197/5‑101 (MQ=38)
                 cgcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCTC.CTTTCACCTTGGc                                                                                  <  1:3822830/97‑1 (MQ=38)
                  gcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCt                                                                                 >  2:860631/4‑101 (MQ=255)
                  gcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCt                                                                                 <  1:2182846/98‑1 (MQ=38)
                   cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTc                                                                                <  1:2051419/99‑1 (MQ=39)
                    acccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCa                                                                               <  1:1962212/100‑1 (MQ=39)
                    acccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCa                                                                               <  2:2966091/100‑1 (MQ=39)
                     cccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCaa                                                                              <  1:1756297/101‑1 (MQ=255)
                     cccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCaa                                                                              <  2:774854/101‑1 (MQ=255)
                       cccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAat                                                                            <  1:139213/101‑1 (MQ=255)
                       cccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAat                                                                            <  2:1447549/101‑1 (MQ=255)
                        ccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAata                                                                           <  1:4899841/101‑1 (MQ=255)
                        ccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAata                                                                           <  2:3919128/101‑1 (MQ=255)
                            cGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATccc                                                                       <  2:1029025/101‑1 (MQ=255)
                             gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCg                                                                      <  1:2093180/101‑1 (MQ=255)
                             gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTAATGGCCCAAGAAACGCGGAA.CCGCCC.C.CTTTCACCTTGGCGCAACTCTCCCg                                                                      <  2:3341899/101‑1 (MQ=255)
                              tATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAACCCC.CC.C.CTTTCACCTTGGCTCAAATATCCCgg                                                                     >  1:2035846/1‑101 (MQ=255)
                              tATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCgg                                                                     >  2:2344725/1‑101 (MQ=255)
                               aTTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCggg                                                                    >  1:2288343/1‑101 (MQ=255)
                                tttATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCgggg                                                                   <  2:2164473/101‑1 (MQ=255)
                                tttATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCgggg                                                                   >  1:218562/1‑101 (MQ=255)
                                tttATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCgggg                                                                   >  1:4195748/1‑101 (MQ=255)
                                 ttATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCggggg                                                                  <  2:2190232/101‑1 (MQ=255)
                                 ttATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCggggg                                                                  >  1:2563483/1‑101 (MQ=255)
                                  tATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGt                                                                 >  2:4036391/1‑101 (MQ=255)
                                  tATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGt                                                                 >  1:1261360/1‑101 (MQ=255)
                                   aTCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTc                                                                <  2:1660950/101‑1 (MQ=255)
                                     cccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAAAGCCGAAGCGG.CC.C.CTTGCACCTTCGCTCAAATATCCCGGGGGTCCg                                                              >  2:849986/1‑101 (MQ=255)
                                      ccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCgg                                                             <  1:2447015/101‑1 (MQ=255)
                                         aGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCggggg                                                          >  1:254682/1‑101 (MQ=255)
                                          gAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGGTGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGc                                                         >  2:1589750/1‑101 (MQ=38)
                                          gAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGc                                                         <  2:2341415/101‑1 (MQ=255)
                                          gAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGc                                                         <  2:3005628/101‑1 (MQ=255)
                                           aaaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCg                                                        <  1:2129171/101‑1 (MQ=255)
                                            aaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCgg                                                       >  2:360495/1‑101 (MQ=255)
                                            aaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCgg                                                       <  2:90817/101‑1 (MQ=255)
                                            aaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAAAGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGGCCGGGGGCgg                                                       >  1:2290097/1‑101 (MQ=38)
                                             aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGa                                                      >  2:474331/1‑101 (MQ=255)
                                              aaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAg                                                     <  1:3063013/101‑1 (MQ=255)
                                               agcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..c                                                  >  2:1374580/1‑101 (MQ=255)
                                                gcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAggc                                                   >  2:3471801/1‑99 (MQ=38)
                                                gcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cc                                                 <  2:2288179/101‑1 (MQ=255)
                                                gcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cc                                                 <  1:3049709/101‑1 (MQ=255)
                                                gcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cc                                                 <  2:419578/101‑1 (MQ=255)
                                                gcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cc                                                 <  1:3748669/101‑1 (MQ=255)
                                                 cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..ccc                                                >  1:1976428/1‑101 (MQ=255)
                                                 cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CC.CGCTTTCACCTTGGCTCAAAGATCCCGGGGGTCCGGGGGCGGAG..ccc                                                >  1:4243165/1‑101 (MQ=37)
                                                  agcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccc                                               <  1:2196920/101‑2 (MQ=37)
                                                  agcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccc                                               <  1:860733/101‑2 (MQ=37)
                                                  agcGGCTTGCCCCGAAGGACTTG.TTTTAATGGCACAGGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGAGGAG..cccc                                               <  2:4150550/101‑2 (MQ=25)
                                                   gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCGAATATCCCGGGGGTCCGGGGGCGGAG..ccccc                                              <  1:4534055/101‑3 (MQ=37)
                                                   gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..ccccc                                              <  1:2344889/101‑3 (MQ=37)
                                                   gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..ccccc                                              <  2:2714406/101‑3 (MQ=37)
                                                   gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGGCCGGGGGCGGAG..ccccc                                              >  2:375640/1‑99 (MQ=37)
                                                        ttGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcgg                                         >  1:1118530/1‑94 (MQ=25)
                                                          gCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcgggg                                       >  2:2166333/1‑92 (MQ=21)
                                                          gCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcgggg                                       >  2:3480721/1‑92 (MQ=21)
                                                           ccccGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcgggcg                                      >  2:576145/1‑91 (MQ=21)
                                                                                           aCGCCGAC.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCCCGGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACAcccc       <  2:4712730/101‑1 (MQ=2)
                                                                                                gAA.CCGCCC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCCCGGCGGGGCGCCGCCGGG.G.AAGGCCGAGGCGGAAAAcccccccgc  >  1:638168/1‑96 (MQ=2)
                                                                                                       |                                                                                                  
GCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAA.CCG.CC.C.CTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..CCCGCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGC  >  NC_014034/2185830‑2186024

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