Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,537 G→T 13.8% P473P (CCG→AAT)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,5370G→T13.8% 96.1 ‑1.7 ‑0.0297P473P (CCG→CCT) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (16/25):  ref base (130/126):  total (146/151)
Fisher's exact test strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.90e-01

Only 100 of 315 total aligned reads displayed.

GCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGC  >  NC_014034/1397437‑1397637
                                                                                                    |                                                                                                       
  tCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGAt                                                                                                       >  2:1136457/1‑101 (MQ=255)
   ccccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATc                                                                                                      <  2:3923732/101‑1 (MQ=255)
     ccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCAt                                                                                                    <  1:1994561/101‑1 (MQ=255)
      cccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                                                                   <  2:3734429/101‑1 (MQ=255)
       ccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATgg                                                                                                  <  1:3766829/101‑1 (MQ=255)
         gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c                                                                                               <  2:666211/101‑1 (MQ=255)
          cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c..c                                                                                            <  1:4287608/101‑1 (MQ=255)
           gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..Cg                                                                                           >  2:3353052/1‑101 (MQ=255)
           gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..Cg                                                                                           >  2:4089654/1‑101 (MQ=255)
             gggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGcc                                                                                         <  2:173657/101‑1 (MQ=255)
               gCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTc                                                                                       >  2:1900110/1‑101 (MQ=255)
                cGTGGATCGAAACACCGTCATCACCACGCACCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATAATGGC.C..CGCCTCa                                                                                      <  2:2894483/101‑1 (MQ=255)
                  tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                                    <  2:4141995/101‑1 (MQ=255)
                      tCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGcgccg                                                                                <  1:135578/101‑1 (MQ=255)
                      tCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGcgccg                                                                                >  1:4275591/1‑101 (MQ=255)
                       cGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGcgccgc                                                                               >  2:2087318/1‑101 (MQ=255)
                        gAAACACCGTCATCACCACGCAACCGTCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCcgcg                                                                              <  1:2693877/101‑1 (MQ=255)
                        gAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCcgcg                                                                              >  1:2937446/1‑101 (MQ=255)
                        gAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCcgcg                                                                              >  1:2818027/1‑101 (MQ=255)
                            cacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGa                                                                          <  1:3182106/101‑1 (MQ=255)
                                gTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTc                                                                      >  1:504215/1‑101 (MQ=255)
                                     caccacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcc                                                                 <  1:1001709/101‑1 (MQ=255)
                                      accacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccc                                                                >  1:692614/1‑101 (MQ=255)
                                       ccacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccc                                                               <  1:577133/101‑1 (MQ=255)
                                        cacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccccc                                                              <  1:1438461/101‑1 (MQ=255)
                                             aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgg                                                         <  2:1944497/101‑1 (MQ=255)
                                               ccGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCGTGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgggg                                                       <  1:3932774/101‑1 (MQ=255)
                                                   ccTGCCCCCCCGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGt                                                   <  1:67039/101‑1 (MQ=255)
                                                    cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg                                                  >  2:1644320/1‑101 (MQ=255)
                                                    cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg                                                  >  2:1306909/1‑101 (MQ=255)
                                                    cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg                                                  >  2:3158082/1‑101 (MQ=255)
                                                      tgccgctgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGa                                                <  1:220126/92‑1 (MQ=2)
                                                      gCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGa                                                <  1:1121442/101‑1 (MQ=255)
                                                      gCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGa                                                <  2:1592518/101‑1 (MQ=255)
                                                       caccacCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGAt                                               <  2:2238892/101‑1 (MQ=255)
                                                       caccacCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGAt                                               <  1:3543289/101‑1 (MQ=255)
                                                          aggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCCCGGGGGCGTGGAtgaca                                           >  2:3612589/6‑97 (MQ=0)
                                                            gatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGGGGGTCGGAAc                                         >  2:2431248/4‑101 (MQ=0)
                                                             cGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                         <  1:1272308/101‑1 (MQ=255)
                                                             cGGTCGCTCCACCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                         >  1:3849619/1‑101 (MQ=255)
                                                             cGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                         <  1:1144305/101‑1 (MQ=255)
                                                              ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                        <  1:647964/101‑1 (MQ=255)
                                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                       <  2:3456483/101‑1 (MQ=2)
                                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                        >  2:2944194/1‑101 (MQ=2)
                                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggt                                     <  2:2295930/101‑5 (MQ=2)
                                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                      >  2:4167502/1‑101 (MQ=255)
                                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                      >  2:1209741/1‑101 (MQ=255)
                                                                 ctcTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcc                                     <  2:1838111/99‑1 (MQ=255)
                                                                 cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                     <  2:2836700/101‑2 (MQ=2)
                                                                 cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                     >  2:4238696/1‑100 (MQ=2)
                                                                 cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                      >  1:4047448/1‑99 (MQ=2)
                                                                 cGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcc                                     <  1:1926276/101‑1 (MQ=255)
                                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaa                                     <  1:2569240/101‑4 (MQ=2)
                                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                    >  1:2597123/1‑101 (MQ=255)
                                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                    >  2:2977803/1‑101 (MQ=255)
                                                                  ccTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                   <  2:578656/100‑7 (MQ=2)
                                                                   ccccccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                    <  2:4411418/99‑5 (MQ=0)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGGGGGTCGaaaaggggcga                                  >  1:1694057/1‑93 (MQ=2)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGATCCCCCCGCGGGGGCGTGGATCGAAAcgggtcga                                  <  1:2133141/101‑8 (MQ=2)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgat                                   <  2:2152850/101‑4 (MQ=2)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                    >  2:3279913/1‑97 (MQ=2)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                    >  2:3063328/1‑97 (MQ=2)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaa                                 <  2:350176/101‑9 (MQ=2)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatg                                  >  1:4346928/1‑97 (MQ=2)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATcgatg                                  >  2:2410476/1‑97 (MQ=2)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaag                                   >  2:879322/1‑96 (MQ=0)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgcg                                  >  2:1980929/1‑101 (MQ=255)
                                                                    cccccccGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaa                                 <  1:2115922/100‑9 (MQ=2)
                                                                    cccccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaag                                   <  1:1654585/100‑6 (MQ=0)
                                                                     ccccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaagc                                  >  1:3327065/1‑95 (MQ=0)
                                                                     ccccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCgg                                 <  2:1109480/101‑1 (MQ=255)
                                                                      cccccGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaaga                               <  1:1524874/101‑11 (MQ=0)
                                                                      cccccGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaaga                               <  2:3531133/101‑11 (MQ=0)
                                                                        cccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatca                              >  1:2726014/1‑93 (MQ=0)
                                                                        cccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCATCCGCGGATc                              <  1:2422240/101‑1 (MQ=255)
                                                                        cccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATc                              <  1:3603453/101‑1 (MQ=255)
                                                                           gcgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGc                           >  1:2385005/1‑101 (MQ=255)
                                                                            cgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCg                          >  1:1776753/1‑101 (MQ=255)
                                                                             gggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGt                         >  1:1595403/1‑101 (MQ=255)
                                                                              ggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                        >  1:849278/1‑101 (MQ=255)
                                                                               gggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGa                       <  2:793693/101‑1 (MQ=255)
                                                                               gggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGa                       <  1:1847947/101‑1 (MQ=255)
                                                                               gggCGTGGATCGAAACAACCCGATCATGGC.C..CCCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGGGGATCGAAAAATCCGCGGATCGGCGTgc                       >  1:1087594/1‑100 (MQ=255)
                                                                                   gTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                   <  1:3729837/101‑1 (MQ=255)
                                                                                   gTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                   <  2:2700299/101‑1 (MQ=255)
                                                                                   gTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                   <  1:981317/101‑1 (MQ=255)
                                                                                    tGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgc                  >  1:3195153/1‑101 (MQ=255)
                                                                                       aTCGAAACACCCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc               <  1:2268745/101‑1 (MQ=255)
                                                                                       aTCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc               <  1:3510730/101‑1 (MQ=255)
                                                                                        tCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccg              >  2:4010154/1‑101 (MQ=255)
                                                                                        tCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccg              <  2:1416260/101‑1 (MQ=255)
                                                                                        tCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccg              <  1:742272/101‑1 (MQ=255)
                                                                                        tCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccc              >  2:4369634/1‑100 (MQ=255)
                                                                                         cGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccgc             <  2:547999/101‑1 (MQ=255)
                                                                                           aaaCAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc           <  2:2391187/101‑1 (MQ=255)
                                                                                              caacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc        >  1:3100135/1‑101 (MQ=255)
                                                                                               aacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCg       >  1:4178426/1‑101 (MQ=255)
                                                                                               aacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCCCCCATCg       >  1:3043434/1‑101 (MQ=255)
                                                                                               aacCCGATCATGCC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCg       <  2:694784/101‑1 (MQ=255)
                                                                                                 cccGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCCCCCATCGGt     >  1:1466391/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
GCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGC  >  NC_014034/1397437‑1397637

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