Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,640,149 Δ1 bp 1.1% intergenic (‑125/+151) thrS ← / ← RCAP_rcc02462 threonyl‑tRNA synthetase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,640,1490C→. 1.1% 2.7 ‑0.5 ‑0.0379intergenic (‑125/+151)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (105/268):  total (107/272)
Fisher's exact test strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.36e-01

Only 100 of 442 total aligned reads displayed.

CTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGA  >  NC_014034/2640049‑2640249
                                                                                                      |                                                                                                        
 tACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.cc                                                                                                         <  2:5671268/101‑1 (MQ=255)
 tACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.cc                                                                                                         >  1:2602012/1‑101 (MQ=255)
  aCAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAAGGGGGGGCGG.cgc                                                                                                        >  1:4043569/1‑99 (MQ=255)
   cAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.cccc                                                                                                       >  1:1043422/1‑101 (MQ=255)
     aaGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.cccccc                                                                                                     <  2:1167668/101‑1 (MQ=255)
      aGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCct                                                                                                    <  2:4797741/101‑1 (MQ=255)
       gcgcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGCCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCctc                                                                                                   <  2:4463792/101‑1 (MQ=255)
         gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGGGCTG.cccccccgtt                                                                                                 >  2:2049367/1‑97 (MQ=255)
         gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.tt                                                                                                <  2:5226121/101‑1 (MQ=255)
         gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.tt                                                                                                <  1:2345415/101‑1 (MQ=255)
         gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.tt                                                                                                <  1:7491114/101‑1 (MQ=255)
         gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.tt                                                                                                <  1:5329755/101‑1 (MQ=255)
         gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAAGGGGGGGCGG.gcccccccct                                                                                                 >  2:4584069/1‑97 (MQ=255)
          cccGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTg                                                                                               <  2:8973214/101‑1 (MQ=255)
          cccGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTg                                                                                               <  2:795459/101‑1 (MQ=255)
            cGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGG.c                                                                                            <  1:5233721/101‑1 (MQ=255)
             ggTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGGGCAGGC.AAAAGGGGGGGGCTG.CCCCCCCC.TTCG.cc                                                                                           >  1:4343990/1‑100 (MQ=255)
              gTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGG.cct                                                                                          >  2:8871341/1‑99 (MQ=255)
                tGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGG.ccttc                                                                                        <  1:8883776/101‑5 (MQ=255)
                 gCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGG.ccttcg                                                                                       <  2:8647203/101‑6 (MQ=255)
                    ggTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGGGCCG.CCCCCCCC.CTGG.ccttccggc                                                                                    >  1:7561282/1‑93 (MQ=255)
                     gTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCa                                                                                    <  2:3403341/101‑1 (MQ=255)
                     gTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCa                                                                                    <  2:2322667/101‑1 (MQ=255)
                       tATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAAt                                                                                  <  1:9105103/101‑1 (MQ=255)
                       tATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAAt                                                                                  <  1:7037882/101‑1 (MQ=255)
                        atgatgCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGGGCTG.CCCCCCTC.TTGGCCTTCGGCCTAtt                                                                                 >  1:3105454/1‑101 (MQ=38)
                         tgatgCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTc                                                                                <  2:8719731/101‑1 (MQ=255)
                          gatgCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGGTG.GCCCCCTC.TTGGCCTTCGGCCAAttta                                                                               >  1:7326520/1‑99 (MQ=38)
                           atgCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCTCTG.CCCCCCTC.TTGGCCTTCGGCCAATTcac                                                                              <  2:6169001/101‑3 (MQ=255)
                             gcgcCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCAccc                                                                            <  1:8741726/101‑1 (MQ=255)
                               gcCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCAccccc                                                                          <  1:1837523/101‑1 (MQ=255)
                               gcCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGCGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCAccccc                                                                          <  1:7915626/101‑1 (MQ=37)
                                ccGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCGcccccc                                                                         <  1:9055418/101‑1 (MQ=255)
                                 cGGGTATGCGCCCCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCg                                                                        <  2:2434124/101‑1 (MQ=25)
                                   ggTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAAGGGGGGGGCGGCCCCCCCCC.TTCG.CCTCGGCCACTTACcccccccc                                                                       >  1:3424303/1‑99 (MQ=255)
                                   ggTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGGGCTG.CCCCCCTC.TCGGCCTTCGGCCAATTCCCCCCCCgcg                                                                      >  1:4507798/1‑99 (MQ=255)
                                     tATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCACCGAGGt                                                                    <  1:3771759/101‑1 (MQ=37)
                                     gaTGCGACGCCGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGt                                                                    <  1:6367536/100‑1 (MQ=37)
                                      aTGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTa                                                                   <  1:4662937/101‑1 (MQ=255)
                                      aTGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTa                                                                   <  1:6192779/101‑1 (MQ=255)
                                       tGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAt                                                                  <  1:396431/101‑1 (MQ=255)
                                        gCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAtt                                                                 <  1:4738093/101‑1 (MQ=255)
                                         cGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAttt                                                                <  2:3735096/101‑1 (MQ=255)
                                         cGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAttt                                                                <  1:5538345/101‑1 (MQ=255)
                                         cGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAttt                                                                <  2:8389229/101‑1 (MQ=255)
                                         cGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAAATCACCCCCCGAGGTAttt                                                                >  2:6302800/1‑101 (MQ=25)
                                          gACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTg                                                               <  2:5414428/101‑1 (MQ=255)
                                          gACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTg                                                               <  1:7135542/101‑1 (MQ=255)
                                            cgcgGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGa                                                             <  2:7305277/101‑1 (MQ=255)
                                             gcggcgGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATNTGGag                                                            >  1:2783970/1‑100 (MQ=255)
                                              cggcggCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGagc                                                           <  1:7801743/101‑3 (MQ=255)
                                              cggcggCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGagc                                                           <  1:3020865/101‑3 (MQ=255)
                                                gcggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCaa                                                         <  2:2789072/101‑1 (MQ=37)
                                                 cggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAg                                                        <  1:215044/101‑1 (MQ=37)
                                                 cggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGA.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAg                                                        <  2:6602884/101‑1 (MQ=25)
                                                     aGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGa                                                    <  2:529478/101‑1 (MQ=37)
                                                        cGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaacg                                                 >  2:1671333/1‑99 (MQ=25)
                                                         gtgtCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaaga                                                <  1:801475/101‑3 (MQ=25)
                                                          tgtCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaagac                                               <  2:7704924/101‑4 (MQ=25)
                                                             cAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGGGCTG.GCCCCCTCTTTGC.CTTCGGCCAAGTCACCCCCCCAGGGATTTTGAGGCAGATTAAAA.ACCAt                                           >  1:6555930/1‑101 (MQ=255)
                                                              aaCCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg                                          <  1:1553225/101‑1 (MQ=25)
                                                               aCCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg..g                                       <  2:749244/101‑1 (MQ=25)
                                                               aCCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg..g                                       <  2:8268564/101‑1 (MQ=25)
                                                               aCCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg..g                                       <  2:189009/101‑1 (MQ=25)
                                                                ccGTGGCGCTTTTGCGCAGGC.AAAACGGGGGCGCTG.CCCCCCTC.TTGGCCCTCCGCCAATTCACCCCCCCAGGCATTTGGAGCAAGATG.AAAGAGCATg..gg                                      >  2:3852988/1‑101 (MQ=21)
                                                                    gccgcTTTGCCGCCGGC.AAGGAGGGGCCGCGC.CGCCCCTC.TTGGCCTTCGCCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCa                                  <  1:4243960/99‑1 (MQ=16)
                                                                     gcgcTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCAt                                 <  2:789557/101‑1 (MQ=25)
                                                                        cTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCtt                              <  1:8549792/101‑1 (MQ=21)
                                                                          tttGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGc                            <  2:7390367/101‑1 (MQ=25)
                                                                          tttGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGc                            <  2:5784892/101‑1 (MQ=25)
                                                                           ttGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGcc                           <  2:6732890/101‑1 (MQ=25)
                                                                           ttGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGcc                           <  2:4710028/101‑1 (MQ=25)
                                                                                cAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCaaa                      <  2:3178298/101‑1 (MQ=25)
                                                                                  cccctaaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAata                    <  2:2124653/96‑1 (MQ=17)
                                                                                   gC.AAAAGGGGGGCCCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAatat                   <  1:3327943/101‑1 (MQ=21)
                                                                                   cccgaaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAatat                   >  2:2694194/5‑101 (MQ=17)
                                                                                    c.AAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATc                  <  2:8102212/101‑1 (MQ=21)
                                                                                       aaaGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATccc                <  1:6121047/101‑1 (MQ=21)
                                                                                       aaaGGGGGGCGCTG.CCACCCTC.TTGGCCTTCGGACAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATccc                <  2:6958684/101‑1 (MQ=25)
                                                                                        aaagggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAATATCCCgg              >  1:1414557/4‑101 (MQ=17)
                                                                                        aaagggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAATATCCCgg              <  1:8172126/98‑1 (MQ=17)
                                                                                          agggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAATATCCCgggg            >  1:99828/2‑101 (MQ=18)
                                                                                          agggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAATATCCCgggg            >  2:1322702/2‑101 (MQ=18)
                                                                                            ggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAATATCCCgggggg          >  1:2741492/1‑101 (MQ=18)
                                                                                            ggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATCCCggggg           >  2:4307349/1‑101 (MQ=25)
                                                                                            ggggCGCTG..CCCCCTC.TTGGCCTTTGGCCAATTC.ACCCCTGGGATATTTCGGGCAAGATGAA...ATCATGCCCCTTCATCTTGCCGCAAATATCCCGGGGGGa         <  2:3525400/101‑1 (MQ=2)
                                                                                             gggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAATATCCCgggggtg         <  1:5077480/101‑3 (MQ=18)
                                                                                              ggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAATATCCCgggggtga        >  1:6834317/1‑98 (MQ=21)
                                                                                               gcgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAATATCCCgggggtgag       <  1:4542479/101‑5 (MQ=21)
                                                                                               gcgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGG..GGTTTATCTTGCCCCAAATATCCCgggggtgag       <  1:3905923/101‑5 (MQ=21)
                                                                                                cgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGGATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGGGt       >  1:5667808/1‑101 (MQ=25)
                                                                                                 gcTG.CCCCCCTC.TTGGCCTTCGGGCCAATCAACCCCCCCGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTTCTCTTCGCTCAAACATTCCGGGGGGAggc      >  2:3021254/1‑99 (MQ=9)
                                                                                                 gcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGATGA....GGTTTATCTTGCCCCAAATATCCCgggggtgaggc     <  2:3269523/101‑7 (MQ=21)
                                                                                                 gcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGAGTc      >  2:4100169/1‑101 (MQ=25)
                                                                                                  cTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGAGTcc     >  1:1119187/1‑101 (MQ=25)
                                                                                                  cTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGAGTcc     >  2:647789/1‑101 (MQ=25)
                                                                                                  cTG.CCCCCCTC.CTGGCCTTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGAGTcc     <  1:5774876/101‑1 (MQ=25)
                                                                                                    g.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg   <  1:1261088/101‑1 (MQ=37)
                                                                                                    g.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg   <  2:550863/101‑1 (MQ=37)
                                                                                                    g.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg   <  1:3640360/101‑1 (MQ=25)
                                                                                                      |                                                                                                        
CTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGC.AAAAGGGGGGCGCTG.CCCCCCTC.TTGG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATG..GGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGA  >  NC_014034/2640049‑2640249

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