Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,015 G→C 6.2% T214T (ACC→ACG RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0150G→C 6.2% 202.2 ‑0.5 ‑0.01792T214T (ACC→ACGRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (56/62):  ref base (833/836):  total (892/900)
Fisher's exact test strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 2295 total aligned reads displayed.

CGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGG  >  NC_014034/221915‑222115
                                                                                                    |                                                                                                    
ccgtgggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcg                                                                                                      >  2:4363393‑M1/42‑101 (MQ=255)
 cgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                                                                     >  1:1455811‑M1/41‑101 (MQ=255)
  gtggggtccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                                                    >  2:3157033‑M1/40‑101 (MQ=255)
  ctggcatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                                    >  2:1278681‑M1/40‑101 (MQ=255)
   gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGTTCGGGCGGGTCCGGGGGATCGGTCGGATCCGTGGGAttggtc                                                                                                   >  2:2276944/1‑101 (MQ=25)
           ggggaatctttCGGGTCCGGGGGATCGGTCGTATCCTTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGt                                                                                           <  2:4109081/91‑1 (MQ=2)
             cggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                                                         <  1:599320‑M1/73‑1 (MQ=255)
               tGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggg                                                                                       >  2:4227518/1‑100 (MQ=255)
                 ttcgtggggtcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                     >  2:4493937‑M1/25‑101 (MQ=255)
                  ccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtcg                                                                                    <  2:472415‑M1/78‑1 (MQ=255)
                     tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                 >  1:4586585/1‑101 (MQ=255)
                       ggatcggtggggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                               >  1:1095037‑M1/19‑101 (MQ=255)
                       ggATGGGTGGGATCCGTGGGGTCGGTCGGGTGTGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                               <  1:5496234/101‑1 (MQ=11)
                       ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                               >  1:2713129/1‑101 (MQ=255)
                           gCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                                                           <  1:1202738/101‑1 (MQ=255)
                           cggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                                           <  2:1043393‑M1/87‑1 (MQ=255)
                             gtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGt                                                                         >  1:4385847‑M1/13‑101 (MQ=255)
                               gggATCGGTCTGGGCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                                                       <  1:5616936/101‑3 (MQ=11)
                                ggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                                                      <  1:2493687‑M1/92‑1 (MQ=255)
                                ggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                                                      <  1:1554660/101‑1 (MQ=255)
                                ggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                                                      >  2:3366721/1‑101 (MQ=255)
                                 gatcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCggcggg                                                                     >  2:699247/3‑96 (MQ=11)
                                  gtcggttggGTCGGGCGGGTCCGGCGGGTCCGTCGGGGCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                                                    <  2:1678628/100‑1 (MQ=11)
                                   gccgtcggATCGGTGGGATCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                   <  1:4306615/98‑1 (MQ=11)
                                    gcgttgggtgggtcgggtccgtcgggtaggtcgggcccgtgggatcCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                  <  2:2862227/91‑1 (MQ=11)
                                    ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                                                  >  1:3840404‑M1/6‑101 (MQ=255)
                                      gtgGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                >  1:4035672‑M1/4‑101 (MQ=255)
                                      gtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggt                                                                <  1:423464‑M1/98‑1 (MQ=255)
                                      gtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATAGGTCGGGTCCGTGGGATCGGTCGGatccgtggggtccgt                                                                <  1:2261373‑M1/98‑7 (MQ=255)
                                      gtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAccggg                                                                >  2:2181156/1‑100 (MQ=11)
                                       tggggtccgtcggatccgtgggatcggtcggatccgtgggatcgGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAcaggtc                                                               >  2:1875737/9‑101 (MQ=11)
                                       tgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCGGTgggggccgtg                                                               >  2:5679725/3‑94 (MQ=11)
                                        gggatcggtcgggtccgtggagtcggtcgggtccgtgGGATCGGTCGGTTCCGTGGGATCGGTCGGTTCGGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                              <  1:3254806/100‑1 (MQ=11)
                                        gGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                              >  2:2997595‑M1/2‑101 (MQ=255)
                                         ggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGAGGGATCGGTCGGATCCGTgggcctggtggg                                                             >  2:2717969/1‑92 (MQ=14)
                                          gatcggttggatccgtggggtcggtcgggtccgtggGGTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgga                                                            <  2:1579820/101‑2 (MQ=11)
                                           gtcggtcggctccgcgggatcggtcgtgtccgtgcgaTCGGCGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                                                           <  2:5242925/100‑1 (MQ=14)
                                           agcgggcgggtctggggggtcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                                                           <  1:2031839/99‑1 (MQ=11)
                                             cggtgggatgcgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                         <  1:3604223/101‑1 (MQ=11)
                                             cggtcgggtccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatcc                                                         <  1:3734374‑M2/101‑5 (MQ=255)
                                             cggtcgggtccgtcggatcggtcggatccgtaggatCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcagatcggaag                                                         >  1:1291522/1‑91 (MQ=11)
                                              ggtcgggtccgtgggatcggtcgggtccgtgggatgGGTCGGGCCCGGGGGATCGGGCGGGGCCCTGGGATGGGTCGGGTCCGGGGGGTCGGTGGGGTCCg                                                        >  1:2611772/1‑101 (MQ=11)
                                              ggtcgggtccgtgggatcggtcggatccgtgggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCg                                                        <  1:2836426/101‑1 (MQ=11)
                                               gtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTTGGTTCCGAGGGATCGGTCGGGTCCgg                                                       >  2:4837497/1‑100 (MQ=11)
                                               gtcggatccgtgggatcggtcgggtccgtgggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgt                                                       >  2:3262772‑M2/1‑95 (MQ=255)
                                               ggcgggtccgggggagcggtcgggtccgtgggatcggtGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                       <  1:827932/99‑1 (MQ=11)
                                                tcgcgtccgtgggatcggtcgggtccgggggatcggTCGGGTCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtg                                                      <  2:5134755/101‑2 (MQ=11)
                                                 ctggtccttgggatcggtcgggtccgttggatcggtcgGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtgg                                                     <  2:3092143/99‑3 (MQ=11)
                                                 cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgg                                                     >  2:1487822‑M2/1‑93 (MQ=255)
                                                 cggatccgtgggatcggtcggatccgtgggatcggtCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgagatcggaa                                                     >  2:802990/1‑92 (MQ=14)
                                                   ggtccgtgggatcggtcgggtccgtcggatcggtcgGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGCCCGTGGGa                                                   >  1:2147917/1‑101 (MQ=14)
                                                   ggtccgtgggatcgctcgtgtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                                   <  2:2397799/101‑1 (MQ=11)
                                                    tcttccgatctcggtcggagccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                  <  1:4113566/91‑1 (MQ=14)
                                                    gaccgtgggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAt                                                  <  2:4586309/99‑1 (MQ=11)
                                                    atccgtgggatcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggac                                                  <  1:3480234‑M2/100‑12 (MQ=255)
                                                     tccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatc                                                 <  1:2345208‑M2/101‑13 (MQ=255)
                                                     tccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatc                                                 >  1:3153473‑M2/1‑89 (MQ=255)
                                                     tccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGCGGGAtg                                                 >  2:4305572/1‑100 (MQ=11)
                                                     tccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                 <  1:1644365/101‑1 (MQ=14)
                                                       cgtggggtcggtcgggtccgtggggtccgtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                               <  1:5005178/101‑1 (MQ=11)
                                                        gtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                              <  1:4671557/101‑1 (MQ=11)
                                                         tgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGACGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                             >  2:3840159/1‑101 (MQ=14)
                                                         tgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCGGTGGGGtcggtg                                             >  2:2066709/1‑100 (MQ=14)
                                                         tgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtc                                             <  2:5363083‑M2/101‑17 (MQ=255)
                                                         tgggatcggtcggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGAtcgggc                                             >  1:3952069/1‑99 (MQ=11)
                                                           ggatcggtcgggtccgtggggtgcgtaggAGCCGGGGGATCGGTCGGATCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                           <  1:715557/101‑1 (MQ=11)
                                                           ggatcgggcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                           <  1:3080770/101‑1 (MQ=11)
                                                            gatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggg                                          >  1:1832947/1‑101 (MQ=11)
                                                             atcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggt                                         <  1:2931929‑M2/101‑21 (MQ=255)
                                                               cgggcgggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatcc                                       <  1:1362309/101‑1 (MQ=11)
                                                                ggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccg                                      >  1:5708758/1‑101 (MQ=11)
                                                                 gtcggatccgtgggttccgtcggACCCGTGGGATCGCTCGGATTCTGGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgt                                     <  2:582978/101‑1 (MQ=11)
                                                                  tcgggtccgttgggtcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtg                                    >  2:5455901‑M2/1‑76 (MQ=255)
                                                                  tcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGAACGgacgggtccggg                                    >  2:3603995/1‑99 (MQ=14)
                                                                   cggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGgtgggggcggtgg                                   >  2:4041900/1‑94 (MQ=14)
                                                                    gggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtagg                                  >  2:5556476‑M2/1‑74 (MQ=255)
                                                                     ggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtggga                                 >  2:3943540‑M2/1‑73 (MQ=255)
                                                                     ggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtggga                                 <  1:5195530‑M2/101‑29 (MQ=255)
                                                                      atccgtgggatcggtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtaggatcggtcgggtccgtaggat                                <  2:4929274‑M2/100‑30 (MQ=255)
                                                                      atccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggt                                <  2:1775045/100‑1 (MQ=11)
                                                                        ccgtgggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                              <  1:100380/101‑1 (MQ=14)
                                                                            gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcg                          <  1:5296944‑M2/101‑36 (MQ=255)
                                                                             gggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTggg                         <  1:490359/95‑1 (MQ=11)
                                                                             ggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggg                         <  1:3366926/101‑1 (MQ=11)
                                                                              ggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGa                        >  1:93188/6‑101 (MQ=11)
                                                                               gtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAt                       <  1:5097013/97‑1 (MQ=11)
                                                                                   gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggt                   <  1:1712459/101‑1 (MQ=14)
                                                                                   gtaggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgt                   <  2:5037147‑M2/101‑43 (MQ=255)
                                                                                     cggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                 >  1:4804995‑M2/1‑57 (MQ=255)
                                                                                      gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgg                >  1:5404887/4‑101 (MQ=14)
                                                                                        gtcCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggat              <  1:195957‑M2/100‑48 (MQ=255)
                                                                                        gtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgggg              >  2:852290/2‑99 (MQ=14)
                                                                                         tcggTGGGTTCGTTCGGGTCCTTGGGATCTGTCGGGTCCGGGGTGGCGGTCGGGTCGGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTc             <  2:1979854/98‑1 (MQ=9)
                                                                                         gccgTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc             >  2:372677‑M2/2‑53 (MQ=255)
                                                                                          ccGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcc            >  2:1155708‑M2/1‑52 (MQ=255)
                                                                                             tGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtc         <  1:450812‑M2/101‑53 (MQ=255)
                                                                                               ggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggaccgtgggatctgtcgg       <  2:2545071‑M2/101‑55 (MQ=255)
                                                                                                  tcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtc    <  2:1341548‑M2/101‑58 (MQ=255)
                                                                                                   cggtctgGTCCGTGGGTTCGGTCGCCTCCGAGGTCTCGGCTCGTTCCCTGCGCTCGGTCGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg   <  2:2645015/101‑1 (MQ=2)
                                                                                                    ggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg  <  2:3477533‑M2/101‑60 (MQ=255)
                                                                                                    |                                                                                                    
CGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGG  >  NC_014034/221915‑222115

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG