Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,734,917 +A 3.8% intergenic (‑170/+23) RCAP_rcc01603 ← / ← RCAP_rcc01604 hemolysin‑III family protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,734,9171.→A 3.8% 31.6 ‑3.8 ‑0.1264intergenic (‑170/+23)RCAP_rcc01603/RCAP_rcc01604hemolysin‑III family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/1):  ref base (131/123):  total (140/124)
Fisher's exact test strand distribution p-value = 2.14e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.66e-01

Only 100 of 300 total aligned reads displayed.

GACGGGCCGGGCAGGGGGGCGCTG.CCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTCCTCGGT  >  NC_014034/1734818‑1735017
                                                                                                      |                                                                                                    
 aCGGGCCGGGCAGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.g                                                                                                     <  2:3273151/101‑1 (MQ=255)
  cgggccgggcAGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGGTATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCCTAAATATCCC.gg                                                                                                    >  1:4849401/1‑101 (MQ=14)
   gggccgggcAGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.ggg                                                                                                   >  2:3026591/1‑101 (MQ=255)
    ggccgggcAGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.gggg                                                                                                  <  1:2549031/101‑1 (MQ=255)
     tgacggacggggggCTCTG.CCCCCCGGCTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTg                                                                                                >  2:2769197/9‑101 (MQ=2)
      gacggacggggggCTCTG.CCCCCCGGCTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTga                                                                                               >  2:2594548/8‑100 (MQ=2)
      gacggacggggggCTCTG.CCCCCCGGCTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTga                                                                                               >  1:4781948/8‑100 (MQ=2)
      gacggacggggggCTCTG.CCCCCCGGCTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTga                                                                                               >  1:4662485/8‑100 (MQ=2)
      ccgggcAGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTg                                                                                                <  2:242560/101‑1 (MQ=255)
       cgggcAGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGc                                                                                               <  2:2325281/101‑1 (MQ=255)
        cggacggggggCTCTG.CCCCCCGGCTTCG..CTCCCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTgaat                                                                                             >  2:4681000/6‑98 (MQ=2)
           cAGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGGGCgggg                                                                                           >  2:1486843/1‑101 (MQ=255)
           cAGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATATTGCCATAAATATCCC.GGGGTGCgggg                                                                                           >  1:3911733/1‑101 (MQ=255)
            aGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGc                                                                                          <  1:3407255/101‑1 (MQ=255)
            aGGGGGGCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGGGCGGGGc                                                                                          >  2:1114058/1‑101 (MQ=255)
               ggggCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCaga                                                                                       >  2:3600413/1‑101 (MQ=255)
               ggggCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCaga                                                                                       <  1:411552/101‑1 (MQ=255)
                gggCGCTGCCCCCCGTCCTGCGGACTCCCCCGGGAATATTGATGGCAAGATGAACCTGC.GGCCTTCATCTTGCCATAAATACCCC.GGGGTGCGGGGCagag                                                                                      <  1:975387/101‑1 (MQ=18)
                gggCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCagag                                                                                      <  2:865399/101‑1 (MQ=255)
                 ggCGCTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGc                                                                                     <  2:4826812/101‑1 (MQ=255)
                   cgcTGCCCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGccc                                                                                   <  1:4903257/101‑1 (MQ=255)
                    gctgccccccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGcccc                                                                                  <  1:2674244/96‑1 (MQ=255)
                    gctgccccccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGcccc                                                                                  >  2:1930343/6‑101 (MQ=255)
                    gctgccccccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGcccc                                                                                  <  1:2003926/96‑1 (MQ=255)
                    gctgccccccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGcccc                                                                                  <  2:199661/96‑1 (MQ=255)
                     ctgccccccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCg                                                                                 >  1:2019543/5‑101 (MQ=255)
                     ctgccccccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCg                                                                                 >  2:3318036/5‑101 (MQ=255)
                      tgccccccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGCGCGGGGCAGAGCCCCGc                                                                                >  2:750858/4‑101 (MQ=255)
                       tccccccTTCCTGCGGACTCCCCCCTGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCt                                                                               <  2:4996964/99‑1 (MQ=255)
                       gccccccGTCCTGCGGACTCCCCCCGGGATATTTATTGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.CGGGGGCGGGGCAGAGCCCCGcc                                                                               >  1:71606/3‑100 (MQ=18)
                        ccccccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTc                                                                              <  2:4271438/100‑1 (MQ=255)
                         cccccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCa                                                                             >  2:706784/1‑101 (MQ=255)
                           cccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAAc                                                                           <  1:2315326/101‑1 (MQ=255)
                            ccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAAcc                                                                          <  2:2034933/101‑1 (MQ=255)
                            ccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAAcc                                                                          <  2:4419471/101‑1 (MQ=255)
                            ccGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAAcc                                                                          >  1:1251987/1‑101 (MQ=255)
                             cGTCCTGCGGACTCCCCCCGGGATATTGATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCg                                                                         <  2:1085547/101‑1 (MQ=255)
                              gTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGa                                                                        >  2:3492325/1‑101 (MQ=255)
                                ccTGCGGACCCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAg                                                                      <  2:1392983/101‑1 (MQ=255)
                                  tGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGAt                                                                    <  1:3461814/101‑1 (MQ=255)
                                   gCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATg                                                                   <  1:2238455/101‑1 (MQ=255)
                                     ggACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGAc                                                                 <  1:2402227/101‑1 (MQ=255)
                                     ggACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGAc                                                                 <  1:1305450/101‑1 (MQ=255)
                                      gACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACa                                                                >  2:779060/1‑101 (MQ=255)
                                      gACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACa                                                                <  2:907400/101‑1 (MQ=255)
                                      gACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACa                                                                <  1:3235838/101‑1 (MQ=255)
                                      gACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTCCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACa                                                                <  1:933511/101‑1 (MQ=255)
                                          ccccccGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGAtcggaagag                                                            >  2:3809246/1‑93 (MQ=255)
                                          ccccccGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGAt                                                            >  2:458631/1‑101 (MQ=255)
                                            ccccGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATcc                                                          <  1:4490161/101‑1 (MQ=255)
                                                 ggATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGc                                                     >  1:3702916/1‑101 (MQ=255)
                                                   atatTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGccc                                                   >  2:2834522/1‑101 (MQ=255)
                                                    tatTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCg                                                  >  1:505298/1‑101 (MQ=255)
                                                    tatTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCg                                                  >  2:4241207/1‑101 (MQ=255)
                                                    tatTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCg                                                  >  2:196058/1‑101 (MQ=255)
                                                    tatTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGCCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCg                                                  >  1:1479410/1‑101 (MQ=255)
                                                    tatTGATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCg                                                  <  2:2637414/101‑1 (MQ=255)
                                                     atTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCgg                                                 <  2:636489/101‑1 (MQ=255)
                                                      tttATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGt                                                >  1:401233/1‑101 (MQ=255)
                                                      tttATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGt                                                >  1:2515515/1‑101 (MQ=255)
                                                             cAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCt                                         >  2:4693612/1‑101 (MQ=255)
                                                              aaGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCtg                                        >  2:132855/1‑101 (MQ=255)
                                                                 aTGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCtgatg                                     <  2:313812/101‑1 (MQ=255)
                                                                 aTGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCtgatg                                     >  1:3055065/1‑101 (MQ=255)
                                                                 aTGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCtgatg                                     <  1:1930480/101‑1 (MQ=255)
                                                                 aTGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCtgatg                                     <  1:3748588/101‑1 (MQ=255)
                                                                  tGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCtgatga                                    >  2:4152610/1‑101 (MQ=255)
                                                                  tGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCtgatga                                    <  1:3191177/101‑1 (MQ=255)
                                                                    aaaCTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGc                                  >  1:1613082/1‑101 (MQ=255)
                                                                     aaCTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCa                                 <  2:3706932/101‑1 (MQ=255)
                                                                      aCTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCaa                                >  2:3371863/1‑101 (MQ=255)
                                                                       cTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAg                               <  1:496018/101‑1 (MQ=255)
                                                                       cTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAg                               <  2:2521513/101‑1 (MQ=255)
                                                                         gC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAgcg                             <  2:2180039/101‑1 (MQ=255)
                                                                            ggCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGcc                           >  1:4734712/1‑101 (MQ=255)
                                                                             gCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGccc                          <  2:387689/101‑1 (MQ=255)
                                                                               cTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGt                        <  1:3347191/101‑1 (MQ=255)
                                                                                ttCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGtt                       >  2:563751/1‑101 (MQ=255)
                                                                                 tCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCTGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTc                      <  1:771612/101‑1 (MQ=255)
                                                                                 tCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTc                      >  1:1510426/1‑101 (MQ=255)
                                                                                   aTCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCgg                    <  2:1603734/101‑1 (MQ=255)
                                                                                     cTTGCCATAAATATCCC.GGGGGGCGGGGCAGAGCCCCGCTCAACCGAAGATGCCAAGATCCGCGGCCCGGTCTTCGCTGATAAGCAAGCGCCCGTTCGGCa                  >  1:2876588/1‑101 (MQ=38)
                                                                                       tGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAgg                >  1:1650829/1‑101 (MQ=255)
                                                                                       tGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAgg                >  1:4527314/1‑101 (MQ=255)
                                                                                         ccATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTc              <  2:1816812/101‑1 (MQ=255)
                                                                                           aTAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGt            >  1:3836353/1‑101 (MQ=255)
                                                                                            tAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGtt           <  1:3600684/101‑1 (MQ=255)
                                                                                             aaaTATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGttt          >  1:1879013/1‑101 (MQ=255)
                                                                                             aaaTATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGttt          <  1:1005683/101‑1 (MQ=255)
                                                                                             aaaTATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGttt          >  1:3280625/1‑101 (MQ=255)
                                                                                             aaaTATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGttt          <  1:3318298/101‑1 (MQ=255)
                                                                                              aaTATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGtttt         <  2:1245522/101‑1 (MQ=255)
                                                                                               atatCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTc        >  2:4620596/1‑101 (MQ=255)
                                                                                               atatCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTc        >  1:2996128/1‑101 (MQ=255)
                                                                                               atatCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTc        <  2:2947900/101‑1 (MQ=255)
                                                                                                tatCCC.GGGGTGCGGGCCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTcc       <  2:1479280/101‑1 (MQ=255)
                                                                                                 atCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTCCt      <  1:1117654/101‑1 (MQ=255)
                                                                                                  tCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCAGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTCCTc     <  1:4270088/101‑1 (MQ=255)
                                                                                                     c.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCGGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTCCTCGGt  >  2:156651/1‑101 (MQ=255)
                                                                                                     c.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCGTATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTCCTCGGt  <  2:2236503/101‑1 (MQ=255)
                                                                                                      |                                                                                                    
GACGGGCCGGGCAGGGGGGCGCTG.CCCCCGTCCTGCGGACTCCCCCCGGGATATTTATGGCAAGATGAAACTGC.GGCCTTCATCTTGCCATAAATATCCC.GGGGTGCGGGGCAGAGCCCCGCTCAACCGAAGATGACAAGATCCGCGGCCCGGTTTTCGCTGATGAGCAAGCGCCCGTTCGGCAGGTCGTTTTCCTCGGT  >  NC_014034/1734818‑1735017

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