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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,4930A→T20.5% 124.6 ‑1.8 ‑1.9210T459S (ACC→TCC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/23):  ref base (99/67):  total (119/91)
Fisher's exact test strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.34e-02

Only 100 of 245 total aligned reads displayed.

TCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1397393‑1397593
                                                                                                       |                                                                                                    
tCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca                                                                                                      >  1:599753/1‑101 (MQ=255)
tCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca                                                                                                      >  2:3674896/1‑101 (MQ=255)
tCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGca                                                                                                      >  1:181836/1‑101 (MQ=255)
 cGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcac                                                                                                     >  2:805669/1‑101 (MQ=255)
  gAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacc                                                                                                    <  1:4737428/101‑1 (MQ=255)
  gAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacc                                                                                                    >  2:3605571/1‑101 (MQ=255)
   aaaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacca                                                                                                   >  2:3007338/1‑101 (MQ=255)
   aaaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacca                                                                                                   <  2:4093732/101‑1 (MQ=255)
    aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaccac                                                                                                  <  1:468067/101‑1 (MQ=255)
    aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaccac                                                                                                  >  1:2127539/1‑101 (MQ=255)
          tCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCg                                                                                            >  2:4205751/1‑101 (MQ=255)
             aTGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTc                                                                                         >  1:2532386/1‑101 (MQ=255)
             aTGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTc                                                                                         >  1:1696328/1‑101 (MQ=255)
              tGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcc                                                                                        >  1:1292300/1‑101 (MQ=255)
               ggATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTccc                                                                                       <  2:866112/101‑1 (MQ=255)
                 agatcgCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTccccc                                                                                     >  2:2919982/6‑101 (MQ=255)
                  tatGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcccccc                                                                                    <  2:4373413/101‑1 (MQ=255)
                  tatGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcccccc                                                                                    >  2:1563305/1‑101 (MQ=255)
                   atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                                                   >  1:3268135/1‑101 (MQ=255)
                   atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                                                   <  1:822844/101‑1 (MQ=255)
                   atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                                                   <  2:2569508/101‑1 (MQ=255)
                   atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                                                   >  2:530240/1‑101 (MQ=255)
                     ggCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TC.TC.AACCCGCAA.CCGCCCCGCACCCCCGGGCGCCTCCCCCCCgg                                                                                >  1:2762949/1‑101 (MQ=21)
                      gCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgg                                                                                <  1:2338710/101‑1 (MQ=255)
                        aaaaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgggg                                                                              >  2:4509547/1‑101 (MQ=255)
                          aaCCGTCCGCGACTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGc                                                                            <  2:229271/101‑1 (MQ=255)
                           aCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCg                                                                           <  2:2504551/101‑1 (MQ=255)
                             cGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTg                                                                         <  1:3556731/101‑1 (MQ=255)
                             cGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTg                                                                         >  2:2003457/1‑101 (MQ=255)
                              gTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTgg                                                                        >  2:2056976/1‑101 (MQ=255)
                               ccgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTg                                                                         <  2:1660860/91‑1 (MQ=0)
                                cgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                        >  2:2274767/10‑101 (MQ=0)
                                ccGCGACATGTCGCTCTCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      <  2:1528627/101‑1 (MQ=255)
                                ccGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      >  1:1539629/1‑101 (MQ=255)
                                ccGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      <  1:1020711/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)
                                   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)
                                      cATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                >  2:6882/1‑101 (MQ=255)
                                      cATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                >  2:2535641/1‑101 (MQ=255)
                                       gcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                 >  1:3124721/3‑101 (MQ=2)
                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaacaac                                                             >  2:4686575/1‑101 (MQ=255)
                                          tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATCGG.CGCTCCCCCCGCGGGGGCGTGGATCGAaacgg                                                              >  2:1251931/1‑99 (MQ=2)
                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtg                                                            >  2:782510/1‑97 (MQ=2)
                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtg                                                            >  2:3405435/1‑97 (MQ=2)
                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCg                                                          <  1:3104803/101‑1 (MQ=255)
                                             cTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaaaggtgt                                                           >  2:238838/1‑95 (MQ=2)
                                             cTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgt                                                           >  2:1868245/1‑96 (MQ=2)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGAt                                                        >  1:4516962/1‑101 (MQ=255)
                                               ccccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcc                                                         >  1:625758/1‑94 (MQ=2)
                                                cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                                                        <  1:2724670/101‑9 (MQ=2)
                                                cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                                                        <  1:3502971/101‑9 (MQ=2)
                                                cccccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                                                      >  2:1119777/1‑101 (MQ=255)
                                                cccccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                                                      >  1:4256316/1‑101 (MQ=255)
                                                  cccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaacaaagatcggaag                                                    >  2:176320/1‑91 (MQ=255)
                                                  cccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                    >  1:4411862/1‑101 (MQ=255)
                                                     gcgGGGGCGTGGATCGAAACACCG.CCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGcc                                                 <  2:2794916/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:1066458/1‑101 (MQ=255)
                                                          ggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCt                                            <  2:2649471/101‑1 (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.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAg                                         >  1:3693404/1‑101 (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)
                                                                 aTCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgcc                                     <  1:1106155/101‑1 (MQ=255)
                                                                  tCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgccg                                    >  2:3086317/1‑101 (MQ=255)
                                                                    gAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCcgcg                                  >  1:1068291/1‑101 (MQ=255)
                                                                    gAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCcgcg                                  >  1:2627197/1‑101 (MQ=255)
                                                                     aaaCACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCgcgc                                 >  2:3320444/1‑101 (MQ=255)
                                                                       acacCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTg                               <  1:2075332/101‑1 (MQ=255)
                                                                         acCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCACCGCCGCGCTGAc                             >  1:3786/1‑101 (MQ=255)
                                                                            g.TCCTC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTc                          <  1:2453368/101‑1 (MQ=255)
                                                                            g.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTc                          >  2:4057934/1‑101 (MQ=255)
                                                                               catc.aCCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGc                        <  1:115198/101‑1 (MQ=255)
                                                                                 tc.aCCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTc                      <  2:4479429/101‑1 (MQ=255)
                                                                                  c.accacGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcc                     >  1:4657170/1‑101 (MQ=255)
                                                                                    accacGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTccc                    <  1:329581/101‑1 (MQ=255)
                                                                                     ccacGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGATGTCGCTcccc                   <  2:4672419/101‑1 (MQ=255)
                                                                                       acGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCGGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcccccc                 <  1:945139/101‑1 (MQ=255)
                                                                                       acGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcccccc                 >  1:2841800/1‑101 (MQ=255)
                                                                                       acGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcccccc                 <  1:1403126/101‑1 (MQ=255)
                                                                                       acGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAAACCCGATCATCGACCGCCTCAGCGGCGCGCTGATGTCGCTtccccc                 >  2:2453389/1‑101 (MQ=37)
                                                                                          cAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgcg              <  2:3751349/101‑1 (MQ=255)
                                                                                          cAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgcg              >  2:3845235/1‑101 (MQ=255)
                                                                                           aa.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgg             >  2:3283618/1‑101 (MQ=255)
                                                                                                gCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGc         <  1:2570571/101‑1 (MQ=255)
                                                                                                  ccTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGt       <  1:2920158/101‑1 (MQ=255)
                                                                                                    tGCCCCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTgg     <  2:2404256/101‑1 (MQ=255)
                                                                                                    tGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTgg     >  1:1253985/1‑101 (MQ=255)
                                                                                                    tGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTgg     <  1:240174/101‑1 (MQ=255)
                                                                                                     gCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGa    >  2:4407981/1‑101 (MQ=255)
                                                                                                       accaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATc  <  1:774593/101‑1 (MQ=255)
                                                                                                       |                                                                                                    
TCGAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1397393‑1397593

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