Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,961 G→C 3.1% T232T (ACC→ACG RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9610G→C 3.1% 26.9 ‑2.4 ‑1.2542T232T (ACC→ACGRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/5):  ref base (255/267):  total (269/273)
Fisher's exact test strand distribution p-value = 8.85e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.90e-02

Only 100 of 1205 total aligned reads displayed.

CGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCG  >  NC_014034/221861‑222061
                                                                                                    |                                                                                                    
ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtcggtgggatcggtcggatccgggggatcggtcgggtccgtgggatcg                                                                                                      <  1:3718662‑M1/6‑1 (MQ=255)
ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                                                      >  1:926703‑M1/96‑101 (MQ=255)
cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                                      >  2:3366552/1‑101 (MQ=255)
 ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                                     >  2:3527668/1‑101 (MQ=255)
  gggggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggt                                                                                                    >  1:1773281‑M1/94‑101 (MQ=255)
   tgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                                   <  2:2839475‑M1/9‑1 (MQ=255)
    gggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg                                                                                                  <  1:3003714‑M1/10‑1 (MQ=255)
         cggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtcc                                                                                             >  1:3855724‑M1/87‑101 (MQ=255)
           gtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgg                                                                                           <  2:1285569‑M1/17‑2 (MQ=255)
           acCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGGGGGCGTCGGAGGCGTGGGCTCGGTCggatcggtcgggtccgt                                                                                           >  1:1353912/1‑101 (MQ=255)
            ccGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtg                                                                                          <  1:3585628/101‑1 (MQ=255)
               ggtcggtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggga                                                                                       >  1:4015616‑M1/81‑101 (MQ=255)
                   cGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgggaggggacgtgggatcgg                                                                                   >  1:1077975/1‑101 (MQ=255)
                     aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                                 >  1:2641239/1‑101 (MQ=255)
                     aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgggggctgggtc                                                                                 >  2:4091209/1‑101 (MQ=255)
                      gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggagcggtcg                                                                                <  2:635753‑M1/28‑1 (MQ=255)
                      gggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcg                                                                                >  1:773111‑M1/74‑101 (MQ=255)
                       ggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgg                                                                               <  2:3457229‑M1/29‑1 (MQ=255)
                         ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                             <  2:3080498/101‑1 (MQ=255)
                         atcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggat                                                                             <  1:2512305‑M1/31‑3 (MQ=255)
                         atcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgagggatcgggcgggg                                                                             <  2:926618‑M1/31‑2 (MQ=255)
                          tcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtc                                                                            <  1:2388599‑M1/32‑1 (MQ=255)
                          tcgggcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtc                                                                            >  1:4083713‑M1/70‑101 (MQ=255)
                             gtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgt                                                                         >  1:2757070‑M1/67‑101 (MQ=255)
                              tcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtg                                                                        <  1:2590032‑M1/36‑1 (MQ=255)
                                  atccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggat                                                                    <  1:3419573‑M1/40‑1 (MQ=255)
                                    tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcggggccgtgggatcg                                                                  >  1:2708375/1‑101 (MQ=255)
                                    ccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtggggtcg                                                                  <  1:3148560‑M1/42‑1 (MQ=255)
                                     cgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcgg                                                                 <  1:1226238‑M1/43‑1 (MQ=255)
                                      gtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgt                                                                <  1:2121169‑M1/44‑7 (MQ=255)
                                        cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                                              <  2:2646024/101‑1 (MQ=255)
                                        cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcggggccgtgggatcggtcg                                                              >  1:1844007/1‑101 (MQ=255)
                                         gggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgg                                                             >  1:4165302‑M1/55‑101 (MQ=255)
                                          ggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa                                                            <  1:3284898‑M1/48‑1 (MQ=255)
                                           gtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                                           >  2:812104‑M1/53‑101 (MQ=255)
                                            gTCGCCGGATCGCCCTCATCGACGGTGGGCGTCTGATGCGTGGGCCCGGTCGGCTCGGCGGGGACCGTGGGAGCGGTCGGGTCCGTGGGCCCGGTGGGGTc                                                          >  2:2105700/1‑101 (MQ=255)
                                             cggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                         >  1:3783316‑M1/51‑101 (MQ=255)
                                               gtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGt                                                       <  1:3926953‑M1/53‑1 (MQ=255)
                                                tcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                                                      <  1:1588498‑M1/54‑1 (MQ=255)
                                                tcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                      <  1:1720877‑M1/54‑1 (MQ=255)
                                                 cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                     >  2:2450178/1‑101 (MQ=255)
                                                  gggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTggg                                                    <  2:2005778‑M1/56‑1 (MQ=255)
                                                  gggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTggg                                                    <  2:2994505‑M1/56‑1 (MQ=255)
                                                   ggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGa                                                   <  1:1005158‑M1/57‑1 (MQ=255)
                                                      ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                                >  1:2807912‑M1/42‑101 (MQ=255)
                                                      ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGCtcg                                                <  1:3099519‑M1/60‑1 (MQ=255)
                                                      ccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGATCCCTGGGAtcg                                                <  1:1445810‑M1/60‑1 (MQ=255)
                                                       cgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcgg                                               <  1:633119‑M1/61‑1 (MQ=255)
                                                       cgtgggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcgg                                               <  1:2266379‑M1/61‑1 (MQ=255)
                                                         tgggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                             >  2:3140854‑M1/39‑101 (MQ=255)
                                                               cgtgcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                       >  1:164226‑M1/33‑101 (MQ=255)
                                                               cggtcgggtccgtcggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                       <  1:114800‑M1/69‑1 (MQ=255)
                                                               ccgtcggatccgtgggatcggtgggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                       >  1:2912483‑M1/33‑101 (MQ=255)
                                                                 gtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatcctt                                     <  2:464862‑M1/71‑3 (MQ=255)
                                                                  tgggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                                    >  1:2314384‑M1/30‑101 (MQ=255)
                                                                  tcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                    <  1:3643450‑M1/72‑1 (MQ=255)
                                                                   cggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                   <  2:1095540‑M1/73‑1 (MQ=255)
                                                                     ggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCGGTGGGATCGGTCGGGACGGTGGGa                                 <  2:3673824‑M1/75‑1 (MQ=255)
                                                                       tccgtggtatcggtcgggtccgtgGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                               >  2:3767083‑M1/25‑101 (MQ=255)
                                                                       tccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                               >  2:585027‑M1/25‑101 (MQ=255)
                                                                       tccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                               >  1:189487‑M1/25‑101 (MQ=255)
                                                                       tccgtgggatcggtccggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCcgtggggtc                               >  2:361193‑M1/25‑98 (MQ=255)
                                                                       ggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCGGTGGGGtc                               >  2:1901930/1‑101 (MQ=255)
                                                                       ggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                               <  1:2052160/101‑1 (MQ=255)
                                                                        ccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg                              >  1:1960310‑M1/24‑101 (MQ=255)
                                                                         cgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                             <  1:3223716‑M1/79‑1 (MQ=255)
                                                                          gtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGatcggg                            <  2:2340448‑M1/80‑2 (MQ=255)
                                                                          gtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                            <  1:2875495‑M1/80‑1 (MQ=255)
                                                                           tgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                           >  2:4006304‑M1/21‑101 (MQ=255)
                                                                           tCGGGTCCGTGGGCTCGGTGGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtccgtg                           <  1:868695/101‑2 (MQ=11)
                                                                            gggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                          >  1:991847‑M1/20‑101 (MQ=255)
                                                                              gatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGatcggtcggg                        <  1:3523021‑M1/84‑2 (MQ=255)
                                                                               atccGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCGGTGGGAtcggcggggt                       >  1:1313022/4‑95 (MQ=14)
                                                                                tgggTGCGCTCCGTCGGATCGCGACGCTCCGTGGAATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATc                      <  2:1212474/98‑1 (MQ=11)
                                                                                tcggtcgtgtcggtgggAGCGGCCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAtt                      <  2:1632173/92‑2 (MQ=11)
                                                                                tGCGTGGGCTCCGTCGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                      <  2:2519909/101‑1 (MQ=255)
                                                                                 gCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                     <  2:3257619/101‑1 (MQ=255)
                                                                                 ccGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGCTCGGTCGGGTCCGTTGGAtggggcgggccc                     >  2:2746760/2‑97 (MQ=14)
                                                                                    gggggtgggtgggATCGGTGGGGTGCGTGGGATGGGGCGGGTCCGTGGCCTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                  <  1:2443565/94‑1 (MQ=11)
                                                                                     cgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                 <  1:669194‑M1/91‑1 (MQ=255)
                                                                                     cggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                 <  2:1038503‑M1/91‑1 (MQ=255)
                                                                                          ccgttgagtcggtcgggtccgtgggagcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg            <  2:1453927/93‑1 (MQ=14)
                                                                                          ccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGAtcg            >  1:3380109/3‑101 (MQ=11)
                                                                                           ctttggctcggtcgggtccgtgtgatcggtcggatccttgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcgg           <  2:3929308/94‑1 (MQ=11)
                                                                                           cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGAtccgtgggggccg           <  1:1522642‑M1/97‑6 (MQ=255)
                                                                                           cgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcgg           <  2:1750911/100‑1 (MQ=11)
                                                                                            gtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt          <  2:738258‑M1/98‑1 (MQ=255)
                                                                                             gggggccggtcggatccgtgggatcggtcggatccgtaggatCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtc         <  1:464919/95‑1 (MQ=11)
                                                                                              gggagcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg        <  2:2906956/100‑1 (MQ=14)
                                                                                                gatcggtcgggtccgtgggatcggtcgggtccgtggGGTCCGTCGGATCCGTGGGATCGGGCGGAGCCGTGGGGTCGGGCCGGTCCGTGGGAACGGGCggg      >  2:1551034/1‑101 (MQ=11)
                                                                                                gatcggtcgggtccgtgggatcggtcgggtccgtcgGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCggg      >  1:3473768/1‑101 (MQ=11)
                                                                                                 atcggtcttatccgtgggatccgtcgggtccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcggat     <  2:2271727/101‑3 (MQ=11)
                                                                                                 atcggtcgggtccgtgggatcggtcggatccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGGCGGGTCCGTGCGGTGTGTCgggg     >  2:1806838/1‑100 (MQ=11)
                                                                                                  tcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATc    <  2:2200306/101‑1 (MQ=14)
                                                                                                   cggtcggatccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc   <  1:2183644/101‑1 (MQ=11)
                                                                                                   cggccggatccgtgggatcggtcggatccgtaggatCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc   <  1:559045/101‑1 (MQ=11)
                                                                                                    ggttggatccgtgggatttgtggttttcgtgggatcTGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATTCGTGGGATCGGTCGGTTCCg  <  2:1343024/101‑1 (MQ=2)
                                                                                                    ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGGACGGTCGGGTTCg  >  2:421174/1‑101 (MQ=11)
                                                                                                    ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGCGATCGGTCGGATcag  >  2:3857848/1‑99 (MQ=14)
                                                                                                    ggtcggggccgtgggctcggtcgggtccgtgggctcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCg  <  2:2586203/101‑1 (MQ=14)
                                                                                                    |                                                                                                    
CGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCG  >  NC_014034/221861‑222061

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