Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,493,971 Δ1 bp 1.4% intergenic (+19/‑223) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,493,9710G→. 1.4% 5.8 ‑0.5 ‑0.0365intergenic (+19/‑223)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (91/269):  total (93/272)
Fisher's exact test strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.11e-01

Only 100 of 395 total aligned reads displayed.

TGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GG.GGGGG..CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAGCAGA  >  NC_014034/2493876‑2494067
                                                                                                   |                                                                                                
tGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..AGCTGCCCCCCGTCCg                                                                                              >  2:790557/1‑101 (MQ=255)
  aTGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCa                                                                                            <  2:4613818/101‑1 (MQ=39)
    ggCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAg                                                                                          >  2:772440/1‑101 (MQ=38)
     gCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGc                                                                                         <  1:342539/101‑1 (MQ=38)
       cGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..GGCTGGCCCCCGTCCGCAAGCgg                                                                                       >  2:1425922/1‑101 (MQ=38)
         ggCTGCCCCCCCGGGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGGAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGAc                                                                                     <  1:8626689/100‑1 (MQ=25)
         aGCTGCCCGCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGAc                                                                                     <  1:4708498/101‑1 (MQ=37)
           cTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTc                                                                                   <  2:3056340/101‑1 (MQ=37)
           cTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTc                                                                                   <  1:4477869/101‑1 (MQ=37)
            tGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTcc                                                                                  <  1:2649921/101‑1 (MQ=37)
             gCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCCAGCGGACTccc                                                                                 >  1:7627945/1‑101 (MQ=37)
             gCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTccc                                                                                 >  2:2715102/1‑101 (MQ=37)
               ccACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCCAGCGGACTccccc                                                                               >  2:1154309/1‑101 (MQ=37)
               ccACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTccccc                                                                               >  2:2583795/1‑101 (MQ=37)
                 ccccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCg                                                                             <  1:2844268/100‑1 (MQ=25)
                 ccccGCGCCCCGCGTCGGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCg                                                                             <  1:2394515/100‑1 (MQ=25)
                   ccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCggg                                                                           >  2:3775742/1‑101 (MQ=25)
                    cgcggcccgcTTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGAAAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGa                                                                          <  1:2595437/95‑1 (MQ=21)
                     gcgcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGat                                                                         <  2:5511536/101‑1 (MQ=25)
                     gcgcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGat                                                                         <  1:2135915/101‑1 (MQ=25)
                     gcgcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGat                                                                         <  1:3505917/101‑1 (MQ=25)
                     gcgcgcCGCGTCGGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGat                                                                         <  2:2831014/101‑1 (MQ=25)
                       gcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGatat                                                                       <  2:1010938/101‑1 (MQ=25)
                       gcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACAACAGGCCGAGTAAGGG.GGGGG..GGCTGCCCCCCGGCCGCAAGCGGACTCCCCCCGGGgtat                                                                       >  2:6489274/1‑101 (MQ=21)
                        cgccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCCAGCGGACTCCCCCCCGGATAtt                                                                      >  2:8497676/1‑101 (MQ=25)
                        cgccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATAtt                                                                      <  1:3747971/101‑1 (MQ=25)
                         gccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATAttt                                                                     <  1:7090858/101‑1 (MQ=25)
                         gccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATAttt                                                                     <  2:5002641/101‑1 (MQ=25)
                          ccccgTCTGGCCCAGCGCCGCGGCGGGCCCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTg                                                                    <  2:4210597/98‑1 (MQ=21)
                             cgTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcag                                                                 <  2:5348972/101‑1 (MQ=25)
                             cgTATGGCCCAGCGCCGCGGCGGGCGGGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcag                                                                 <  2:5752692/101‑1 (MQ=21)
                              gTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGGCGCCCCCCCTCCGCCAGCGGACTCCCCCCGGGGTATTTgcagc                                                                >  1:2222050/1‑101 (MQ=25)
                               tCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                               <  2:8910033/101‑1 (MQ=25)
                                cTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCaa                                                              >  2:3597667/1‑101 (MQ=25)
                                cTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCaa                                                              >  2:6027240/1‑101 (MQ=25)
                                  ggCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGa                                                            <  2:2872178/101‑1 (MQ=21)
                                   gCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATAATTGCAGCAAGAt                                                           <  1:6573133/101‑1 (MQ=21)
                                     ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                                         <  2:8083763/101‑1 (MQ=21)
                                     ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                                         <  2:4697996/101‑1 (MQ=21)
                                      cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaa                                                        >  1:4257417/1‑101 (MQ=21)
                                      cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGTGGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATTaa                                                        <  2:9146321/101‑1 (MQ=21)
                                       aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                       <  1:8468030/101‑1 (MQ=21)
                                       aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                       <  2:8838211/101‑1 (MQ=21)
                                       aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                       >  1:1291021/1‑101 (MQ=21)
                                        gcgcCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAc                                                      <  1:5079510/101‑1 (MQ=21)
                                         cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                     <  1:6872812/101‑1 (MQ=21)
                                         cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                     <  2:5256757/101‑1 (MQ=21)
                                         cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                     <  1:6363469/101‑1 (MQ=21)
                                          gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                    <  2:7460995/101‑1 (MQ=21)
                                          gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                    >  1:7704678/1‑101 (MQ=21)
                                           ccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCgg                                                   <  2:4003316/101‑1 (MQ=21)
                                           ccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCgg                                                   <  1:935108/101‑1 (MQ=21)
                                             gcggcgGGGGCGCGGCACGACAGGCCGGGTAAGGG.GGGGG..CGCTGCGCCCCCTCCCGAAGCGGCCTTCCCCCGGGGAAGTTGCAGGACGATCAAACCGGag                                                 >  2:162524/1‑100 (MQ=16)
                                              cggcggGCGCGCGGCACGACAGGCCGGGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGGAAGATGAAACCGGATg                                                >  2:5961040/1‑101 (MQ=18)
                                              cggcggGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATg                                                <  2:8683934/101‑1 (MQ=21)
                                              cggcggGCGCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATg                                                <  2:2095322/101‑1 (MQ=21)
                                               ggcgggcgCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGGGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGc                                               >  1:5295295/1‑101 (MQ=21)
                                                  gggcgCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                            >  1:5995945/1‑101 (MQ=21)
                                                   ggcgCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                           <  1:1945162/101‑1 (MQ=21)
                                                   ggcgCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                           <  1:5678555/101‑1 (MQ=21)
                                                   ggcgCGCGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                           <  2:6888974/101‑1 (MQ=21)
                                                    gcgcgcGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAGATGAAACCGGCTGCCTTTc                                          >  1:1802607/1‑101 (MQ=18)
                                                    gcgcgcGGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                          <  1:786711/101‑1 (MQ=21)
                                                     cgcgcgGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGCTATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                         <  1:6468131/101‑1 (MQ=21)
                                                     cgcgcgGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                         >  1:4839243/1‑101 (MQ=21)
                                                     cgcgcgGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCa                                         <  2:3501035/101‑1 (MQ=21)
                                                         cgGCACGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                     <  2:6857187/101‑1 (MQ=21)
                                                             aCGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCCTCTTGccc                                 >  2:1895580/1‑101 (MQ=18)
                                                             aCGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                 >  1:314690/1‑101 (MQ=21)
                                                              cGACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                <  2:3357381/101‑1 (MQ=21)
                                                               gACAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                               >  2:7073967/1‑101 (MQ=21)
                                                                aCAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaa                              <  2:379302/101‑1 (MQ=21)
                                                                aCAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaa                              <  1:3255690/101‑1 (MQ=21)
                                                                 cAGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                             >  2:7990327/1‑101 (MQ=21)
                                                                 cAGGCCGAGTAA.GG.GGGGGT.CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                             <  1:3587970/101‑1 (MQ=21)
                                                                  aGGCCGAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGCAat                            >  2:6214733/1‑101 (MQ=18)
                                                                    gCCGGGTAA.GG.GGGGGGCCCCTCCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                           <  1:8946637/101‑1 (MQ=18)
                                                                       gAGTAAGGG.GGGGG..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                       >  1:1152592/1‑101 (MQ=18)
                                                                          taaggg.ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                    >  2:7402294/5‑101 (MQ=255)
                                                                             ggg.ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                 <  2:6346311/100‑1 (MQ=255)
                                                                             ggg.ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGAGGAAACCGGATGCCTTTCATCTTGCCCGAAATATCACggcggg                 >  1:8301384/2‑101 (MQ=18)
                                                                             ggg.ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCACGGGCTATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                 >  1:4038957/2‑101 (MQ=18)
                                                                              gg.ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggggg                >  1:5378965/1‑100 (MQ=255)
                                                                              gg.ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggggg                >  2:4399036/1‑100 (MQ=255)
                                                                               g.ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg               <  1:880170/101‑1 (MQ=255)
                                                                                 ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa              <  1:1052772/101‑1 (MQ=255)
                                                                                 ggggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATAGCCCGGGGGGTGa              <  1:8155471/101‑1 (MQ=255)
                                                                                   ggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAAt            <  2:7193470/101‑1 (MQ=255)
                                                                                   ggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAAt            >  2:5189108/1‑101 (MQ=255)
                                                                                   ggg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAAt            <  2:8176758/101‑1 (MQ=255)
                                                                                    gg..CGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATg           <  2:3220547/101‑1 (MQ=255)
                                                                                        cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGc         >  1:4154401/1‑101 (MQ=255)
                                                                                        cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGc         >  1:734792/1‑101 (MQ=255)
                                                                                        cgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAATGGc         >  1:6159433/1‑101 (MQ=255)
                                                                                         gcTGCCCCCCGTCCGCAAGCGGACTCCCCGCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCg        <  2:2362552/101‑1 (MQ=255)
                                                                                            gCCCCCCGTCCGCAAGCGGAGTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGCGGGTGAATGGCGaag     <  1:7251021/101‑3 (MQ=21)
                                                                                            gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGaag     >  1:5765528/1‑99 (MQ=255)
                                                                                             cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGaagcc   <  1:9107171/101‑5 (MQ=255)
                                                                                             cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGaagcc   <  1:8102985/101‑5 (MQ=255)
                                                                                              ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGaagcca  >  2:5877657/1‑96 (MQ=255)
                                                                                                   |                                                                                                
TGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GG.GGGGG..CGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAGCAGA  >  NC_014034/2493876‑2494067

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