Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,022 C→G 12.5% intergenic (+834/‑97) RCAP_rcc01306 → / → RCAP_rcc01307 Cas2 family CRISPR‑associated protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,0220C→G12.5% 99.0 0.0 ‑0.4339intergenic (+834/‑97)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (19/23):  ref base (134/162):  total (154/185)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.90e-01

Only 100 of 361 total aligned reads displayed.

GATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395922‑1396122
                                                                                                      |                                                                                                       
 aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACg                                                                                                        >  2:4832365/1‑101 (MQ=255)
   cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGcc                                                                                                      <  1:4301024/101‑1 (MQ=255)
   cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGcc                                                                                                      <  2:2056968/101‑1 (MQ=255)
    gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGccc                                                                                                     >  2:5213506/1‑100 (MQ=255)
    gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCt                                                                                                     >  2:4903752/1‑101 (MQ=255)
     aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTg                                                                                                    >  1:4166521/1‑101 (MQ=255)
        cGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.cg                                                                                                <  1:238145/101‑1 (MQ=255)
         gATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.cgc                                                                                               <  1:2132675/101‑1 (MQ=255)
           tGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..Gt                                                                                           <  2:3759686/101‑1 (MQ=255)
            gCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTc                                                                                          >  2:3604996/1‑101 (MQ=255)
             cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCg                                                                                         <  2:2095655/101‑1 (MQ=255)
                     gCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTTCC..TCACCCCTGC.CGC..GTCGCTcccccc                                                                                 >  1:5117047/1‑101 (MQ=255)
                      ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCcg                                                                                <  1:4464022/101‑1 (MQ=255)
                      ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCcg                                                                                >  2:3854736/1‑101 (MQ=255)
                       cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCcgc                                                                               >  2:3286711/1‑101 (MQ=255)
                             cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGc                                                                         >  1:3645473/1‑101 (MQ=255)
                              gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCg                                                                        >  2:748416/1‑101 (MQ=255)
                               atctgTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCCCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGt                                                                       <  1:1444162/97‑1 (MQ=255)
                                 cAGTCGAGTTGTCGCTCCCCCCGCGGGGGCGTGGACCGAAACTCGCCCCCGCTCCGCGCTGGCC..TCACCCCTGC.CGC..GTCTCTCCCCCCGCGGGGGCGTgg                                                                     <  2:3723268/101‑1 (MQ=25)
                                  aggggccTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAACCTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGa                                                                    <  1:1007181/95‑1 (MQ=38)
                                   gTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   <  2:1247863/101‑1 (MQ=255)
                                   gTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTAGAt                                                                   >  2:1444069/1‑101 (MQ=255)
                                    tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  <  1:5530928/101‑1 (MQ=255)
                                    tGGACGGGTCGCGCCCCCCGCGGGGTCGTGGATCTAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  <  2:870029/101‑1 (MQ=38)
                                     ttcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  <  1:1350890/96‑1 (MQ=9)
                                     tgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  <  2:2231940/96‑1 (MQ=2)
                                     ggcgcggTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  <  1:630716/96‑1 (MQ=2)
                                       aCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                               <  2:3102602/101‑1 (MQ=255)
                                       aCTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                               <  1:2493161/101‑1 (MQ=255)
                                       aCTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                               <  1:4371265/101‑1 (MQ=255)
                                        gcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                               <  2:5585802/99‑1 (MQ=2)
                                        cTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                              <  2:2434526/101‑1 (MQ=255)
                                         tGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCNGGGGCGTGGATCGAAAc                                                             >  2:2564574/1‑101 (MQ=255)
                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAg                                                           >  1:3410985/1‑101 (MQ=255)
                                            tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActt                                                           <  2:1959601/101‑3 (MQ=2)
                                            tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActt                                                           <  2:2321253/101‑3 (MQ=2)
                                            tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActt                                                           <  1:1981962/101‑3 (MQ=2)
                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgc                                                         <  1:474198/101‑1 (MQ=255)
                                               gcccccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCGGTCAGCCCCGTATTGCGCGC.TG...GTCGCTCCCCCCGCGGGGGCGGGGATCGAAActtgat                                                        >  1:4106998/3‑96 (MQ=2)
                                               cTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgcgc                                                       >  2:2703023/1‑101 (MQ=255)
                                               cTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgcgc                                                       >  1:5278945/1‑101 (MQ=255)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                       >  1:5663087/1‑95 (MQ=2)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                       >  1:4916508/1‑95 (MQ=2)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                                      <  1:4348258/101‑1 (MQ=255)
                                                 ccccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatga                                                      <  1:161132/101‑8 (MQ=2)
                                                 ccccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCC..TC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatga                                                      >  1:1508549/1‑94 (MQ=2)
                                                 ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                                     >  1:2998106/1‑101 (MQ=255)
                                                  cccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATc                                                    >  2:50999/1‑101 (MQ=255)
                                                   ccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCt                                                   >  1:797114/1‑101 (MQ=255)
                                                    cccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTg                                                  <  2:5357209/101‑1 (MQ=255)
                                                        cgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCT.CCCCCGCGGGGGCGTGGGTCGAAACAGGCGCATCTGTTGta                                             >  2:4005173/1‑101 (MQ=255)
                                                          ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                            <  1:457218/101‑1 (MQ=255)
                                                          ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                            >  1:5413820/1‑101 (MQ=255)
                                                           gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                           <  2:4209748/101‑1 (MQ=255)
                                                           gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                           <  2:3674332/101‑1 (MQ=255)
                                                            gatctGGATCGAAACTCCCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTatat                                          <  2:2926023/97‑1 (MQ=255)
                                                             gCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATg                                         <  1:323689/101‑1 (MQ=255)
                                                              cGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGc                                        >  2:4989369/1‑101 (MQ=255)
                                                              cGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGc                                        >  1:329171/1‑101 (MQ=255)
                                                              cGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGc                                        >  2:652216/1‑101 (MQ=255)
                                                                tGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGGCGCATCTGATGTATATGccc                                      >  2:4335156/1‑101 (MQ=255)
                                                                tGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGccc                                      >  2:5501937/1‑101 (MQ=255)
                                                                tGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGccc                                      <  1:3760206/101‑1 (MQ=255)
                                                                   aTCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGt                                   >  1:54862/1‑101 (MQ=255)
                                                                   aTCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGt                                   <  2:2396051/101‑1 (MQ=255)
                                                                    tCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTg                                  <  1:4832646/101‑1 (MQ=255)
                                                                      caaaCCCGCCCTCGATCCGCGCTTGCC..TCCCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                                <  1:1075237/100‑1 (MQ=255)
                                                                       aaaCTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgaga                               <  2:5784638/101‑1 (MQ=255)
                                                                       aaaCTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgaga                               <  2:3338963/101‑1 (MQ=255)
                                                                        aaCTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGAc                              >  1:3885628/1‑101 (MQ=255)
                                                                        aaCTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGAc                              <  1:1312119/101‑1 (MQ=255)
                                                                         aCTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACg                             >  1:1160518/1‑101 (MQ=255)
                                                                         aCTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACg                             >  1:2060887/1‑101 (MQ=255)
                                                                            cccccTCGATCCGCGCTTGCC..CCACGCCTGC.CGC..GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCt                          <  1:3078861/99‑1 (MQ=255)
                                                                            cGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCt                          <  2:1582880/101‑1 (MQ=255)
                                                                              cccTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGt                        <  1:1868232/101‑1 (MQ=255)
                                                                               ccTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTc                       <  2:5634478/101‑1 (MQ=255)
                                                                                cTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTca                      <  2:1642209/101‑2 (MQ=255)
                                                                                  cGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt                    <  1:2557373/101‑1 (MQ=255)
                                                                                  cGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt                    >  2:3070969/1‑101 (MQ=255)
                                                                                    aTCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcc                  <  1:2680595/101‑1 (MQ=255)
                                                                                     tCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccc                 >  2:3440071/1‑101 (MQ=255)
                                                                                      ccGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccc                >  2:4435855/1‑101 (MQ=255)
                                                                                       cgcgCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccccc               >  1:2752515/1‑101 (MQ=255)
                                                                                        gcgcTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc              <  1:888500/101‑1 (MQ=255)
                                                                                        gcgcTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc              <  2:5164008/101‑1 (MQ=255)
                                                                                        gcgcTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc              <  2:1729951/101‑1 (MQ=255)
                                                                                        gcgcTTGCC..TCACGCCTGC.CAC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc              <  2:4069357/101‑1 (MQ=255)
                                                                                         cgcTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcg             <  1:2900941/101‑1 (MQ=255)
                                                                                         cgcTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcg             <  1:4067869/101‑1 (MQ=255)
                                                                                           cTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcgcg           >  2:1763149/1‑101 (MQ=255)
                                                                                             tGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCggg         >  2:52321/1‑101 (MQ=255)
                                                                                             tGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGAGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCggg         >  2:1384590/1‑101 (MQ=255)
                                                                                                c..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGc      <  2:1074996/101‑1 (MQ=255)
                                                                                                c..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGc      >  2:1306994/1‑101 (MQ=255)
                                                                                                c..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGc      >  1:225422/1‑101 (MQ=255)
                                                                                                c..TCACGCCTGC.CGC..GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGc      <  1:3471539/101‑1 (MQ=255)
                                                                                                    cACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGt    >  1:5190397/1‑101 (MQ=255)
                                                                                                     ccGCCTGC.CGC..GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTg   <  2:5425361/100‑1 (MQ=255)
                                                                                                     aCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTg   <  1:5791677/101‑1 (MQ=255)
                                                                                                      |                                                                                                       
GATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCC..TCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395922‑1396122

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