Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,493 A→T 15.4% T459S (ACC→TCG)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,4930A→T15.4% 103.2 ‑1.9 ‑1.8233T459S (ACC→TCC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (16/20):  ref base (115/82):  total (131/102)
Fisher's exact test strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.54e-02

Only 100 of 267 total aligned reads displayed.

TCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1397393‑1397593
                                                                                                       |                                                                                                    
 cGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcac                                                                                                     >  2:577065/1‑101 (MQ=255)
  gAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaca                                                                                                    >  1:4227759/1‑100 (MQ=255)
   aaaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacca                                                                                                   >  1:4324662/1‑101 (MQ=255)
   aaaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacca                                                                                                   <  2:3224688/101‑1 (MQ=255)
   aaaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacca                                                                                                   >  1:1944626/1‑101 (MQ=255)
    aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaccac                                                                                                  <  2:2409425/101‑1 (MQ=255)
    aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaccac                                                                                                  <  2:2579823/101‑1 (MQ=255)
    aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaccac                                                                                                  >  2:2550605/1‑101 (MQ=255)
     aCTNTNCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCaccacc                                                                                                 <  1:1489565/101‑1 (MQ=255)
      ctctTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCg                                                                                                >  2:3684923/1‑101 (MQ=255)
        ctTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGt                                                                                              >  1:793775/1‑101 (MQ=255)
        ctTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGt                                                                                              <  2:844419/101‑1 (MQ=255)
        ctTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCACCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGt                                                                                              >  2:2279656/1‑101 (MQ=255)
        ctTCAATGGATAAGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGt                                                                                              >  1:2419396/1‑101 (MQ=255)
           cAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGc                                                                                           >  1:3475342/1‑101 (MQ=255)
             aTGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTc                                                                                         >  1:3891165/1‑101 (MQ=255)
             aTGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTc                                                                                         <  2:1412342/101‑1 (MQ=255)
              tGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcc                                                                                        >  1:1761681/1‑101 (MQ=255)
               ggATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTccc                                                                                       >  2:1822211/1‑101 (MQ=255)
                   atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                                                   <  2:951813/101‑1 (MQ=255)
                     ggCACAACCGTCCGCGACATGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcgcg                                                                                 <  2:2001456/101‑1 (MQ=255)
                      gCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgg                                                                                >  1:1988468/1‑101 (MQ=255)
                       cAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCggg                                                                               >  2:3922258/1‑101 (MQ=255)
                       cAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCggg                                                                               <  2:1111086/101‑1 (MQ=255)
                        aaaaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgggg                                                                              >  2:3112576/1‑101 (MQ=255)
                          aaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGc                                                                            >  2:1665025/1‑101 (MQ=255)
                          aaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGc                                                                            >  1:2036230/1‑101 (MQ=255)
                           aCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCg                                                                           >  1:3299984/1‑101 (MQ=255)
                              gTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTgg                                                                        >  1:4075263/1‑101 (MQ=255)
                                 gataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGGGGa                                                                       >  1:3149864/9‑101 (MQ=0)
                                 cgcgACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                     <  1:4258147/101‑1 (MQ=255)
                                 cgcgACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                     <  1:599647/101‑1 (MQ=255)
                                   taccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                     <  2:3740394/95‑1 (MQ=2)
                                   taccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGGGGATc                                                                     >  2:1365416/7‑101 (MQ=0)
                                    accgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGGGGATCg                                                                    >  1:2096492/6‑101 (MQ=0)
                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacg                                                               >  1:2208601/1‑100 (MQ=2)
                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaacaac                                                             <  1:2550768/101‑1 (MQ=255)
                                          tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacgg                                                              <  2:2725824/101‑3 (MQ=2)
                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggg                                                             >  2:3565193/1‑98 (MQ=2)
                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAccc                                                           >  2:2884051/1‑101 (MQ=255)
                                             cTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggggt                                                           >  1:1884655/1‑96 (MQ=2)
                                             cTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGa                                                         >  1:548051/1‑101 (MQ=255)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGAt                                                        >  2:1136457/1‑101 (MQ=255)
                                                cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                                                        <  1:3026767/101‑9 (MQ=2)
                                                 ccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccag                                                       <  2:2083408/101‑10 (MQ=2)
                                                 ccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccag                                                       <  2:4188379/101‑10 (MQ=2)
                                                 ccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccag                                                       <  1:3681798/101‑10 (MQ=2)
                                                 ccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccag                                                       <  1:3561874/101‑10 (MQ=2)
                                                  cccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccagc                                                      <  1:2759917/101‑11 (MQ=2)
                                                  cccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                    >  1:4280316/1‑101 (MQ=255)
                                                   ccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATgg                                                   <  1:3766829/101‑1 (MQ=255)
                                                   ccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATgg                                                   >  1:2391337/1‑101 (MQ=255)
                                                    cgcgGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc                                                  <  1:2809044/101‑1 (MQ=255)
                                                     gcgGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGcc                                                 <  2:4051855/101‑1 (MQ=255)
                                                      cgGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGccc                                                <  2:2895306/101‑1 (MQ=255)
                                                       gggggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGcccc                                               >  2:2181596/1‑100 (MQ=255)
                                                       gggggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCg                                               >  2:1272215/1‑101 (MQ=255)
                                                       gggggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCg                                               >  2:3353052/1‑101 (MQ=255)
                                                         gggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGcc                                             <  2:173657/101‑1 (MQ=255)
                                                         gggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGcc                                             <  2:2426866/101‑1 (MQ=255)
                                                          ggCGTGGATCGAAACACCG.TCATC.CCCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCt                                            <  2:1960397/101‑1 (MQ=255)
                                                          ggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCt                                            <  2:1602395/101‑1 (MQ=255)
                                                           gCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTc                                           >  2:1900110/1‑101 (MQ=255)
                                                            cGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCa                                          >  2:484441/1‑101 (MQ=255)
                                                              tGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgc                                        >  1:1465181/1‑101 (MQ=255)
                                                              tGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgc                                        >  1:762079/1‑101 (MQ=255)
                                                               ggATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgcg                                       >  1:579157/1‑101 (MQ=255)
                                                                gATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgcgc                                      >  1:1908538/1‑101 (MQ=255)
                                                                 aTCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTNNNngnc                                     >  1:740612/1‑101 (MQ=255)
                                                                  tCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgccg                                    <  1:135578/101‑1 (MQ=255)
                                                                   cGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgccgc                                   <  1:1099134/101‑1 (MQ=255)
                                                                   cGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgccgc                                   >  1:1676487/1‑101 (MQ=255)
                                                                   cGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgccgc                                   >  1:2829438/1‑101 (MQ=255)
                                                                    gAAACACCG.TCATC.ACCACGCAA.CCGTCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCcgcg                                  <  1:2693877/101‑1 (MQ=255)
                                                                    gAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCcgcg                                  >  1:2937446/1‑101 (MQ=255)
                                                                      aaCACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCt                                >  2:666550/1‑101 (MQ=255)
                                                                      aaCACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCt                                >  1:2161093/1‑101 (MQ=255)
                                                                         acCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGAc                             <  1:889485/101‑1 (MQ=255)
                                                                           cG.TCTTC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCTGGGGCGTGGATCGAAACAACCCGACCATGGCCCGCCTCAGGGCCGCGCTGACga                           >  2:712878/1‑100 (MQ=38)
                                                                           cG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGt                           >  1:165736/1‑101 (MQ=255)
                                                                                atc.aCCACGCAA.CCGCCCTTCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCt                       <  1:4428433/101‑1 (MQ=255)
                                                                                 tc.aCCACGCAA.CCGCCGTGCCCCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTc                      <  2:3257767/101‑1 (MQ=255)
                                                                                  c.accacGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcc                     <  1:1001709/101‑1 (MQ=255)
                                                                                      cacGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTccccc                  <  1:1438461/101‑1 (MQ=255)
                                                                                      cacGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTccccc                  <  1:131714/101‑1 (MQ=255)
                                                                                       acGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcccccc                 >  2:4002436/1‑101 (MQ=255)
                                                                                        cGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcg                <  2:3296020/101‑1 (MQ=255)
                                                                                              ccGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCGTGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgggg           <  1:3932774/101‑1 (MQ=255)
                                                                                              ccGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgggg           >  2:3024173/1‑101 (MQ=255)
                                                                                              ccGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgggg           <  1:1563843/101‑1 (MQ=255)
                                                                                               cGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGTGCCTCGCTGACGTCGCTCCCCCCGCggggg          >  2:1468803/1‑101 (MQ=255)
                                                                                                gCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGc         <  2:4260715/101‑1 (MQ=255)
                                                                                                  cccgCACCACCGGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGt       <  1:484493/98‑1 (MQ=255)
                                                                                                  ccTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGt       >  2:3495317/1‑101 (MQ=255)
                                                                                                   cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg      >  2:3158082/1‑101 (MQ=255)
                                                                                                   cTGCACCACCGGTCGCTCCCCCCGCGAGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg      <  1:3733532/101‑1 (MQ=255)
                                                                                                    tGCCCCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTgg     <  1:3099114/101‑1 (MQ=255)
                                                                                                    tGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTgg     >  1:318928/1‑101 (MQ=255)
                                                                                                    tGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCAGATCATGGCCCGCCTCAGCGCCGCGCTTACGTCGCTCCCCCCGCGCGGGCGTgg     >  2:4115048/1‑101 (MQ=255)
                                                                                                     gCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGa    <  1:1121442/101‑1 (MQ=255)
                                                                                                       |                                                                                                    
TCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1397393‑1397593

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG