Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,532 A→T 7.6% N472Y (AAC→TCT)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,5320A→T 7.6% 51.6 ‑0.7 ‑1.8275N472Y (AAC→TAC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/13):  ref base (121/131):  total (129/146)
Fisher's exact test strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.49e-02

Only 100 of 295 total aligned reads displayed.

ATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCG  >  NC_014034/1397432‑1397632
                                                                                                    |                                                                                                       
atgccccTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaaca                                                                                                         <  1:157763/95‑1 (MQ=255)
aTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaaca                                                                                                         >  2:1682912/1‑101 (MQ=255)
aTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaaca                                                                                                         <  2:3738419/101‑1 (MQ=255)
aTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaaca                                                                                                         >  2:3924702/1‑101 (MQ=255)
     gCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCg                                                                                                    >  1:168306/1‑101 (MQ=255)
      cTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGa                                                                                                   <  1:1213951/101‑1 (MQ=255)
      cTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGa                                                                                                   >  1:4109396/1‑101 (MQ=255)
      cTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGa                                                                                                   >  2:3986581/1‑101 (MQ=255)
        ccccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATc                                                                                                 <  2:1579259/101‑1 (MQ=255)
        ccccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATc                                                                                                 >  2:2533829/1‑101 (MQ=255)
          ccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCAt                                                                                               <  2:3118327/101‑1 (MQ=255)
            ccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATgg                                                                                             >  1:2076417/1‑101 (MQ=255)
              gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c                                                                                          >  1:1641126/1‑101 (MQ=255)
               cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c..c                                                                                       >  1:1106602/1‑101 (MQ=255)
                gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..Cg                                                                                      <  2:3382273/101‑1 (MQ=255)
                  gggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGcc                                                                                    <  1:3106458/101‑1 (MQ=255)
                    gCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTc                                                                                  <  2:1917865/101‑1 (MQ=255)
                     cGTGGATCGAAACACCGTCATCCCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCa                                                                                 <  2:4072528/101‑1 (MQ=255)
                      gTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAg                                                                                >  1:2949694/1‑101 (MQ=255)
                      gTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAg                                                                                >  1:2745492/1‑101 (MQ=255)
                       tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                               >  1:3903885/1‑101 (MQ=255)
                       tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                               <  1:139921/101‑1 (MQ=255)
                       cgatctgAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                               <  1:1847999/95‑1 (MQ=255)
                               aaCCCCGTCATCACCACGCACCCGCCCTTCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCt                                                                       <  1:1897870/101‑1 (MQ=255)
                                acacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTg                                                                      <  2:3291178/101‑1 (MQ=255)
                                 cacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGa                                                                     >  2:2397761/1‑101 (MQ=255)
                                    cGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGt                                                                  >  1:2712480/1‑101 (MQ=255)
                                    cGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGt                                                                  >  1:1728192/1‑101 (MQ=255)
                                          caccgcGCACCCGTCCTCCCCCTCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcc                                                            <  2:2744553/101‑1 (MQ=37)
                                          caccacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcc                                                            >  1:1758253/1‑101 (MQ=255)
                                           accccGCACCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccc                                                           <  2:601023/101‑1 (MQ=255)
                                             cacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccccc                                                         <  2:3066163/101‑1 (MQ=255)
                                             cacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccccc                                                         <  2:201182/101‑1 (MQ=255)
                                             cacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccccc                                                         <  2:942914/101‑1 (MQ=255)
                                              ccGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                        <  2:759481/100‑1 (MQ=255)
                                               cGCAACCGCCCTGCACCACCGGTCGCTCCCCCCATGGCGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcg                                                       <  1:3986866/101‑1 (MQ=255)
                                                gCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgc                                                      >  2:1776701/1‑101 (MQ=255)
                                                 cAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgcg                                                     >  2:4573429/1‑101 (MQ=255)
                                                  aaCCGCCCTGCGCCGCCGGTCCCACCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgg                                                    <  2:234324/101‑1 (MQ=38)
                                                  aaCCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgg                                                    <  1:2206861/101‑1 (MQ=255)
                                                     cGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCggggg                                                 >  1:3172526/1‑101 (MQ=255)
                                                      cccgtGCGCCACCGGTCGCCCCCCCGTCGGGGGCGTGGATCGAAACAGCCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCGCCCCCCGCGGGGGc                                                <  1:3548910/97‑1 (MQ=37)
                                                         cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg                                             <  1:2022339/101‑1 (MQ=255)
                                                           gccccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGa                                           <  1:3856181/98‑1 (MQ=255)
                                                           gCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGa                                           <  2:3517195/101‑1 (MQ=255)
                                                              ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                        <  2:4293502/101‑1 (MQ=255)
                                                               caccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGGGTGGATCGa                                       >  2:2897363/1‑101 (MQ=255)
                                                               caccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                       >  2:1060705/1‑101 (MQ=255)
                                                                ctgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                      >  1:4180787/5‑101 (MQ=2)
                                                                 tgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                     <  2:2065067/98‑1 (MQ=2)
                                                                   cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                    <  1:3808232/100‑1 (MQ=2)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                  <  1:1962536/101‑1 (MQ=255)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                  >  1:3719664/1‑101 (MQ=255)
                                                                    gTCGCGCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                  <  1:4581766/101‑1 (MQ=2)
                                                                     tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                 <  1:3999578/101‑1 (MQ=255)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCCCGGGGGCGGGGATCGaaaaggggc                               >  1:884633/1‑95 (MQ=2)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                <  2:4330274/101‑2 (MQ=2)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGGGGATCGACACATcg                                >  2:2297310/1‑100 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                              <  2:851291/101‑7 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                               <  2:1866339/101‑1 (MQ=255)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgc                              >  1:3344850/1‑101 (MQ=255)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgc                              >  2:280529/1‑101 (MQ=255)
                                                                         tCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaag                              >  1:3495139/1‑96 (MQ=0)
                                                                         tCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgcg                             >  2:262560/1‑101 (MQ=255)
                                                                         cccccccGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaa                            <  2:3617962/100‑9 (MQ=2)
                                                                         cccccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaag                              <  1:2410790/100‑6 (MQ=0)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaag                           <  1:72623/101‑10 (MQ=2)
                                                                          ccccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatga                            >  1:4088146/1‑96 (MQ=2)
                                                                            ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatc                          <  1:2553121/101‑8 (MQ=2)
                                                                            ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatc                          <  1:2628110/101‑8 (MQ=2)
                                                                            ccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatc                          <  2:2188897/101‑8 (MQ=2)
                                                                             cccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGCTCGAAACATcgatgatca                         >  2:521830/1‑93 (MQ=0)
                                                                             cccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatca                         <  1:3424577/101‑9 (MQ=0)
                                                                             cccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATc                         <  1:2833394/101‑1 (MQ=255)
                                                                             cccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATc                         <  1:2766908/101‑1 (MQ=255)
                                                                             cccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATc                         <  2:4483735/101‑1 (MQ=255)
                                                                               cggggggtCGTGGATCGAAACAACCCTATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCgg                       <  1:437389/98‑1 (MQ=255)
                                                                               cgcgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCgg                       <  1:2486039/101‑1 (MQ=255)
                                                                               cgcgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCgg                       <  2:4185211/101‑1 (MQ=255)
                                                                                gcgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGc                      <  2:4115329/101‑1 (MQ=255)
                                                                                 cgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCg                     <  2:1828725/101‑1 (MQ=255)
                                                                                  gggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGt                    >  1:3048734/1‑101 (MQ=255)
                                                                                   ggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                   <  1:3593312/101‑1 (MQ=255)
                                                                                   ggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                   <  2:1106513/101‑1 (MQ=255)
                                                                                    gggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTgc                  >  2:2104938/1‑100 (MQ=255)
                                                                                    gggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGa                  <  1:632945/101‑1 (MQ=255)
                                                                                     ggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGACGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAc                 >  2:714505/1‑101 (MQ=255)
                                                                                      gCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcc                <  1:2873334/101‑1 (MQ=255)
                                                                                       cGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAccc               <  2:3457423/101‑1 (MQ=255)
                                                                                        gTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg              <  1:2132300/101‑1 (MQ=255)
                                                                                          ggATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcc            >  1:588676/1‑101 (MQ=255)
                                                                                           gATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGaagatcgg           >  1:3381590/1‑94 (MQ=255)
                                                                                            aTCGAACCAACCCTATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTTGATCGAAACATCCCCGGATCTGCGTGAcccgcccc          <  2:1224584/101‑1 (MQ=38)
                                                                                            aTCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc          >  2:3249656/1‑101 (MQ=255)
                                                                                                aaaCAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc      >  2:1474120/1‑101 (MQ=255)
                                                                                                   caacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc   >  2:134150/1‑101 (MQ=255)
                                                                                                   caacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc   >  1:4302263/1‑101 (MQ=255)
                                                                                                   caacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc   >  1:132750/1‑101 (MQ=255)
                                                                                                   caacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGGGACCCGCCCCGCCCATc   >  1:666989/1‑101 (MQ=255)
                                                                                                    aacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCg  >  1:352202/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
ATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCG  >  NC_014034/1397432‑1397632

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