Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,469 +A 18.9% coding (1351/1578 nt) RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,4691.→A18.9% 184.6 ‑2.5 ‑0.5227coding (1351/1578 nt)RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (21/22):  ref base (117/67):  total (138/89)
Fisher's exact test strand distribution p-value = 8.43e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.11e-01

Only 100 of 250 total aligned reads displayed.

CGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTC  >  NC_014034/1397370‑1397569
                                                                                                    |                                                                                                      
cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.t                                                                                                       >  1:2972547/1‑101 (MQ=255)
  cTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.tca                                                                                                     >  1:3046095/1‑101 (MQ=255)
    ccccccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCT.tcatc                                                                                                   <  2:1788498/101‑1 (MQ=255)
    ccccccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.tcatc                                                                                                   <  1:853592/101‑1 (MQ=255)
     cccccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.tcatc.a                                                                                                 >  1:3751584/1‑101 (MQ=255)
       cccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATc.acc                                                                                               <  2:2516896/101‑1 (MQ=255)
        ccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATc.acca                                                                                              <  2:4568029/101‑1 (MQ=255)
        ccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATc.acca                                                                                              <  1:2440066/101‑1 (MQ=255)
          gcgGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACACCG.TCATC.ACCACg                                                                                            >  1:4093974/1‑101 (MQ=255)
          gcgGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACg                                                                                            <  2:2040105/101‑1 (MQ=255)
            gggggCGTGGATCGAAACTCTTCAATGGATATGGCAAACCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCa                                                                                          <  2:1618078/101‑1 (MQ=255)
            gggggCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCa                                                                                          >  2:3556512/1‑101 (MQ=255)
                 cGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGc                                                                                    >  2:1403032/1‑101 (MQ=255)
                 cGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGc                                                                                    <  1:3253397/101‑1 (MQ=255)
                 cGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGc                                                                                    <  1:1265888/101‑1 (MQ=255)
                    ggATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCt                                                                                 >  2:2570419/1‑101 (MQ=255)
                    ggATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCt                                                                                 >  1:3658974/1‑101 (MQ=255)
                    ggATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCt                                                                                 >  2:3104603/1‑101 (MQ=255)
                    ggATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCt                                                                                 >  2:468023/1‑101 (MQ=255)
                    ggATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCt                                                                                 <  1:1862253/101‑1 (MQ=255)
                     gATCGAACCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTg                                                                                <  2:4249523/101‑1 (MQ=255)
                     gATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTg                                                                                >  2:3684293/1‑101 (MQ=255)
                       tagAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca                                                                              <  1:4403631/99‑1 (MQ=255)
                       tCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca                                                                              >  1:599753/1‑101 (MQ=255)
                       tCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca                                                                              <  2:2989498/101‑1 (MQ=255)
                        cGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcac                                                                             >  1:1695934/1‑101 (MQ=255)
                          aaaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacca                                                                           >  1:3271817/1‑101 (MQ=255)
                          aaaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacca                                                                           <  2:4093732/101‑1 (MQ=255)
                           aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaccac                                                                          >  1:2127539/1‑101 (MQ=255)
                              tctTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCgg                                                                       >  1:2504695/1‑101 (MQ=255)
                                ttCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTc                                                                     >  1:1143913/1‑101 (MQ=255)
                                ttCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTc                                                                     >  1:2795081/1‑101 (MQ=255)
                                  cAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGc                                                                   >  2:2481955/1‑101 (MQ=255)
                                     tGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcc                                                                >  1:1292300/1‑101 (MQ=255)
                                       gATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcccc                                                              <  2:4143128/101‑1 (MQ=255)
                                       gATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcccc                                                              >  1:1024407/1‑101 (MQ=255)
                                        agatcgCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTccccc                                                             >  2:2919982/6‑101 (MQ=255)
                                         tatGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcccccc                                                            >  2:1563305/1‑101 (MQ=255)
                                          atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                           <  2:2750708/101‑1 (MQ=255)
                                               aaaaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgggg                                                      >  2:88896/1‑101 (MQ=255)
                                                 aaCCGTCCGCGACTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGc                                                    <  2:229271/101‑1 (MQ=255)
                                                 aaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGc                                                    <  2:2449850/101‑1 (MQ=255)
                                                    cGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTg                                                 >  1:1944511/1‑101 (MQ=255)
                                                     gTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTgg                                                >  2:2056976/1‑101 (MQ=255)
                                                      ccgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTg                                                 <  2:1096182/91‑1 (MQ=0)
                                                      ccgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTg                                                 <  2:155169/91‑1 (MQ=0)
                                                       cgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTgg                                                >  2:3767464/10‑101 (MQ=0)
                                                       cgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTgg                                                >  2:2274767/10‑101 (MQ=0)
                                                       cgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGGgg                                                >  1:2719598/10‑101 (MQ=0)
                                                       ccgcaccATGTCGCTCACCCCGCGGGGGCGTGGATCGGAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                              <  2:1900126/95‑1 (MQ=37)
                                                       ccGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGGGGAt                                              >  2:3663524/1‑101 (MQ=255)
                                                       ccGCGACATGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                              <  2:1589370/101‑1 (MQ=255)
                                                        gataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGa                                               >  1:2759379/9‑101 (MQ=0)
                                                        gataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGGGGa                                               >  1:479829/9‑101 (MQ=0)
                                                        cgcgACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                             >  2:960085/1‑101 (MQ=255)
                                                         gcgACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                            <  2:4492626/101‑1 (MQ=255)
                                                         ataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCTCCCGCGGGGGCGTGGAt                                              >  2:3783012/8‑101 (MQ=0)
                                                          ttccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                             <  2:4720600/95‑1 (MQ=2)
                                                          cgACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                           >  1:1000454/1‑101 (MQ=255)
                                                            aCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                         >  1:4672675/1‑101 (MQ=255)
                                                             cgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                          >  2:817919/4‑101 (MQ=2)
                                                              gcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                         <  1:2342061/99‑1 (MQ=2)
                                                               cttccgatctcccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaacaa                                      <  1:176345/91‑1 (MQ=255)
                                                               cgtccgtCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                        <  2:79404/95‑1 (MQ=2)
                                                               cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                        >  1:4422136/2‑101 (MQ=2)
                                                               cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                        <  2:4048450/100‑1 (MQ=2)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATCGG.CGCTCCCCCCGCGGGGGCGTGGATCGAaacgg                                      >  2:1251931/1‑99 (MQ=2)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggt                                     <  1:883732/101‑4 (MQ=2)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAccc                                   >  1:1470355/1‑101 (MQ=255)
                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtg                                    >  2:3405435/1‑97 (MQ=2)
                                                                    ccccccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgt                                   <  2:439335/99‑6 (MQ=2)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaaaggtgt                                   >  2:238838/1‑95 (MQ=2)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGa                                 <  2:599689/101‑1 (MQ=255)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGa                                 <  1:2535786/101‑1 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggggtc                                  >  2:3624334/1‑95 (MQ=2)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGAt                                >  1:4516962/1‑101 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGAt                                >  1:3246560/1‑101 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGAt                                >  1:10871/1‑101 (MQ=255)
                                                                       cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                                <  1:3502971/101‑9 (MQ=2)
                                                                       cccccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                              >  1:4256316/1‑101 (MQ=255)
                                                                        ccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccag                               <  1:720084/101‑10 (MQ=2)
                                                                        ccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccag                               <  2:2553251/101‑10 (MQ=2)
                                                                        ccccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCAt                             >  2:726026/1‑101 (MQ=255)
                                                                         cccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                            <  1:2091823/101‑1 (MQ=255)
                                                                              gggggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCg                       >  2:197033/1‑101 (MQ=255)
                                                                                 ggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGNNCGCCt                    >  1:2216197/1‑101 (MQ=255)
                                                                                 ggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCt                    <  2:2649471/101‑1 (MQ=255)
                                                                                  gCGTGGATCGAAACACCG.CCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTc                   >  1:29844/1‑101 (MQ=255)
                                                                                   cGTGGATCGAAACACCG.TCATC.CCCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCa                  <  1:2275203/101‑1 (MQ=255)
                                                                                    gTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAg                 >  1:1718644/1‑101 (MQ=255)
                                                                                    gTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAg                 >  2:3874729/1‑101 (MQ=255)
                                                                                     tGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgc                >  2:1004692/1‑101 (MQ=255)
                                                                                      ggATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgcg               >  1:337130/1‑101 (MQ=255)
                                                                                         tCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgccg            >  2:3086317/1‑101 (MQ=255)
                                                                                           gAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCcgcg          >  1:2627197/1‑101 (MQ=255)
                                                                                            aaaCACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCgcgc         >  2:3320444/1‑101 (MQ=255)
                                                                                             aaCACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCt        >  2:1734791/1‑101 (MQ=255)
                                                                                               cacCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGa      <  1:4712097/101‑1 (MQ=255)
                                                                                                acCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCACCGCCGCGCTGAc     >  1:3786/1‑101 (MQ=255)
                                                                                                  cG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGt   >  1:310501/1‑101 (MQ=255)
                                                                                                    |                                                                                                      
CGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTC  >  NC_014034/1397370‑1397569

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG