BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012967	3,174,256	0	T	G	24.0%	310.7 / 23.8	150	T175P (ACC→CCC)	ygjV	YgjV family protein
Reads supporting (aligned to +/- strand): ref base T (23/91); new base G (36/0); total (59/91)
Fisher's exact test for biased strand distribution p-value = 2.39e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 227 total aligned reads displayed.

AAGACGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAGGCCATTCA  >  NC_012967/3174185‑3174327

aagaCGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGt                                                                            <  2:2047755/72‑1 (MQ=255)

aagaCGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGt                                                                            <  1:6451826/72‑1 (MQ=255)

 agaCGGTCGTCTNGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGtt                                                                           <  1:1074645/72‑1 (MQ=255)

 agaCGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGtt                                                                           <  1:9856890/72‑1 (MQ=255)

 agaCGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGtt                                                                           >  1:6650675/1‑72 (MQ=255)

 aaaCGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGtt                                                                           <  2:2399860/70‑1 (MQ=255)

  gaCGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATTGCGGGA‑GGGGttt                                                                          >  1:9783044/1‑72 (MQ=255)

   aCGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGG‑GGGGgttt                                                                         >  2:8146477/1‑72 (MQ=255)

   aCGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGgttt                                                                         >  2:937562/1‑72 (MQ=255)

   aCGGTCGTCTGGCAAATAGCGACCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGtttt                                                                         <  1:5748974/72‑1 (MQ=255)

      gtcgtcTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGC‑GGGGTTTtttc                                                                      >  2:4612346/1‑72 (MQ=255)

      gtcgtcTGGCAAATAGCGCCCTCCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTtctc                                                                      >  1:9534764/1‑72 (MQ=255)

      gtcgtcTCGCAGCCAGCGCCCGGCATAATTAATTGCTTAGCCCCTTC‑T‑GCGGATGGCGGGA‑GGGGTTTtctc                                                                      <  2:4051083/72‑1 (MQ=255)

        cgtcTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGG‑GGGGTTTTTTCT‑a                                                                   >  1:77952/1‑72 (MQ=255)

        cgtcTGGAAAATAGCGCCCTGCGTAATTAATTGCTTACCCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑a                                                                   <  2:6727489/72‑1 (MQ=255)

         gtcTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑TGGGGTTTCTCT‑Ac                                                                  >  1:4859268/1‑72 (MQ=255)

         gtcTGGAAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTctcttc                                                                   >  2:7865860/1‑70 (MQ=255)

          tcTGACAAATAGCGCCGTGCTTAATTAATTGCTTACCCCTTTC‑C‑GTCTAGGGGGGGG‑GGGGTTTTCTCT‑AC‑t                                                                <  1:4409132/72‑1 (MQ=255)

            tGGCAAATCGCGCCCTGCGTAATTATTTGCTTACCCCCTTT‑C‑GTCTATGGGGGGA‑GGGGTTTTCTCT‑AC‑ttt                                                              <  1:6932387/72‑1 (MQ=255)

            tGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGATTTCTCT‑At‑ttt                                                              >  2:7516830/1‑72 (MQ=255)

            tGGAAAATGGCGCCCTGCTTAATTATTTCCTTCCCCCCTTT‑C‑GTCTGGGGGGGGG‑GGGGTTTTCTCT‑AC‑ttt                                                              <  1:9059580/72‑1 (MQ=255)

             ggCAAATAGCGCCCTGCGTAATTAATTGCTTACCCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTa                                                             <  1:6885327/72‑1 (MQ=255)

             ggCAAATAGCGCCCTGCGTAATTAATTGCTTACCCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTa                                                             <  1:4259066/72‑1 (MQ=255)

                aaaTAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGNGGTTTCTCT‑AC‑CTTTAAc                                                          >  1:3291263/1‑72 (MQ=255)

                aaaTAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTCTTT‑TC‑TTTTAAc                                                          >  2:8799921/1‑72 (MQ=255)

                 aaTAGNGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGAGGGGGTTTTTTTT‑CC‑TTTaaag                                                          >  1:942287/1‑71 (MQ=255)

                 aaTAGCGCCTTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTAGGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACg                                                         <  2:5095536/72‑1 (MQ=255)

                 aaTAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACg                                                         >  2:9493651/1‑72 (MQ=255)

                   tAGCGCCTTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGa                                                       <  2:7378118/72‑1 (MQ=255)

                   gaGCGCCCTGCGTAATGAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGa                                                       <  1:2494086/71‑1 (MQ=255)

                    ggcgcCCGTCGTAATTAATTTC‑TAACCCCTTTGC‑TTGTAGGGGGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAt                                                      <  2:2247994/71‑1 (MQ=255)

                     gcgcCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑c                                                    <  1:9623951/72‑1 (MQ=255)

                     gcccccTGCGTATTTAATTGCTTACCCCCTTT‑C‑GTCTGGGGGGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑c                                                    <  2:4389564/69‑1 (MQ=255)

                      cgcACTGCATATTTATTTGCATACCCCCTTT‑C‑GTCTAGGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑Ca                                                   <  1:2360946/72‑1 (MQ=255)

                       gcCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTCT‑TC‑TTTAAACGGAA‑Caa                                                  >  1:9007222/1‑72 (MQ=255)

                        cccTGCGTAATTAATTGCTTACCCCCTTT‑C‑GTCTATGGGGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAAt                                                 <  1:2099594/72‑1 (MQ=255)

                         ccTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTCTAAC‑TTTAAA‑GGAT‑CAAtt                                                >  2:6216967/1‑72 (MQ=255)

                          cTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTTTCT‑NC‑TTTAAAAGGGT‑AAATTc                                               >  1:3802271/1‑72 (MQ=255)

                          cTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTCT‑AC‑TTTAAACGGAT‑CAATTc                                               >  2:9086355/1‑72 (MQ=255)

                          cTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTCTCT‑AC‑TTTAAACGGAT‑CAATTc                                               >  1:7635655/1‑72 (MQ=255)

                          cTGCGTAA‑TGATTGCCTACCCCCCGT‑CGGTCAAGGGCGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTc                                               <  2:5461206/72‑1 (MQ=255)

                           tGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTAGGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTcc                                              <  1:9169252/72‑1 (MQ=255)

                            gCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTCT‑TC‑TTTTAACGGAT‑CAATTtcc                                             >  2:8694683/1‑72 (MQ=255)

                            gCGTAATTAATTCCTTACCCCCTTT‑C‑GTCTAGGGGGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTccc                                             <  1:3211775/72‑1 (MQ=255)

                             cGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTCT‑TC‑TTTAAAAGGGT‑CAAATcccc                                            >  2:6517620/1‑72 (MQ=255)

                              gTAATTCATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGG‑GGGGGTTTTTCT‑AC‑TTTTAAAGGAT‑CAATTccccc                                           >  1:1732728/1‑71 (MQ=255)

                               taattaatTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTTTCT‑TC‑TTTTAACGGAT‑CAATTtcccct                                          >  1:3468517/1‑70 (MQ=255)

                                aattaattGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTCT‑TC‑TTTAAAAGGATACATTTccccct                                          >  2:723692/1‑70 (MQ=255)

                                 attaattGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GTGGTTTTTTCT‑AC‑TTTTAAAGGAT‑CAAATCCCCtttt                                        >  1:7000805/1‑72 (MQ=255)

                                 attaattGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCtttt                                        <  1:7474986/72‑1 (MQ=255)

                                 attaattGCTTAACCCCTTT‑C‑GTCTATGCCGGGA‑GGGGGTTTCTCT‑AC‑TTTAAAAGGAT‑CAAATTCCCtttt                                        >  2:5912477/1‑72 (MQ=255)

                                   taattTCTTAACCCCTTT‑C‑GTCTATGGGGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCt                                      <  1:4314749/72‑1 (MQ=255)

                                   taattGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCCTTTCt                                      >  2:3409419/1‑72 (MQ=255)

                                    aattGCTTCCCCCCTTT‑C‑GTCTGTGGGGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTg                                     <  1:1757929/72‑1 (MQ=255)

                                    aattGCTTACCCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTg                                     <  2:7078782/72‑1 (MQ=255)

                                     gttGCTTACCCCCTTT‑C‑GTCTAGGGGGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGc                                    <  1:3019344/71‑1 (MQ=255)

                                     attGCTTACCCCCTTT‑C‑GTCGAGGGGGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGc                                    <  1:4792127/72‑1 (MQ=255)

                                     attGCTTAACCCCTTT‑C‑GTCTGGGGGGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGc                                    <  2:5048085/72‑1 (MQ=255)

                                      ttGCTCACCCCTTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCa                                   <  2:8043052/72‑1 (MQ=255)

                                         cTTACCCCCTTT‑G‑TTCTGGGGGGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATcc                                <  2:3966175/72‑1 (MQ=255)

                                         cTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTTT‑AC‑TTTTAAAGGAT‑CAATTTCCCCTTTCTGCAttc                                >  1:1282477/1‑70 (MQ=255)

                                         cTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTCTCT‑TC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCCTcc                                >  2:7675082/1‑72 (MQ=255)

                                           tAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTTTCT‑AC‑TTTAAACGGCA‑CAATTCCCCTTTTCTGCATCCGc                              >  1:7678847/1‑72 (MQ=255)

                                           tAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGc                              <  1:8100729/72‑1 (MQ=255)

                                            ccccccTTG‑C‑GTCTATGGGGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGcc                             <  1:2668209/70‑1 (MQ=255)

                                            acccccTTT‑CGGTGTAGGGGGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGc                              <  1:3671031/70‑1 (MQ=255)

                                            acccccTTT‑C‑GTTTAGGGGGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGcc                             <  2:6920019/70‑1 (MQ=255)

                                            aaCCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGcc                             >  1:5423485/1‑72 (MQ=255)

                                             aCCCCTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTCT‑AC‑TTTTAAAGGAT‑CAATTCCCCCTTTCTGCATACCCCa                            >  1:724723/1‑72 (MQ=255)

                                              ccccTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTCT‑TCTTTTAAAGGAAT‑AAATCCCCCTTTTTTTCTTCCCCCa                            >  2:2904315/1‑72 (MQ=255)

                                              ccccTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTTTCT‑AA‑TTTAAACGGAT‑CAAATTCCCCTTTCTTCATCCGcccg                           >  2:3899083/1‑70 (MQ=255)

                                              ccccTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAg                           >  1:8875951/1‑72 (MQ=255)

                                               cccTTT‑C‑GTCTATGGCGGGG‑GGGGGTTTTTCT‑AC‑TTTTAAAGGGT‑CAAATCCCCTTTTCTTCATCCGCCCGa                          >  1:8607787/1‑72 (MQ=255)

                                               cccTTT‑C‑GTCTATGGCGGGA‑GGGGGTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTTCATCCGCCAGa                          >  2:4902555/1‑72 (MQ=255)

                                               cccTTC‑C‑GTCAGGGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGa                          <  1:409406/72‑1 (MQ=255)

                                                ccTTT‑C‑GTCTATGGGGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCAAGaa                         <  1:7039979/72‑1 (MQ=255)

                                                  ttt‑C‑GTCTATGGCGGGA‑GGGATTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAAc                       <  1:2646793/72‑1 (MQ=255)

                                                        gTCTATGGCGGGG‑GGGGGTTTTTCT‑TC‑TTTAAACGGAT‑CAAATCCCCTTTTCTTCATCCGCCAGAAACGgac                   >  2:8395864/1‑71 (MQ=255)

                                                        gTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGgat                   <  1:5148158/72‑1 (MQ=255)

                                                        gTCTATGGCGGGA‑GGGGGTTTTTCT‑AC‑TTTAAAAGGAT‑CAATTCCCCCTTTCTGCATTCCCCACAAACgggn                   >  2:3666264/1‑70 (MQ=255)

                                                          cTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGACAGAAACGgatga                 <  1:1914703/72‑1 (MQ=255)

                                                          cTATGGCGGGA‑GGGGGTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCCTTTCTTCATCCCCCCGAAACGggtga                 >  1:2602626/1‑72 (MQ=255)

                                                              ggCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAAAGGAT‑CAATTCCCCCTTTCTGCATCCGCCAGAAAACGCTGATAtt             >  2:213361/1‑72 (MQ=255)

                                                               gggggg‑ggGTTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTc            <  1:5532562/70‑1 (MQ=255)

                                                               gCGGGG‑GGGGTTTTCTCT‑AC‑TTTAAACGGAG‑CAATTCCCCTTTTCTGCAGCCGCCAGAAACGGATGATATTc            <  2:3801889/72‑1 (MQ=255)

                                                               gCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACTGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTc            <  2:6879077/72‑1 (MQ=255)

                                                                ggggT‑GGTGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCa           <  1:6748435/71‑1 (MQ=255)

                                                                cGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCa           >  1:8297306/1‑72 (MQ=255)

                                                                 gggA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAg          >  2:5677904/1‑72 (MQ=255)

                                                                 gggA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAg          >  1:9298908/1‑72 (MQ=255)

                                                                  ggA‑GGGTTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAgg         <  1:5686500/72‑1 (MQ=255)

                                                                  ggA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAgg         >  2:9837144/1‑72 (MQ=255)

                                                                    a‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAGGcc       >  1:339880/1‑72 (MQ=255)

                                                                       gggTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAGGCCAt     >  2:6239582/1‑72 (MQ=255)

                                                                        gtttttCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAGGCCAtt    <  1:3144526/70‑1 (MQ=255)

                                                                        ggTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAGGCCAtt    <  1:8143109/72‑1 (MQ=255)

                                                                         gTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAGGCCATTc   <  1:7401827/72‑1 (MQ=255)

                                                                         gTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAGGCCATTc   >  2:3213063/1‑72 (MQ=255)

                                                                          ttttCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAGGCCATTCa  <  1:6015663/72‑1 (MQ=255)

                                                                          ttttCTCT‑AC‑TTAAACCGGAC‑CAATCCCCTTTTTCGGCACCCGCCAGAAACGGATGATATTCAGGCCATTCa  <  2:1259076/72‑1 (MQ=255)

AAGACGGTCGTCTGGCAAATAGCGCCCTGCGTAATTAATTGCTTAACCCCTTT‑C‑GTCTATGGCGGGA‑GGGGTTTTCTCT‑AC‑TTTAAACGGAT‑CAATTCCCCTTTTCTGCATCCGCCAGAAACGGATGATATTCAGGCCATTCA  >  NC_012967/3174185‑3174327

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 17 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`