Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,536 C→A 9.3% P473Q (CCG→AAT)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,5360C→A 9.3% 74.1 ‑1.4 ‑0.3335P473Q (CCG→CAG) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/20):  ref base (141/155):  total (155/180)
Fisher's exact test strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.62e-01

Only 100 of 361 total aligned reads displayed.

CGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCG  >  NC_014034/1397436‑1397636
                                                                                                     |                                                                                                       
cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAc.cc                                                                                                         <  1:721314/101‑1 (MQ=255)
 cccccccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCg                                                                                                        <  1:1199049/98‑1 (MQ=255)
    ccccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATc                                                                                                     >  1:3847344/1‑101 (MQ=255)
    ccccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATc                                                                                                     >  1:3712912/1‑101 (MQ=255)
     cccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCa                                                                                                    >  1:3822330/1‑101 (MQ=255)
     cccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCa                                                                                                    <  1:5096106/101‑1 (MQ=255)
     cccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCa                                                                                                    >  1:5469949/1‑101 (MQ=255)
       cccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATg                                                                                                  >  2:302400/1‑101 (MQ=255)
       cccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGACCATg                                                                                                  >  1:4420752/1‑101 (MQ=255)
          gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c                                                                                              >  1:3202834/1‑101 (MQ=255)
          gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c                                                                                              <  1:1256653/101‑1 (MQ=255)
          gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c                                                                                              <  1:776627/101‑1 (MQ=255)
           cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c..c                                                                                           <  1:5223908/101‑1 (MQ=255)
           cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c..c                                                                                           >  2:5324094/1‑101 (MQ=255)
           cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGc.c..c                                                                                           <  2:3362940/101‑1 (MQ=255)
             ggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGc                                                                                         <  2:3999705/101‑1 (MQ=255)
              gggCGTGGATGGAAACACCGTCATCACCACGCAACCGCCCTGCACCGGCGGTAGGTCACCCCGCGGGGGCGTGGATCGAAACAAA.CCGATCATGGC.C..CGcc                                                                                        >  2:4809728/1‑101 (MQ=37)
              gggCGTGGATCGAAACCCCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGcc                                                                                        <  2:5794635/101‑1 (MQ=255)
               ggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCt                                                                                       <  2:326086/101‑1 (MQ=255)
                 cGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCa                                                                                     <  2:1292565/101‑1 (MQ=255)
                    ggATCGAAACACCGTCATCACTACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAgcg                                                                                  >  1:3292788/1‑101 (MQ=255)
                    ggATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAgcg                                                                                  <  1:1385918/101‑1 (MQ=255)
                     gATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAgcgc                                                                                 >  2:25384/1‑101 (MQ=255)
                       tCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGcgccg                                                                               >  2:1021237/1‑101 (MQ=255)
                           gctcttccgatctCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCt                                                                           <  2:410508/92‑1 (MQ=38)
                           aaCACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCt                                                                           <  1:236997/101‑1 (MQ=255)
                            acacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTg                                                                          <  1:1104576/101‑1 (MQ=255)
                               ccGTCATCACCACGCACCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACg                                                                       <  2:3767670/101‑1 (MQ=255)
                                 gTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACgcc                                                                     >  2:2835101/1‑99 (MQ=255)
                                    atcaCCACGCAACCGCCCTGCACCCCCGGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCt                                                                  <  1:1481508/101‑1 (MQ=255)
                                    atcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCt                                                                  >  1:2607650/1‑101 (MQ=255)
                                    atcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCt                                                                  <  1:2270630/101‑1 (MQ=255)
                                     tcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTc                                                                 <  2:2056374/101‑1 (MQ=255)
                                      cactacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcc                                                                >  2:4322439/1‑101 (MQ=255)
                                          acGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                            <  1:5364708/101‑1 (MQ=255)
                                            gCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgc                                                          <  2:2355602/101‑1 (MQ=255)
                                                    ccTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGt                                                  >  2:3666177/1‑101 (MQ=255)
                                                      ggaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGa                                               <  1:1010202/91‑1 (MQ=2)
                                                       gCACCACCGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGa                                               <  2:5194685/101‑1 (MQ=255)
                                                        caccacCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGAt                                              <  1:1679812/101‑1 (MQ=255)
                                                          tcgagcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                             <  2:5621194/95‑1 (MQ=0)
                                                          ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                            <  1:2222844/101‑1 (MQ=255)
                                                           cgagcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                            <  1:3618269/96‑1 (MQ=0)
                                                           ccccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                           <  1:596845/99‑1 (MQ=255)
                                                           accgcgTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCGATGATCAAATC.C..GCAATCGGCCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                            <  1:1553335/96‑1 (MQ=0)
                                                             tgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                         <  2:3320313/98‑1 (MQ=2)
                                                              cGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                        >  1:5588896/1‑101 (MQ=255)
                                                               ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCGCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                       <  1:929925/101‑1 (MQ=255)
                                                               cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                        <  1:82644/100‑1 (MQ=2)
                                                               cgTCGCCCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                        <  2:4198967/100‑1 (MQ=2)
                                                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                      <  1:4666627/101‑1 (MQ=2)
                                                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                      >  2:1663266/1‑101 (MQ=2)
                                                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                      <  2:4417526/101‑1 (MQ=255)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                     <  1:338403/101‑1 (MQ=2)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                     >  1:1904913/1‑101 (MQ=255)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                    <  2:137027/101‑2 (MQ=2)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                     >  2:4682133/1‑99 (MQ=2)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                     <  1:2082591/101‑3 (MQ=2)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                     <  1:3035318/101‑3 (MQ=2)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                     >  2:1224800/1‑99 (MQ=2)
                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                  <  1:4781901/101‑7 (MQ=2)
                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                    <  2:3507299/101‑4 (MQ=2)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATcgat                                  >  2:3565456/1‑98 (MQ=2)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACACC.AATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                   <  2:5069185/101‑5 (MQ=2)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgc                                  <  1:4905232/101‑1 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACTAT.GCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaa                                <  2:1792929/101‑9 (MQ=2)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatg                                 >  1:2855559/1‑97 (MQ=2)
                                                                     cccccccGCGGGGGCGTGGATCGAATCAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgcg                                 <  1:4164889/100‑1 (MQ=255)
                                                                     cccccccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatg                                 <  2:5013856/100‑5 (MQ=2)
                                                                      ccccccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCgg                                >  1:3733821/1‑101 (MQ=255)
                                                                       cccccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgat                               <  1:3565673/101‑7 (MQ=2)
                                                                       cccccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgat                               >  1:3913391/1‑95 (MQ=2)
                                                                        ccccGCGGGGGCGTTGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAt                              <  2:3852165/101‑1 (MQ=255)
                                                                        ccccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAt                              <  1:2780411/101‑1 (MQ=255)
                                                                         cccGCGGGGGCGTGGATCGAAACACC.CCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatca                             <  1:5667739/101‑9 (MQ=0)
                                                                         cccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATc                             <  2:550667/101‑1 (MQ=255)
                                                                          ccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCg                            <  2:2187822/101‑1 (MQ=255)
                                                                          ccGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCg                            >  2:4692683/1‑101 (MQ=255)
                                                                            gcgGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGc                          <  1:1745408/101‑1 (MQ=255)
                                                                             cgGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCg                         >  1:4616252/1‑101 (MQ=255)
                                                                              gggggCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGt                        >  2:5447002/1‑101 (MQ=255)
                                                                               ggggCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                       <  2:1941575/101‑1 (MQ=255)
                                                                                    ggggATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                  <  2:1492902/99‑1 (MQ=255)
                                                                                    gTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                  <  2:4953732/101‑1 (MQ=255)
                                                                                    gTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                  <  2:1714259/101‑1 (MQ=255)
                                                                                     tGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgc                 >  1:4111341/1‑101 (MQ=255)
                                                                                      ggATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccccc                >  2:723943/1‑101 (MQ=255)
                                                                                        ctCGAAACAAC.CCGGTCATGGC.C..CGCCTCAGCGCCGCGCTGACGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc              <  2:2874002/100‑1 (MQ=255)
                                                                                        aTCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc              <  1:2158737/101‑1 (MQ=255)
                                                                                            aaaCAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCGCCGACCGGCGTGACCcgccccgccc          >  2:3159004/1‑101 (MQ=255)
                                                                                              acacc.CCGATCAGGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCAt        <  1:2523282/101‑1 (MQ=255)
                                                                                              acaac.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgcccccccccc        >  1:4264022/1‑99 (MQ=255)
                                                                                              acaac.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCAt        <  2:3822088/101‑1 (MQ=255)
                                                                                              acaac.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCAt        >  2:1541827/1‑101 (MQ=255)
                                                                                                aac.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCg      <  2:3129213/101‑1 (MQ=255)
                                                                                                aac.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCg      >  2:4622675/1‑101 (MQ=255)
                                                                                                 ac.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCgg     >  2:2436555/1‑101 (MQ=255)
                                                                                                 ac.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCgg     <  2:734497/101‑1 (MQ=255)
                                                                                                     cGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCg  >  2:2663564/1‑101 (MQ=255)
                                                                                                     cGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAGCGTCGTGGCACGCCACGCCCCTCGGtag  >  2:888457/1‑99 (MQ=25)
                                                                                                     |                                                                                                       
CGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAC.CCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCG  >  NC_014034/1397436‑1397636

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG