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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0720C→G32.7% 291.0 ‑4.0 ‑0.1334P195P (CCG→CCCRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (51/58):  ref base (136/87):  total (188/146)
Fisher's exact test strand distribution p-value = 1.83e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.65e-01

Only 100 of 660 total aligned reads displayed.

TGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGAC  >  NC_014034/221972‑222172
                                                                                                    |                                                                                                    
tgggattggtcggatccgtgggatcggtcgtGTCCGGAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtc                                                                                                      <  1:4690306/101‑1 (MQ=11)
tgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                                                                                      >  2:2752726/1‑101 (MQ=11)
tgggatcgggcgggtccgggggagcggtcggGGCCGTCGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                                      <  1:2160132/101‑1 (MQ=14)
   gatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCGgtcggg                                                                                                   >  2:3419398/1‑101 (MQ=11)
   gatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcggg                                                                                                   <  1:2286383‑M1/57‑1 (MQ=255)
     tccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtc                                                                                                 >  1:1859743‑M1/43‑101 (MQ=255)
      cggtcgcgtccatgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtcc                                                                                                <  2:1241095/101‑1 (MQ=14)
      ccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtcc                                                                                                <  1:1553634‑M1/60‑1 (MQ=255)
       ggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGAACGGTCGGGTCCGTGGGATCGGTCGGGTCGGTCGGATCGgtcgggtccg                                                                                               >  2:3159523/1‑101 (MQ=11)
         tcggatgggtgggatcggtcggATCCGTAGGGTCGGTGGGGGCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtc                                                                                             <  2:4364003/93‑1 (MQ=11)
         tcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtc                                                                                             <  2:390780/101‑1 (MQ=11)
         tcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtc                                                                                             <  2:4053921/101‑1 (MQ=11)
         gcgggtccgtggtctcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtg                                                                                             <  2:187022/100‑2 (MQ=11)
          cgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGCGGGATCGgtcgggtccgcgg                                                                                            >  1:3265328/1‑98 (MQ=11)
          cggatccgtaggatcggtaggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTTCGTCGGAACCgtggggttcttcg                                                                                            >  2:3703934/1‑101 (MQ=2)
           gggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcgg                                                                                           <  1:1334422/101‑1 (MQ=11)
           ggatcggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                                                           <  1:451098‑M1/65‑1 (MQ=255)
           ggatccgtgggatcggtgggATCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGgtcgggtccgtagg                                                                                           <  1:506296/98‑1 (MQ=11)
            ggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcgga                                                                                          >  1:2496514/1‑100 (MQ=11)
            ggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCAGTGGGATCGGTCGGGGCCGTGGGGTCGgtcgggtccgtgggg                                                                                          >  2:1368415/1‑101 (MQ=11)
            gatcggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggg                                                                                          <  2:4597711‑M1/66‑1 (MQ=255)
              tccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggggtc                                                                                        >  1:589776/1‑101 (MQ=11)
                ggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGAtcggtcgggtccgtcgggtgcg                                                                                      <  2:1545107‑M1/70‑4 (MQ=255)
                ggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGtcgggtccgtcgggtccg                                                                                      <  1:1146189‑M1/70‑1 (MQ=255)
                 gtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgt                                                                                     >  2:3469792/1‑101 (MQ=14)
                  tggggtccgtcggATCCGTGGGATCGGTCGGATCCGTTGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtc                                                                                    <  1:4052109/101‑2 (MQ=11)
                  tgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTg                                                                                    >  1:4106899/1‑101 (MQ=14)
                   ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtcg                                                                                   >  2:2038911/1‑99 (MQ=11)
                   gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTgg                                                                                   >  2:1910374/1‑101 (MQ=11)
                   cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCTGTCGGGTCCGTGGGATCGGTCGGGTCTGTAGGATCGGTCGCGTCCGTCGGGTCCGGgg                                                                                   >  1:3930411/11‑101 (MQ=11)
                    gggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTggg                                                                                  <  2:3009424‑M1/74‑1 (MQ=255)
                    gggtccgtgggatcggtcggatccgtgGGGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTggg                                                                                  >  1:1473050‑M1/28‑101 (MQ=255)
                     ggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGa                                                                                 >  2:4139590/6‑101 (MQ=11)
                      gtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGTATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAt                                                                                <  2:4106657/97‑1 (MQ=11)
                       ttggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtc                                                                               <  2:1130226/99‑1 (MQ=11)
                       tcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtc                                                                               <  2:3522244/101‑1 (MQ=11)
                       tccgtgggatcggtggggtccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtc                                                                               >  1:363228‑M1/25‑101 (MQ=255)
                       ccgtccggATCGGGTGGATCGGCCGGATTCGTGGGGTCGGTCGGATCTGGGGGATAGGTCGGATCTGTGGGATCGGTCGGATCCGTGGGATCGGTCGGAtc                                                                               <  2:3353066/96‑1 (MQ=9)
                        ctgttgcaTCGGTCGGATCGTTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcg                                                                              <  2:4304411/94‑1 (MQ=14)
                        ccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtcc                                                                              <  1:2628884/99‑2 (MQ=11)
                         cgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccg                                                                             >  2:2267088/2‑99 (MQ=11)
                           tgggctccgtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtc                                                                           >  2:626364‑M1/21‑101 (MQ=255)
                           tcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtc                                                                           <  1:1259275/101‑1 (MQ=14)
                           tcggGTCCGTAGGATCCGTCGGGGCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtc                                                                           <  2:1719252/101‑1 (MQ=14)
                            cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtgg                                                                          >  2:683660/1‑99 (MQ=11)
                             gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                                                         >  2:569602/4‑101 (MQ=14)
                             ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtggg                                                                         >  2:1835161/1‑101 (MQ=11)
                             ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtggg                                                                         >  1:2719566/1‑101 (MQ=11)
                              ggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcggg                                                                        >  1:4356859/3‑100 (MQ=14)
                               aTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAt                                                                       >  1:299378/1‑101 (MQ=11)
                                tcggttggATCCGTCGAGGCCGTCAGATCGGTCGGATCCGTGGGATCTGTCGGGTCCGTGGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTc                                                                      <  2:759521/95‑1 (MQ=14)
                                tcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATc                                                                      <  1:1049475‑M1/86‑1 (MQ=255)
                                tcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATc                                                                      <  1:2330835‑M1/86‑1 (MQ=255)
                                  cGTGGGATCGGGCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                                    <  1:3788163/101‑1 (MQ=14)
                                  cGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTGGGGTCCg                                                                    >  2:2178612/1‑101 (MQ=14)
                                     cggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                                 <  2:2098688‑M1/91‑1 (MQ=255)
                                      ggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGGggg                                                                <  1:3353093‑M1/92‑1 (MQ=255)
                                        gtcggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                                              <  1:2204917/100‑1 (MQ=14)
                                        aTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGCGGGGTCCGGGGGATGGGTGGGGGCCGCGGGag                                                              >  2:4472268/1‑100 (MQ=11)
                                         tccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                                             <  2:365049‑M1/95‑1 (MQ=255)
                                          cggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                            <  2:4732831/101‑1 (MQ=14)
                                          cggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                            <  2:2719412/101‑1 (MQ=14)
                                           ggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTGGGGTCCGTCGGGTCCGGGGGAtggg                                                           >  2:10076/1‑101 (MQ=11)
                                            gtgggggcggTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                          <  1:198549/92‑1 (MQ=17)
                                            gtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                          <  1:4266295‑M1/98‑1 (MQ=255)
                                            gtcggGTCCGTCGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggt                                                          <  1:3704164/101‑1 (MQ=14)
                                              cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTGGGATCCGTGGGATCGGTGGGAtcggggg                                                        >  2:499940/1‑98 (MQ=14)
                                                  tCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                    <  1:3159723/101‑1 (MQ=255)
                                                  tCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTc                                                    <  2:2905214/101‑1 (MQ=17)
                                                  ctcttccgATCTGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                                                    <  1:848243/94‑1 (MQ=11)
                                                    cGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCg                                                  >  2:3823815/1‑101 (MQ=17)
                                                      tcggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCAGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTg                                                >  1:2426455/3‑101 (MQ=17)
                                                        ggcgcggtcggGTCCGTGGGGTCGGTCGTGTCGGTGGGCTCGGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTagg                                              <  2:4656988/97‑1 (MQ=11)
                                                          gtcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                                            <  2:420774/100‑3 (MQ=11)
                                                            cggtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtg                                          >  2:3322213/1‑101 (MQ=255)
                                                            cggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGGggggtcg                                          >  2:2541955/1‑98 (MQ=17)
                                                             ggtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtgg                                         >  1:3406829/1‑101 (MQ=255)
                                                              gtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCCGTCGGGTCCGTGGGGCtggt                                        >  2:2385483/1‑101 (MQ=255)
                                                              gtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGCTCGTCGGATCGGTCGGGGCCGTgggacaggc                                        >  1:3512767/1‑95 (MQ=14)
                                                                cggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggttg                                      >  1:2082038/1‑93 (MQ=14)
                                                                 gggtCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTggg                                     >  2:871123/4‑101 (MQ=255)
                                                                  ggtCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGa                                    >  1:1000596/3‑101 (MQ=255)
                                                                   gtCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAt                                   <  1:4597964/100‑1 (MQ=255)
                                                                        tGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGGtcggtc                              >  1:16866/1‑101 (MQ=255)
                                                                        tGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtc                              >  1:2801446/1‑101 (MQ=255)
                                                                          ggATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAgaggtcgg                            >  2:1047824/1‑101 (MQ=255)
                                                                              cggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCg                        <  2:1156804/101‑1 (MQ=255)
                                                                                 tcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTg                     <  2:1968942/101‑1 (MQ=255)
                                                                                  cggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGt                    >  1:3683266/1‑101 (MQ=255)
                                                                                   gggTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTc                   >  1:2938419/1‑101 (MQ=255)
                                                                                   gggTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTc                   >  1:1738001/1‑101 (MQ=255)
                                                                                    ggTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGCCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACAGCTGTcg                  >  2:1344377/1‑101 (MQ=255)
                                                                                      tCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgcc                >  1:2049397/1‑101 (MQ=255)
                                                                                          tAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgccg            <  2:3097627/101‑1 (MQ=255)
                                                                                              aTCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGCCGCCGCCGCCgg        >  2:1464031/1‑101 (MQ=255)
                                                                                                cggtcggGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGAc      >  2:2474479/1‑101 (MQ=255)
                                                                                                cggtcggGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGAc      <  1:1904082/101‑1 (MQ=255)
                                                                                                 ggtcggGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACa     <  2:2330702/101‑1 (MQ=255)
                                                                                                 ggtcggGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACa     >  2:2630131/1‑101 (MQ=255)
                                                                                                    cgggtccgtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGAc  >  2:3749083/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
TGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGAC  >  NC_014034/221972‑222172

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