Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,018 C→G 0.4% P213P (CCG→CCC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0180C→G 0.4% 4.2 ‑2.0 0.03577P213P (CCG→CCCRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (7/16):  ref base (1675/1873):  total (1684/1893)
Fisher's exact test strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 4275 total aligned reads displayed.

GCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCG  >  NC_014034/221918‑222118
                                                                                                    |                                                                                                    
 gggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                                                                     <  1:2441284‑M1/64‑1 (MQ=255)
      tggtgggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                                                >  1:3933001‑M1/33‑101 (MQ=255)
      cggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                                <  2:5633791‑M1/69‑1 (MQ=255)
      cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                                                >  2:1175358/1‑101 (MQ=255)
      cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                                <  2:4376802/101‑1 (MQ=255)
      cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                                <  1:6860532/101‑1 (MQ=255)
        gtcggttccgtgggatcggtcgggtccgtgGGATCGGGCGGGTCCGTGGGATCGGTTGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                                                                                              >  1:8008227‑M1/31‑101 (MQ=255)
           gTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATGGGGCGGGTcgggtgg                                                                                           >  2:2309242/1‑95 (MQ=18)
           gTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                                                                                           <  1:3166605/101‑1 (MQ=255)
             gtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGgac                                                                                         <  1:5850923‑M1/76‑2 (MQ=255)
              tccgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                        >  2:1515169‑M1/25‑101 (MQ=255)
                cgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                      >  1:1627480‑M1/23‑101 (MQ=255)
                 gtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                                     <  1:1137202‑M1/80‑1 (MQ=255)
                  gcGGATGCGTGGGATGCGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                    <  2:6675672/100‑1 (MQ=11)
                   gggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGCCGGGTCCGTGGGAtcggtcg                                                                                   <  2:8816303‑M1/82‑1 (MQ=255)
                   gggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGAtcggtcg                                                                                   <  1:5640369‑M1/82‑1 (MQ=255)
                    ggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccgtcgg                                                                                  >  1:7278780‑M1/19‑101 (MQ=255)
                    ggatcggtcgggtccgtgGGATCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                  >  2:5325282‑M1/19‑101 (MQ=255)
                          gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGt                                                                            >  2:3192336/1‑101 (MQ=11)
                            cgggtccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                                          >  1:8278974/8‑101 (MQ=11)
                            cggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                                                          >  1:7657677/11‑101 (MQ=11)
                             gggtcgggcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTcgg                                                                         <  1:1000386/98‑1 (MQ=11)
                               gtcgggtggatgggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                                                       <  2:4518045/94‑1 (MQ=11)
                               cTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggt                                                                       <  1:8152441‑M2/101‑5 (MQ=255)
                                 cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                                     >  1:8603419/1‑101 (MQ=255)
                                 ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATccgtggggtcc                                                                     >  1:5527941‑M1/6‑97 (MQ=255)
                                   gtcgtATCGGGAGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                   <  2:7681288/101‑1 (MQ=11)
                                   gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGAGGGATCGGTCGCGTCCGCGGGCtcggg                                                                   >  1:1330125/1‑100 (MQ=11)
                                    tgggatcgggcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                                                  <  2:4024489/99‑1 (MQ=11)
                                    tggcgtcggccgggtccggggggtccgtcgggtccgtcggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                  <  2:2624823/96‑1 (MQ=11)
                                     gggatcggtcgggtccgtgggatcggtcgggtccgtggGATCGGTCGGATCCGTGAGATCGGTCGGGTCCGTGGGAACGGTCGGATCCGTGgggtccgtgg                                                                 <  2:3550143‑M2/100‑11 (MQ=255)
                                     gggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCGGTCGGATCCGTGGGATCGGTCGGGGCCGTGGGGGCGGGCGGCTGCGTGGGGtccgcgg                                                                 >  2:9031657/2‑96 (MQ=11)
                                     gGGATAGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCggatccgtggggtccgtgg                                                                 <  2:5995903‑M1/100‑9 (MQ=255)
                                      ggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGCGGGAgcggtggg                                                                >  2:6973487/1‑101 (MQ=11)
                                      ggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGACGGGTGCGGGGGAtggagagg                                                                >  2:8048978/1‑94 (MQ=11)
                                       gatcggtcgggtccgtgggatcggtcgggtccgtggGATCGGTCGGGTGCGTGGGATCGGGGGGGTACGGGGGATCGGTCGGTGGCGTGGGTGCGGTCggg                                                               >  1:8929578/1‑101 (MQ=11)
                                       gatcggtcggatccgtgggatcggtcgggtccgtggGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtacggcggg                                                               >  2:8485001/1‑93 (MQ=11)
                                        atcggtcgggtccgtgggatcggtcgggtccgtgggGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGACCCCGGGATCGGTCggga                                                              >  1:3008448/1‑100 (MQ=14)
                                        atcggtcgggtccgtgggatcggtcgggtccgtcggATCGGTGGGATCCGTGGGGTCGGGGGGGACCGTGGGATCGGTCGTGGCCGTCGGATCGGTGgggg                                                              >  2:5498876/1‑100 (MQ=11)
                                         tcggtcgggtccgtgggatcggtcgggtccgtcggaTCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTc                                                             <  1:6698765/101‑1 (MQ=11)
                                          cggtcgggtccgtgggatcggtcgggtccgtggggtCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                            >  1:3285669/1‑101 (MQ=11)
                                          cggtcgggtccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGGtcggtcggatag                                                            >  2:8186561/1‑97 (MQ=11)
                                           ggtgggatccgtggggtcggtcggacgcgtgggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                                           <  1:1488659/101‑1 (MQ=11)
                                           gctcgggtccgtgggatcggtcgggtccgtggggtccgTCGGATCCGTGGGATCGGGCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                                           <  1:4998485/99‑1 (MQ=14)
                                            gtcgggtccgtgggatcggtcgggtccgtgggatcgGTCGGGGCCGTGGGGTCCGTCGGGTCCGTGGGGTCGGTCGGTCCCGGGGGATGGGTTGGGGCCgg                                                          >  2:3351083/1‑100 (MQ=11)
                                              cgggtccgtgggatcggtcggctccgtgggatcggtCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcg                                                        >  2:3868436‑M2/1‑82 (MQ=255)
                                              cggatccgtgggatcggtcggatccgtaggatcggtCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtcg                                                        <  1:2124602/101‑3 (MQ=11)
                                               gggtcggggggatcggtcggatccgggggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGgg                                                       <  1:3470823/101‑1 (MQ=11)
                                               gggtccttgggagcggtcgggtgcgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtggc                                                       <  2:6533018/101‑4 (MQ=14)
                                               gggtccgtggggtcgttcgggtgggtgggattggtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCTGGTCCGTGGGATCGGTCGGGTCCGTGgg                                                       <  1:5216718/101‑1 (MQ=11)
                                               gggtccgtgggatcggtcgggtccgtgggatcggtcGGTTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGgg                                                       <  2:1330968/101‑1 (MQ=11)
                                               gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                                       <  2:866299/101‑1 (MQ=11)
                                               gggtccgtgggatcggtcgggtccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgg                                                       <  2:6203430‑M2/101‑21 (MQ=255)
                                               gggtccgtgggatcggtcggggccgtcggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGgg                                                       <  1:7489547/101‑1 (MQ=14)
                                                 gtccgtggggtcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                     <  2:7279231/101‑1 (MQ=14)
                                                 gtccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCGGTGGGGTCGGGAGGGTCGGTGgggg                                                     >  2:4196510/1‑99 (MQ=14)
                                                  tccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGCCGGGTCCGTGGGACCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatc                                                    <  1:7137609‑M2/101‑7 (MQ=255)
                                                  tccgtgggatcggtcggatccgtgggatcggtcgtGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                    >  1:2795121/1‑101 (MQ=11)
                                                   tcgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                                   >  2:2810942‑M2/2‑94 (MQ=255)
                                                     tggggatcggtcggctccgtgggatcggtcgggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggt                                                 >  2:4790625‑M2/3‑92 (MQ=255)
                                                      tgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTGGGGCCCGGCGGCtcggcg                                                >  2:4680104/1‑99 (MQ=14)
                                                      tgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcgggc                                                <  1:6410844‑M2/101‑11 (MQ=255)
                                                       ggggtcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                               <  1:2295003/101‑1 (MQ=11)
                                                       gggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAtcagatc                                               >  2:6375800/1‑96 (MQ=11)
                                                           tcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatc                                           >  2:7138644‑M2/1‑69 (MQ=255)
                                                            cggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggggcc                                          <  1:5355614‑M2/101‑17 (MQ=255)
                                                            cggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatcc                                          >  1:261796‑M2/1‑85 (MQ=255)
                                                            cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatcc                                          <  2:2475111‑M2/101‑17 (MQ=255)
                                                              gtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGggggccgtgggatcggtcgggtccgtgggatcggg                                        <  2:5176302‑M2/101‑36 (MQ=255)
                                                                tgggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccgtgg                                      <  1:4826098/100‑3 (MQ=1)
                                                                 gggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccagatc                                     >  2:2171589/1‑96 (MQ=11)
                                                                 ggatgcgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGgtcgggtccttagg                                     <  2:305285/96‑6 (MQ=11)
                                                                   gtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggat                                   <  2:1847765‑M2/101‑24 (MQ=255)
                                                                      cgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgggggatcgg                                >  1:6978756/1‑92 (MQ=11)
                                                                       gtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggggtccgtcgggtccgt                               <  1:4108365/98‑1 (MQ=14)
                                                                          ggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggggccgtgggatcggtcgg                            <  1:1772459‑M2/101‑31 (MQ=255)
                                                                            agcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAt                          <  2:2033184/99‑1 (MQ=14)
                                                                              cggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtcc                        <  1:9107317/101‑2 (MQ=11)
                                                                               ggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgcgtccgtgggatcggtcgggtccg                       <  2:7503026‑M2/101‑36 (MQ=255)
                                                                               cgtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccg                       >  2:5828835‑M2/2‑66 (MQ=255)
                                                                                 gcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtaggatcggtc                     >  2:7848315‑M2/2‑47 (MQ=255)
                                                                                  cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtgg                    >  1:363708/1‑99 (MQ=11)
                                                                                    ggtccGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtggga                  <  1:456069‑M2/99‑41 (MQ=255)
                                                                                       ccGTGGGATCGGGCGGGTCCGTGGGGTCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtaggatcggtcgggtcc               >  1:4713269‑M2/1‑41 (MQ=255)
                                                                                       ccGTGGGATCGGGCGGGTCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTcc               <  1:4268562/101‑1 (MQ=14)
                                                                                        cGTGGGATCGGTTGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcgg              <  1:6538060‑M2/101‑45 (MQ=255)
                                                                                        cGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcgg              >  2:3707080‑M2/1‑57 (MQ=255)
                                                                                          tGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc            >  2:1007399‑M2/1‑55 (MQ=255)
                                                                                          tGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcggatccgtggggtcggtg            <  1:6059748‑M2/101‑47 (MQ=255)
                                                                                          tGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtc            <  1:8752304‑M2/101‑47 (MQ=255)
                                                                                          ggggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtc            >  1:1957909‑M2/2‑55 (MQ=255)
                                                                                           gggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcg           <  1:7813016‑M2/101‑48 (MQ=255)
                                                                                           cggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg           <  1:6228562‑M2/100‑48 (MQ=255)
                                                                                            ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgggggatcggtcgg          <  2:6848005‑M2/101‑49 (MQ=255)
                                                                                              aTCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggat        >  2:1677564‑M2/1‑51 (MQ=255)
                                                                                              aTCGGTCGGGTCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggat        >  2:4936755‑M2/1‑51 (MQ=255)
                                                                                                cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcggatcc      <  2:832998‑M2/101‑53 (MQ=255)
                                                                                                cggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggttcggtcgggtccgtgggatcggtcggatcc      <  1:4424893‑M2/101‑53 (MQ=255)
                                                                                                 ggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccg     <  1:377‑M2/101‑54 (MQ=255)
                                                                                                  gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgt    <  1:8430095‑M2/101‑55 (MQ=255)
                                                                                                    |                                                                                                    
GCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCG  >  NC_014034/221918‑222118

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