Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,538 A→G 14.7% M474V (ATC→GAC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,5380A→G14.7% 101.0 ‑0.5 ‑0.3278M474V (ATC→GTC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (21/20):  ref base (111/126):  total (132/146)
Fisher's exact test strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.92e-01

Only 100 of 299 total aligned reads displayed.

CTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCT  >  NC_014034/1397438‑1397638
                                                                                                    |                                                                                                         
cTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGa                                                                                                           <  1:1213951/101‑1 (MQ=255)
cTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGa                                                                                                           >  1:4109396/1‑101 (MQ=255)
  ccccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATc                                                                                                         >  2:2533829/1‑101 (MQ=255)
   cccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                                                                                                        <  2:3272670/101‑1 (MQ=255)
   cccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                                                                                                        <  1:1747084/101‑1 (MQ=255)
    ccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCAt                                                                                                       <  2:4182587/101‑1 (MQ=255)
     cccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCACGGGGGCGTGGATCGAAACAACCCGATCATg                                                                                                      >  1:3452440/1‑101 (MQ=255)
      ccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCAGgg                                                                                                     >  1:4556523/1‑101 (MQ=255)
        gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c                                                                                                  <  2:3093492/101‑1 (MQ=255)
        gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c                                                                                                  >  2:1962401/1‑101 (MQ=255)
         cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c..c                                                                                               >  1:1106602/1‑101 (MQ=255)
         cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c..c                                                                                               <  2:3237578/101‑1 (MQ=255)
         cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c..c                                                                                               <  1:620070/101‑1 (MQ=255)
            gggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGcc                                                                                            <  1:3106458/101‑1 (MQ=255)
              gCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTc                                                                                          <  2:2175556/101‑1 (MQ=255)
                gTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAg                                                                                        >  1:2976212/1‑101 (MQ=255)
                 tGGATCGAAACCCCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                                       <  1:3689834/101‑1 (MQ=255)
                 tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                                       >  1:3903885/1‑101 (MQ=255)
                 tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                                       <  1:139921/101‑1 (MQ=255)
                  ggATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgcg                                                                                      >  1:1866465/1‑101 (MQ=255)
                  ggATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgcg                                                                                      >  2:4555664/1‑101 (MQ=255)
                      cGAAAAACCGTCATCACCACGCAACCGCCCTGCACCACCGGCCGCTCCCCCCGCGGGGGCGGGGATCGAACCAACCCGACAATGCC.C..CCCCTAAGcgccgc                                                                                  >  2:2292805/1‑101 (MQ=37)
                        aaaCACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCgcgc                                                                                >  2:1902633/1‑101 (MQ=255)
                           cacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.Ga                                                                            >  2:2397761/1‑101 (MQ=255)
                              cGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGt                                                                        >  1:2712480/1‑101 (MQ=255)
                                     accacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTccc                                                                 <  1:3533878/101‑1 (MQ=255)
                                       cacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTccccc                                                               <  2:942914/101‑1 (MQ=255)
                                        ccGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTcccccc                                                              <  2:759481/100‑1 (MQ=255)
                                        acGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTcccccc                                                              <  2:2076277/101‑1 (MQ=255)
                                         cGCAACCGCCCTGCACCACCGGTCGCTCCCCCCATGGCGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCcg                                                             <  1:3986866/101‑1 (MQ=255)
                                          gCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCcgc                                                            >  2:1776701/1‑101 (MQ=255)
                                          gCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCcgc                                                            >  2:2873118/1‑101 (MQ=255)
                                            aaCCGCCCTGCGCCGCCGGTCCCACCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCgg                                                          <  2:234324/101‑1 (MQ=38)
                                            aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCgg                                                          <  1:2206861/101‑1 (MQ=255)
                                             aCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCggg                                                         <  1:1902768/101‑1 (MQ=255)
                                               cGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCggggg                                                       >  1:3172526/1‑101 (MQ=255)
                                                gCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGc                                                      >  1:3362416/1‑101 (MQ=255)
                                                   cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTg                                                   >  1:1805772/1‑101 (MQ=255)
                                                    tGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTTGCTCCCCCCGCGGGGGCGTgg                                                  >  2:425579/1‑101 (MQ=255)
                                                    tGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTgg                                                  >  1:2096380/1‑101 (MQ=255)
                                                    ggaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCT.AT.CGTCGCTCCCCCCGCGGGGGCGTGGa                                                 <  1:1662507/91‑1 (MQ=2)
                                                     gccccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGa                                                 <  1:3856181/98‑1 (MQ=255)
                                                     gCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGa                                                 <  2:3517195/101‑1 (MQ=255)
                                                      ccccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGAt                                                <  2:1592273/99‑1 (MQ=255)
                                                       atcgagcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGTCGCTCCCCCCGCGGGGGCGTGGAt                                                >  1:718833/8‑101 (MQ=0)
                                                       atcgagcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGTCGCTCCCCCCGCGGGGGCGTGGAt                                                >  1:4271121/8‑101 (MQ=0)
                                                       accaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATc                                               >  2:2986491/1‑101 (MQ=255)
                                                        gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGATGACCA.GAT.C..CGCAATCG.GCCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                              <  2:572672/96‑1 (MQ=0)
                                                        gcgctggTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTTGATGACCA.GAT.C..CGCAATCG.GCCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                              <  1:2837856/96‑1 (MQ=0)
                                                        cgctgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCg                                              <  1:1588754/95‑1 (MQ=2)
                                                          gagcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGTCGCTCCCCCCGCGGGGGCGGGGATCGa                                             >  1:786080/5‑101 (MQ=0)
                                                          ctgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            >  1:4180787/5‑101 (MQ=2)
                                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCT.AT.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggt                                      <  1:2166338/101‑5 (MQ=2)
                                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaa                                        >  2:3937079/1‑100 (MQ=2)
                                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                       >  1:436994/1‑101 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCT.AT.CGTCGCTCCCCCCCCGGGGGCGTGGGTCGAAAcggggc                                     >  1:3562191/1‑96 (MQ=2)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                      <  1:2927772/101‑2 (MQ=2)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                       >  1:2747654/1‑99 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCT.AT.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                    <  2:851291/101‑7 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcga                                     >  2:2015063/1‑99 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                      >  1:2784935/1‑98 (MQ=2)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgat                                    >  2:279151/1‑98 (MQ=2)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgat                                    >  1:202806/1‑98 (MQ=2)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                     <  1:3500734/101‑5 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtagatg                                   >  2:1847543/1‑96 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatg                                   >  2:1088950/1‑97 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaag                                    >  1:3407149/1‑96 (MQ=0)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgcg                                   >  2:262560/1‑101 (MQ=255)
                                                                   cccccccGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCT.AT.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaa                                  <  2:3617962/100‑9 (MQ=2)
                                                                   cccccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaag                                    <  1:2410790/100‑6 (MQ=0)
                                                                    ccccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatga                                  >  1:4088146/1‑96 (MQ=2)
                                                                    ccccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCgg                                  <  2:4109110/101‑1 (MQ=255)
                                                                     cccccGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCT.AT.CGTCGCTCCCCCCCCGGGGGCGGGGGTCGAAAcggggcgaaga                                >  1:1981652/1‑91 (MQ=0)
                                                                     cccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCT.GC.GGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaagcc                                  <  1:3729375/101‑8 (MQ=0)
                                                                      ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATCC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCCAAACATcgatgatc                                <  1:4617109/101‑8 (MQ=2)
                                                                      ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatc                                <  1:2553121/101‑8 (MQ=2)
                                                                      ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatc                                <  1:2628110/101‑8 (MQ=2)
                                                                      ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatc                                <  2:2188897/101‑8 (MQ=2)
                                                                      ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatc                                >  1:1722554/1‑94 (MQ=2)
                                                                      ccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAt                                <  2:2712284/101‑1 (MQ=255)
                                                                      ccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAt                                <  2:4043819/101‑1 (MQ=255)
                                                                       cccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATAC.CG.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatca                               <  2:4324849/101‑9 (MQ=0)
                                                                       cccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGACGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGagc                               >  2:4637007/1‑99 (MQ=255)
                                                                         cgcgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCgg                             <  1:2486039/101‑1 (MQ=255)
                                                                            gggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGt                          >  1:3048734/1‑101 (MQ=255)
                                                                              ccgatctgATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGa                        <  2:3381340/94‑1 (MQ=255)
                                                                               ggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGACGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAc                       >  2:714505/1‑101 (MQ=255)
                                                                                 cGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAccc                     <  2:3457423/101‑1 (MQ=255)
                                                                                  gtcgATCAAAACACCCCGATCCTGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                    <  2:4452179/98‑1 (MQ=38)
                                                                                   tGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGTGCCGAGCT.GA.CGTCGCTGCCCCCGCTGGGGCGTGGATCGAAACATCCGCGGAACGGCGTGAccccc                   >  2:3833732/1‑99 (MQ=38)
                                                                                   tGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgc                   >  1:655630/1‑101 (MQ=255)
                                                                                    ggATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcc                  >  1:3815906/1‑101 (MQ=255)
                                                                                     gATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGTCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccc                 <  1:4000622/101‑1 (MQ=255)
                                                                                      aTCGAACCAACCCTATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTTGATCGAAACATCCCCGGATCTGCGTGAcccgcccc                <  2:1224584/101‑1 (MQ=38)
                                                                                          aaaCAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc            >  2:1474120/1‑101 (MQ=255)
                                                                                             caacCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc         >  2:134150/1‑101 (MQ=255)
                                                                                             caacCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc         >  1:132750/1‑101 (MQ=255)
                                                                                                 ccGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTc     >  2:1537908/1‑101 (MQ=255)
                                                                                                 ccGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTc     <  2:102441/101‑1 (MQ=255)
                                                                                                    aTCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCt  <  2:1125639/101‑1 (MQ=255)
                                                                                                    |                                                                                                         
CTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCT.GA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCT  >  NC_014034/1397438‑1397638

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