Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 243,644 T→C 1.7% intergenic (+73/+47) mogA → / ← ubiG molybdopterin biosynthesis mog protein/3‑demethylubiquinone‑9 3‑O‑methyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034243,6440T→C 1.7% 2.4 ‑0.4 ‑0.0296intergenic (+73/+47)mogA/ubiGmolybdopterin biosynthesis mog protein/3‑demethylubiquinone‑9 3‑O‑methyltransferase
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (172/119):  total (176/120)
Fisher's exact test strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

Only 100 of 343 total aligned reads displayed.

CGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCGACGACAG  >  NC_014034/243544‑243744
                                                                                                         |                                                                                                    
cgcgGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGTATATTTGAGGCAAGAT.AAAGGAGCCTGtttt                                                                                                      <  2:3187639/101‑1 (MQ=255)
 gcggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTa                                                                                                     >  2:2593805/1‑101 (MQ=255)
  cggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTAt                                                                                                    <  1:4160124/101‑1 (MQ=255)
   ggcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATc                                                                                                   <  2:4270120/101‑1 (MQ=255)
    gcgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCt                                                                                                  >  2:3956992/1‑101 (MQ=255)
     cgTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCtt                                                                                                 <  1:1226854/101‑1 (MQ=255)
      gTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTg                                                                                                >  2:3226546/1‑101 (MQ=255)
      gTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTg                                                                                                >  1:5036223/1‑101 (MQ=255)
       tttCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGc                                                                                               <  2:3653525/101‑1 (MQ=255)
        ttCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGcc                                                                                              >  2:2730504/1‑101 (MQ=255)
           gCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCa                                                                                           <  1:3295137/101‑1 (MQ=255)
            cccGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCaa                                                                                          >  2:3351887/1‑101 (MQ=255)
              cGAAGACGAAGTAAAAGGGGGCGCTTCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAAt                                                                                        <  1:4590877/101‑1 (MQ=255)
              cGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAAt                                                                                        <  1:1749010/101‑1 (MQ=255)
              cGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAAt                                                                                        <  2:3426508/101‑1 (MQ=255)
              cGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAAt                                                                                        <  2:4059990/101‑1 (MQ=255)
               gAAGACGAGGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.GTTCACCCCCGGGTTATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATg                                                                                       <  2:3023047/101‑1 (MQ=255)
                  gACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTcc                                                                                    >  1:2578809/1‑101 (MQ=255)
                    cGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCg                                                                                  <  2:3271971/101‑1 (MQ=255)
                    cGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCg                                                                                  <  1:1074510/101‑1 (MQ=255)
                     gAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCgg                                                                                 <  1:4630242/101‑1 (MQ=255)
                     gAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAACGTCCCgg                                                                                 >  1:4409625/1‑101 (MQ=255)
                      aaGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCggg                                                                                >  1:3810552/1‑101 (MQ=255)
                       aGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCgggg                                                                               <  1:306068/101‑1 (MQ=255)
                        gTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCggggg                                                                              <  2:4867084/101‑1 (MQ=255)
                         tAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGc                                                                             >  1:2432827/1‑101 (MQ=255)
                         tAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGc                                                                             >  1:4274503/1‑101 (MQ=255)
                         gggcaGGGGGCGCTGCCCCCTCTTGGCCTTTGGCCAATTCACCCCTGGGATATTTCGGGCAAGATGAAATCATGCCCCTTCATCTTGCCGCAAATATCCCg                                                                                  >  2:284437/5‑101 (MQ=0)
                         gggcaGGGGGCGCTGCCCCCTCTTGGCCTTTGGCCAATTCACCCCTGGGATATTTCGGGCAAGATGAAATCATGCCCCTTCATCTTGCCGCAAATATCCCg                                                                                  >  2:2059290/5‑101 (MQ=0)
                          aaaaGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGcc                                                                            <  2:1077508/101‑1 (MQ=255)
                           aaaGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCt                                                                           >  2:1641171/1‑101 (MQ=255)
                            aaGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTg                                                                          >  2:436296/1‑101 (MQ=255)
                               ggggCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTgggg                                                                       <  1:3387210/101‑1 (MQ=255)
                               ggggCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTgggg                                                                       >  1:1901732/1‑101 (MQ=255)
                                gggCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAAGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTggggg                                                                      >  2:5034214/1‑101 (MQ=255)
                                 ggCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGc                                                                     <  2:1400266/101‑1 (MQ=255)
                                 ggCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGc                                                                     >  1:2173616/1‑101 (MQ=255)
                                     cTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAgcg                                                                 >  1:3900464/1‑101 (MQ=255)
                                        cccccTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGccc                                                              >  2:3842010/1‑101 (MQ=255)
                                         ccccTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGcccc                                                             >  1:907412/1‑101 (MQ=255)
                                          cccTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGccccc                                                            >  2:1577741/1‑101 (MQ=255)
                                           ccTC.TGGC..TCCGCC.ATTCCCCCCGGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAATTGCCCCGGGGGCCTGGGGGCAGCGCCCCCa                                                           <  2:891065/101‑1 (MQ=38)
                                              c.tGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGcc                                                        >  1:4127317/1‑101 (MQ=255)
                                              c.tGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTAGGGGCAGCGCCCCCAGcc                                                        >  1:3473997/1‑101 (MQ=255)
                                                tGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGGAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCt                                                       >  2:4181296/1‑101 (MQ=255)
                                                cggC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCt                                                       <  2:4953482/100‑1 (MQ=255)
                                                 ggC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCtt                                                      >  1:2083339/1‑101 (MQ=255)
                                                 ggC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCta                                                      >  2:4905072/1‑100 (MQ=255)
                                                           c.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTc                                              >  2:4588565/1‑101 (MQ=255)
                                                             aTTCACCCCCGGGATATTTGAGGCAAGTT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGCGCGCCTGGGGGCAGCGCCCCCAGCCTTAGCGCTTCa                                             >  2:3594251/1‑101 (MQ=25)
                                                             aTTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCa                                             >  1:3333716/1‑101 (MQ=255)
                                                              ttCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCAc                                            >  2:2442437/1‑101 (MQ=255)
                                                              ttCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCAc                                            <  2:3748741/101‑1 (MQ=255)
                                                                   ccccGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAgg                                       <  2:3895018/101‑1 (MQ=255)
                                                                    cccGGGATATTTGAGGCNAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAggc                                      <  2:2376489/101‑1 (MQ=255)
                                                                       gggATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAggcggc                                   <  1:1141734/101‑1 (MQ=255)
                                                                          atatTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAggcggcggc                                >  2:626667/1‑101 (MQ=255)
                                                                           tatTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGgcggcggcg                               >  2:293074/1‑101 (MQ=255)
                                                                            atTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGt                              <  2:761916/101‑1 (MQ=255)
                                                                             tttGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTa                             <  1:2914116/101‑1 (MQ=255)
                                                                               tGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGt                           >  2:1681203/1‑101 (MQ=255)
                                                                               tGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGt                           >  2:2552052/1‑101 (MQ=255)
                                                                                gAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTACtt                          <  1:2891802/101‑1 (MQ=255)
                                                                                 aGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTc                         <  2:1192159/101‑1 (MQ=255)
                                                                                 aGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTc                         >  1:4950209/1‑101 (MQ=255)
                                                                                 aGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCCTAGCGTTTCACGCagg.ggcggccggagttc                        >  1:1023533/1‑93 (MQ=37)
                                                                                  ggCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCa                        >  1:2513474/1‑101 (MQ=255)
                                                                                   gaaaGAG.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCAc                       <  1:3627962/99‑1 (MQ=255)
                                                                                    cAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACg                      >  2:70505/1‑101 (MQ=255)
                                                                                    cAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACg                      >  1:365907/1‑101 (MQ=255)
                                                                                     aaGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGc                     >  2:4734193/1‑101 (MQ=255)
                                                                                       gAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTg                   <  1:3113135/101‑1 (MQ=255)
                                                                                        aT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTga                  >  2:2312638/1‑101 (MQ=255)
                                                                                        aT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTga                  >  1:3567301/1‑101 (MQ=255)
                                                                                        aT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTga                  >  1:696156/1‑101 (MQ=255)
                                                                                         t.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTgag                 >  1:1269810/1‑101 (MQ=255)
                                                                                           aaaGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTgaga                >  1:542576/1‑101 (MQ=255)
                                                                                            aaGTAGCTTGTTTTATCTGTGCGGAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGAt               <  1:2268696/101‑1 (MQ=37)
                                                                                             aGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATc              >  2:4784213/1‑101 (MQ=255)
                                                                                              ggAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATcg             <  1:4354138/101‑1 (MQ=255)
                                                                                              ggAGCCTGTTTTATCTTGCCTCAAATGCCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATcg             <  1:4758042/101‑1 (MQ=255)
                                                                                               gAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATcgc            <  1:1407282/101‑1 (MQ=255)
                                                                                               gAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTAAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATcgc            >  2:1995329/1‑101 (MQ=255)
                                                                                                aGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATcgcg           <  1:2974397/101‑1 (MQ=255)
                                                                                                  ccTTTTTTATCTTTCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATcgcgcg         <  1:1479370/101‑1 (MQ=255)
                                                                                                  ccTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATcgcgcg         <  1:4243220/101‑1 (MQ=255)
                                                                                                   cTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGcga        <  1:3842295/101‑1 (MQ=255)
                                                                                                   cTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGcga        >  1:816242/1‑101 (MQ=255)
                                                                                                   cTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGcga        >  2:4415691/1‑101 (MQ=255)
                                                                                                      ttttATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGcgacga     >  2:1280342/1‑101 (MQ=255)
                                                                                                       tttATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCgacgac    <  1:2355072/101‑1 (MQ=255)
                                                                                                       tttATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCgacgac    >  1:4932237/1‑101 (MQ=255)
                                                                                                       tttATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCgacgac    >  2:2854655/1‑101 (MQ=255)
                                                                                                       tttATCTTGCCTAAAAGGCCCCGGGGGCCTGGGGGCAGCCCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCgacgac    <  2:3567031/101‑1 (MQ=255)
                                                                                                        ttATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCGACGACa   >  1:1951366/1‑101 (MQ=255)
                                                                                                        ttATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCGACGACa   >  1:1603688/1‑101 (MQ=255)
                                                                                                        ttATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCGACGACa   >  1:4007075/1‑101 (MQ=255)
                                                                                                        ttATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCGACGACa   <  1:4051661/101‑1 (MQ=255)
                                                                                                         tATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCGACGACAg  >  1:1781475/1‑101 (MQ=255)
                                                                                                         tATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGGAGTTCACGCTGAGATCGCCCGACGACAg  >  2:2187829/1‑101 (MQ=255)
                                                                                                         |                                                                                                    
CGCGGCGTTTCGCCCGAAGACGAAGTAAAAGGGGGCGCTGCCCCCTC.TGGC..TTCGCC.ATTCACCCCCGGGATATTTGAGGCAAGAT.AAAGGAGCCTGTTTTATCTTGCCTCAAATGTCCCGGGGGCCTGGGGGCAGCGCCCCCAGCCTTAGCGTTTCACGCAGGCGGCGGCGTAGTTCACGCTGAGATCGCGCGACGACAG  >  NC_014034/243544‑243744

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG