Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,467 C→T 18.3% T450I (ACC→ATC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,4670C→T18.3% 118.1 ‑1.8 ‑1.4234T450I (ACC→ATC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (22/21):  ref base (120/70):  total (142/92)
Fisher's exact test strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.42e-02

Only 100 of 247 total aligned reads displayed.

GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGAC  >  NC_014034/1397368‑1397566
                                                                                                   |                                                                                                      
gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAcc                                                                                                       >  1:3287716/1‑101 (MQ=255)
  cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.t                                                                                                    >  1:2972547/1‑101 (MQ=255)
  cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.t                                                                                                    >  1:229296/1‑101 (MQ=255)
   gCTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.tc                                                                                                   >  1:4492876/1‑101 (MQ=255)
    cTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.tca                                                                                                  <  1:3266670/101‑1 (MQ=255)
    cTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.tca                                                                                                  >  1:3046095/1‑101 (MQ=255)
     tCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.tcat                                                                                                 >  1:2841805/1‑101 (MQ=255)
      ccccccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCT.tcatc                                                                                                <  2:1788498/101‑1 (MQ=255)
      ccccccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.tcatc                                                                                                >  1:1785411/1‑101 (MQ=255)
         cccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATc.acc                                                                                            <  2:2516896/101‑1 (MQ=255)
          ccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATc.acca                                                                                           <  1:2440066/101‑1 (MQ=255)
          ccGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATc.acca                                                                                           <  2:4568029/101‑1 (MQ=255)
            gcgGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACg                                                                                         <  2:2040105/101‑1 (MQ=255)
              gggggCGTGGATCGAAACTCTTCAATGGATATGGCAAACCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCa                                                                                       <  2:1618078/101‑1 (MQ=255)
              gggggCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCa                                                                                       >  2:3758301/1‑101 (MQ=255)
               ggggCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCaa                                                                                      <  2:1002235/101‑1 (MQ=255)
                gtgCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.c                                                                                    >  2:3872733/3‑101 (MQ=255)
                 ggCGTCGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.cc                                                                                   <  2:1557134/101‑1 (MQ=255)
                   cGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGc                                                                                 <  1:1265888/101‑1 (MQ=255)
                    gTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGcc                                                                                >  1:1086981/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                                                                              >  2:2570419/1‑101 (MQ=255)
                      ggATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCt                                                                              <  2:3943861/101‑1 (MQ=255)
                      ggATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCt                                                                              >  1:2649625/1‑101 (MQ=255)
                      ggATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCt                                                                              <  1:1862253/101‑1 (MQ=255)
                      ggATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCt                                                                              <  1:2159478/101‑1 (MQ=255)
                       gATCGAACCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTg                                                                             <  2:4249523/101‑1 (MQ=255)
                       gATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTg                                                                             >  2:3684293/1‑101 (MQ=255)
                        aTCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGc                                                                            >  2:3501087/1‑101 (MQ=255)
                         tCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca                                                                           <  2:2989498/101‑1 (MQ=255)
                          cGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcac                                                                          >  2:329547/1‑101 (MQ=255)
                            aaaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacca                                                                        <  2:4093732/101‑1 (MQ=255)
                             aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaccac                                                                       <  1:468067/101‑1 (MQ=255)
                                  ttCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTc                                                                  >  1:1143913/1‑101 (MQ=255)
                                    cAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGc                                                                >  2:2481955/1‑101 (MQ=255)
                                      aTGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTc                                                              >  1:2532386/1‑101 (MQ=255)
                                      aTGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTc                                                              >  1:1696328/1‑101 (MQ=255)
                                        ggATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTccc                                                            <  2:2888277/101‑1 (MQ=255)
                                         gATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcccc                                                           >  1:1024407/1‑101 (MQ=255)
                                         gATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcccc                                                           <  2:4143128/101‑1 (MQ=255)
                                            atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                        <  2:2750708/101‑1 (MQ=255)
                                            atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                        <  1:822844/101‑1 (MQ=255)
                                              ggCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TC.TC.AACCCGCAA.CCGCCCCGCACCCCCGGGCGCCTCCCCCCCgg                                                     >  1:2762949/1‑101 (MQ=21)
                                                 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)
                                                       gTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTgg                                             >  2:2056976/1‑101 (MQ=255)
                                                        ccgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTg                                              <  2:1096182/91‑1 (MQ=0)
                                                         cgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTgg                                             >  2:2274767/10‑101 (MQ=0)
                                                         ccgcaccATGTCGCTCACCCCGCGGGGGCGTGGATCGGAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                           <  2:1900126/95‑1 (MQ=37)
                                                         ccGCGACATGTCGCTCTCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                           <  2:1528627/101‑1 (MQ=255)
                                                         ccGCGACATGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                           <  2:1589370/101‑1 (MQ=255)
                                                          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)
                                                            taccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                          >  1:567312/7‑101 (MQ=2)
                                                            cgACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                        >  1:1000454/1‑101 (MQ=255)
                                                              ccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                        >  1:2046598/5‑101 (MQ=2)
                                                              aCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      >  1:4672675/1‑101 (MQ=255)
                                                              aCATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGGGGATCGaaa                                      >  2:3638223/1‑101 (MQ=255)
                                                               cgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                       >  2:817919/4‑101 (MQ=2)
                                                               cATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                     >  2:2535641/1‑101 (MQ=255)
                                                                gcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      <  1:2342061/99‑1 (MQ=2)
                                                                gcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      >  1:3124721/3‑101 (MQ=2)
                                                                 tGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaacaa                                   >  1:4381244/1‑101 (MQ=255)
                                                                 tGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaacaa                                   >  2:1282712/1‑101 (MQ=255)
                                                                 cgtccgtCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                     <  2:79404/95‑1 (MQ=2)
                                                                 cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                     >  1:4422136/2‑101 (MQ=2)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggt                                  <  1:883732/101‑4 (MQ=2)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAccc                                >  1:1470355/1‑101 (MQ=255)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtg                                 >  1:2550345/1‑97 (MQ=2)
                                                                     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)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggggtc                               >  2:3624334/1‑95 (MQ=2)
                                                                         cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                             <  1:3502971/101‑9 (MQ=2)
                                                                         cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                             <  2:2626139/101‑9 (MQ=2)
                                                                         cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                             <  1:2724670/101‑9 (MQ=2)
                                                                         cccccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                           >  1:4256316/1‑101 (MQ=255)
                                                                         cccccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                           >  2:1119777/1‑101 (MQ=255)
                                                                          ccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccag                            <  1:720084/101‑10 (MQ=2)
                                                                          ccccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCAt                          >  2:726026/1‑101 (MQ=255)
                                                                           cccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                         <  1:2091823/101‑1 (MQ=255)
                                                                           cccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                         <  1:3247942/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                 >  1:1320637/1‑101 (MQ=255)
                                                                                    gCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCCCCGGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTc                <  1:3901285/101‑1 (MQ=255)
                                                                                    gCGTGGATCGAAACACCG.CCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTc                >  1:29844/1‑101 (MQ=255)
                                                                                      gTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAGACAACCCGATCATGGCCCGCCTCAg              >  2:3040432/1‑101 (MQ=255)
                                                                                      gTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAg              >  1:1718644/1‑101 (MQ=255)
                                                                                      gTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAg              >  1:3693404/1‑101 (MQ=255)
                                                                                        ggATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgcg            >  1:337130/1‑101 (MQ=255)
                                                                                               aaCACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCt     >  2:1734791/1‑101 (MQ=255)
                                                                                                  acCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCACCGCCGCGCTGAc  >  1:3786/1‑101 (MQ=255)
                                                                                                   |                                                                                                      
GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGAC  >  NC_014034/1397368‑1397566

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