Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,154 G→A 3.1% A12A (GCG→GTA RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,1540G→A 3.1% 9.1 ‑0.8 ‑0.6222A12A (GCG→GCARCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/5):  ref base (98/117):  total (100/122)
Fisher's exact test strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.76e-01

Only 100 of 291 total aligned reads displayed.

GGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTCCCCCCGCGGGGGC  >  NC_014034/1396054‑1396254
                                                                                                      |                                                                                                    
  aTCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGcgcc                                                                                                    <  2:150239/101‑1 (MQ=255)
   tCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGACACTTTTCGTC.ATCTGG.CAGGCTCGcgccg                                                                                                   >  2:1113026/1‑101 (MQ=255)
   tCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCgcgcgc                                                                                                   >  2:544457/1‑99 (MQ=255)
     gAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCgccgcc                                                                                                 >  1:1000279/1‑101 (MQ=255)
     gAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCgccgcc                                                                                                 >  1:5048238/1‑101 (MQ=255)
     gAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCCCgccgcc                                                                                                 >  2:3701388/1‑101 (MQ=255)
     gAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCCCgccgcc                                                                                                 >  2:1645187/1‑101 (MQ=255)
           ggCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTg                                                                                           <  2:3933573/101‑1 (MQ=255)
            gcgcATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCCCGCTGc                                                                                          >  2:4833388/1‑101 (MQ=255)
              gcATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGt                                                                                        <  2:4414796/101‑1 (MQ=255)
               cATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCCCCCTGCGTc                                                                                       >  2:4518433/1‑101 (MQ=255)
               cATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTc                                                                                       >  1:3977461/1‑101 (MQ=255)
               cATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTc                                                                                       <  1:3856607/101‑1 (MQ=255)
                 tCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGc                                                                                     <  1:994134/101‑1 (MQ=255)
                 tCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCCCCGCCACGCCGCggtgc                                                                                     >  1:4157421/1‑97 (MQ=255)
                  cTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCt                                                                                    >  1:856290/1‑101 (MQ=255)
                   tGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTc                                                                                   >  2:865227/1‑101 (MQ=255)
                        tataTGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTcccccc                                                                              <  1:2210708/101‑1 (MQ=255)
                        tataTGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTcccccc                                                                              <  2:3664271/101‑1 (MQ=255)
                        tataTGCACGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTcccccc                                                                              <  1:4266241/101‑1 (MQ=255)
                         atatGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTccccccc                                                                             >  1:3770631/1‑100 (MQ=255)
                           atGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACCCTGCGCCGCTCCCCCcgcg                                                                           >  1:1373032/1‑101 (MQ=255)
                            tGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCgg                                                                          >  1:2013274/1‑101 (MQ=255)
                              cccGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCgggg                                                                        <  2:2479700/101‑1 (MQ=255)
                              cccGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCgggg                                                                        <  2:2059802/101‑1 (MQ=255)
                               ccGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCggggg                                                                       >  2:3719429/1‑101 (MQ=255)
                                cGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGc                                                                      <  2:2107058/101‑1 (MQ=255)
                                cGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGAGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGGGTCGCACCCCCCCCGGGGGc                                                                      >  1:651189/1‑101 (MQ=38)
                                  gTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCgg                                                                    >  1:3400467/1‑100 (MQ=255)
                                  gTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGt                                                                    <  1:3712601/101‑1 (MQ=255)
                                   tgtctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCGGTCAGCCTCGTATTG.CGCGCTG.GTCGCTCCCCCCGCGGGGGCGTg                                                                   <  2:4218159/91‑1 (MQ=9)
                                   tgtctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCGGTCAGCCTCGTATTG.CGCGCTG.GTCGCTCCCCCCGCGGGGGCGTg                                                                   <  1:939925/91‑1 (MQ=9)
                                   tgtctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCGGTCAGCCTCGTATTG.CGCGCTG.GTCGCTCCCCCCGCGGGGGCGTg                                                                   >  1:3125486/11‑101 (MQ=9)
                                    ttctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCGGTCAGCCTCGTATTG.CGCGCTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                  <  2:4080269/92‑1 (MQ=2)
                                    gtctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCGGTCAGCCTCGTATTG.CGCGCTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                  <  2:3446395/92‑1 (MQ=2)
                                    gagaCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCCCGCTGCGTCGCTCCCCCCCCGGGGGCGGgg                                                                  >  1:4461839/1‑101 (MQ=255)
                                    gagaCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTgg                                                                  <  1:3828455/101‑1 (MQ=255)
                                    gagaCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTgg                                                                  >  2:4659187/1‑101 (MQ=255)
                                     agaCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCCCGGGGGCGTggg                                                                 >  1:3178220/1‑100 (MQ=255)
                                      gaCGGCTGACGCTCTCCCCGCGGGGGCGTGGATGGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                <  1:4363601/101‑1 (MQ=255)
                                       gccgCTTTCGCTCCCCCCGCGGGGGCGTGGATCGAACCTTTGCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCCTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  2:3026100/98‑1 (MQ=38)
                                       aCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  2:637646/101‑1 (MQ=255)
                                       aCGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  2:1218781/101‑1 (MQ=255)
                                        cGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              >  2:5046488/1‑101 (MQ=255)
                                        cGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              <  2:2558721/101‑1 (MQ=255)
                                            tGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          <  2:1751270/101‑1 (MQ=255)
                                              tCGCTACCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAccc                                                        >  2:35579/1‑101 (MQ=255)
                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACCTTTCGTC.ATCTGG.CAGGCTCGGGCCGCCACGCTGCGTCGCCCCCCCCGCGGGGGGGGCGCTCGAAACCCgg                                                      >  2:2336487/1‑101 (MQ=37)
                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCCCGGTGGCTGCCTCCCCCCGCGGGGGGGGGGACCATAACCCAGc                                                     >  2:4578890/1‑101 (MQ=21)
                                                  tCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCg                                                    >  1:2325966/1‑101 (MQ=255)
                                                   ccccccGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGa                                                   <  1:2674245/101‑1 (MQ=255)
                                                    cccccGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGaa                                                  >  2:3217798/1‑101 (MQ=255)
                                                      cccGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGaaaa                                                <  2:911137/101‑1 (MQ=255)
                                                      cccGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGaaaa                                                <  1:4792169/101‑1 (MQ=255)
                                                       ccGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGaaaaa                                               <  1:3650592/101‑1 (MQ=255)
                                                        cgcgGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAAc                                              <  1:250222/101‑1 (MQ=255)
                                                         gcgGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACg                                             <  2:4976597/101‑1 (MQ=255)
                                                         gcgGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACg                                             >  1:408899/1‑101 (MQ=255)
                                                           gggggCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGc                                           >  1:3019027/1‑101 (MQ=255)
                                                            ggggCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAAAGGcc                                          >  2:2245377/1‑101 (MQ=255)
                                                             gggCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCa                                         <  1:4423984/101‑1 (MQ=255)
                                                             gggCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCa                                         <  2:29925/101‑1 (MQ=255)
                                                             gggCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCa                                         >  2:4308628/1‑101 (MQ=255)
                                                               gCGGGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATg                                       <  1:556978/101‑1 (MQ=255)
                                                                cGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAACCACGGCCATGc                                      >  2:2465263/1‑101 (MQ=255)
                                                                cGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGc                                      <  2:2932260/101‑1 (MQ=255)
                                                                 gTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGcc                                     <  2:188795/101‑1 (MQ=255)
                                                                  tGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGccc                                    >  1:3251092/1‑100 (MQ=255)
                                                                  tGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCt                                    <  2:2767448/101‑1 (MQ=255)
                                                                  tGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAACCGGCGAAAAACGGCCATGccc                                    >  2:1784672/1‑100 (MQ=255)
                                                                   ggATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCCTGCCTc                                   >  2:3742562/1‑101 (MQ=255)
                                                                   ggATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTc                                   >  1:2858925/1‑101 (MQ=255)
                                                                    gATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAATGGCCATGCCTcc                                  <  1:2746338/101‑1 (MQ=255)
                                                                    gATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTcc                                  >  1:1305065/1‑101 (MQ=255)
                                                                     aTCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCg                                 <  2:82691/101‑1 (MQ=255)
                                                                      tCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCgg                                <  1:754927/101‑1 (MQ=255)
                                                                        gAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCgggg                              >  1:4487926/1‑100 (MQ=255)
                                                                        gAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCcggg                              >  2:4880280/1‑100 (MQ=255)
                                                                             ttttCGTT.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTg                         <  1:1113140/101‑1 (MQ=255)
                                                                              tttCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTgg                        >  1:103811/1‑101 (MQ=255)
                                                                               ttCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGc                       >  1:2761999/1‑101 (MQ=255)
                                                                                tCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCt                      >  1:3336901/1‑101 (MQ=255)
                                                                                 cGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTg                     <  1:4659510/101‑1 (MQ=255)
                                                                                  gTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGt                    <  2:734053/101‑1 (MQ=255)
                                                                                   tC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTc                   >  2:2711821/1‑101 (MQ=255)
                                                                                   tC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTc                   <  2:1927243/101‑1 (MQ=255)
                                                                                    c.ATCTGG.CAGGCGCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCg                  <  1:1685185/101‑1 (MQ=255)
                                                                                      aTCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGc                 >  1:4678377/1‑101 (MQ=255)
                                                                                        cTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCCCTc               >  2:3116294/1‑101 (MQ=255)
                                                                                        cTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGGGTCCGGCTGTCGCTc               >  1:1711377/1‑101 (MQ=255)
                                                                                         tGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCAGGGTGTCTGGCTGTCGCTcc              <  1:3455641/101‑1 (MQ=255)
                                                                                             cAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTccccc           >  1:500605/1‑101 (MQ=255)
                                                                                             cAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTccccc           <  1:3683922/101‑1 (MQ=255)
                                                                                              aGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTcccccc          <  1:1168801/101‑1 (MQ=255)
                                                                                              aGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGAGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTcccccc          >  2:2580742/1‑101 (MQ=255)
                                                                                               ggCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTCCCCCcg         >  1:1839494/1‑101 (MQ=255)
                                                                                                gCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTCCCCCcgc        >  2:3317709/1‑101 (MQ=255)
                                                                                                   cgcgCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTCCCCCCGCggg     >  1:1781379/1‑101 (MQ=255)
                                                                                                     cgccgcCACGCTGCGTCGCTCCCCCTGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTCCCCCCGCggggg   >  2:4151125/1‑101 (MQ=255)
                                                                                                     cgccgcCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTCCCCCCGCggggg   >  1:4152734/1‑101 (MQ=255)
                                                                                                      |                                                                                                    
GGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTTTCGTC.ATCTGG.CAGGCTCGCGCCGCCACGCTGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGGCGAAAAACGGCCATGCCTCCGGGTGTCTGGCTGTCGCTCCCCCCGCGGGGGC  >  NC_014034/1396054‑1396254

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG