BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,166,239	A→G	5.7%	Q271R (CAG→CGG)	murI →	glutamate racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,166,239	0	A	G	5.7%	328.5 / 6.3	105	Q271R (CAG→CGG)	murI	glutamate racemase
Reads supporting (aligned to +/- strand): ref base A (46/53); new base G (3/3); total (49/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

ACCCGGCTGGTGGATTCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGC  >  NC_000913/4166091‑4166385

aCCCGGCTGGTGGCTTCTGGCGCAGCGATTGCTCGCCGAACGGCCAGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACa                                                                                                                                                    <  6:196346/149‑1 (MQ=255)

 cccGGCTGGTGGATTCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAg                                                                                                                                                   >  6:357598/1‑149 (MQ=255)

         gtggATTCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTTTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGATTTACAGCGTTACgg                                                                                                                                           <  7:174707/149‑1 (MQ=255)

         gtggATTCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACgg                                                                                                                                           >  8:230204/1‑149 (MQ=255)

          tggATTCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGc                                                                                                                                          >  8:147021/1‑149 (MQ=255)

          tggATTCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGc                                                                                                                                          >  3:330362/1‑149 (MQ=255)

          tggATTCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGc                                                                                                                                          >  3:335969/1‑149 (MQ=255)

             aTTCTGGCGCAGCGATTGCTCGCCGAACGGCTTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGCCGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTc                                                                                                                                       >  1:145302/1‑149 (MQ=255)

             aTTCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTc                                                                                                                                       <  6:210273/149‑1 (MQ=255)

              ttCTGGCGCAGCGATTGCTCGCCGATCGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACGGCGTTACGGCTTc                                                                                                                                       >  7:57295/1‑148 (MQ=255)

              ttCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACGGCGTTACGGCTTc                                                                                                                                       >  7:64143/1‑148 (MQ=255)

              ttCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACGGCGTTACGGCTTc                                                                                                                                       >  7:57292/1‑148 (MQ=255)

                 tGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGaaa                                                                                                                                   >  4:275850/1‑149 (MQ=255)

                 tGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGaaa                                                                                                                                   >  4:236864/1‑149 (MQ=255)

                    cgcAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGc                                                                                                                                <  6:85569/149‑1 (MQ=255)

                     gcAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCt                                                                                                                               <  7:147021/149‑1 (MQ=255)

                     gcAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCt                                                                                                                               >  1:76587/1‑149 (MQ=255)

                          gATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGaaa                                                                                                                          <  5:357598/149‑1 (MQ=255)

                              gCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACt                                                                                                                      >  6:260886/1‑149 (MQ=255)

                               cTCGCCGAACGGCCTGGTTGTTAGAACATGACGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACt                                                                                                                      <  6:39643/148‑1 (MQ=255)

                               cTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTg                                                                                                                     <  2:391452/149‑1 (MQ=255)

                               cTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTg                                                                                                                     >  4:140140/1‑149 (MQ=255)

                                tCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTgg                                                                                                                    >  7:290827/1‑149 (MQ=255)

                                cccccGAACGGCCGGGGAGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTgg                                                                                                                    <  4:308985/146‑1 (MQ=255)

                                  gCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCa                                                                                                                  <  1:396189/149‑1 (MQ=255)

                                      aaCGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGttt                                                                                                              <  6:130308/149‑1 (MQ=255)

                                      aaCGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGttt                                                                                                              <  7:230204/149‑1 (MQ=255)

                                      aaCGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGttt                                                                                                              >  4:21509/1‑149 (MQ=255)

                                       aCGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCTATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGtttt                                                                                                             <  2:62515/149‑1 (MQ=255)

                                        cGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTa                                                                                                            <  1:153287/149‑1 (MQ=255)

                                        cGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTa                                                                                                            <  6:101062/149‑1 (MQ=255)

                                         ggCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTCTATGGCAATGACGCCAGGAGCTGAACAAGTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTCGCAGTTTTAg                                                                                                           >  5:303745/1‑149 (MQ=255)

                                          gCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAgg                                                                                                          >  6:149826/1‑149 (MQ=255)

                                                      gAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTAGAAAAACTGGCAGTTTTAGGCTGATTTGGtt                                                                                               >  1:243336/1‑148 (MQ=255)

                                                        aCATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGaa                                                                                            <  3:108483/149‑1 (MQ=255)

                                                          aTGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATg                                                                                          >  7:193884/1‑149 (MQ=255)

                                                             aaGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGtt                                                                                        <  3:21509/148‑1 (MQ=255)

                                                               gCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTgcg                                                                                     <  8:14771/149‑1 (MQ=255)

                                                                ccccGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTgcgc                                                                                    <  3:275850/149‑1 (MQ=255)

                                                                  ccGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCgg                                                                                  >  1:176719/1‑149 (MQ=255)

                                                                       tGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGa                                                                             <  5:260886/149‑1 (MQ=255)

                                                                          aaaaTCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGaaaa                                                                          >  6:90595/1‑149 (MQ=255)

                                                                          aaaaTCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGaaaa                                                                          >  3:63145/1‑149 (MQ=255)

                                                                                  ccGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAattat                                                                      <  4:335969/145‑1 (MQ=255)

                                                                                  ccGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAattat                                                                      <  4:330362/145‑1 (MQ=255)

                                                                                  ccGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTAtttt                                                                   <  2:38644/148‑1 (MQ=255)

                                                                                        cGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTc                                                            <  4:248915/149‑1 (MQ=255)

                                                                                         gAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGGAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTcc                                                           >  4:1887/1‑149 (MQ=255)

                                                                                         gAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTcc                                                           >  7:121328/1‑149 (MQ=255)

                                                                                         gAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTcc                                                           >  3:213905/1‑149 (MQ=255)

                                                                                           atatTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCctc                                                         <  4:332647/149‑1 (MQ=255)

                                                                                              ttGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCtt                                                       <  2:76587/148‑1 (MQ=255)

                                                                                                   ttttGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAgg                                                 >  7:394523/1‑149 (MQ=255)

                                                                                                     ttGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGc                                                <  3:140140/148‑1 (MQ=255)

                                                                                                      tGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAGCGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCg                                              >  2:140708/1‑149 (MQ=255)

                                                                                                       gTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCg                                              >  2:220715/1‑148 (MQ=255)

                                                                                                       gTATGGCAATGACGCCAGGAGCGGAACATTTATTGCCCGTGTTACAGCGTTACTGCTTCGAAACGCTCGAAAAACGGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCg                                              <  8:106417/148‑1 (MQ=38)

                                                                                                        tATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACGGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCgg                                             <  8:64143/148‑1 (MQ=255)

                                                                                                        tATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACGGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCgg                                             <  8:57295/148‑1 (MQ=255)

                                                                                                        tATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACGGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCgg                                             <  8:57292/148‑1 (MQ=255)

                                                                                                        tATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCct                                                          >  7:26917/1‑135 (MQ=255)

                                                                                                        tATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCct                                                          <  8:26917/135‑1 (MQ=255)

                                                                                                         aTGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGaa                                           >  1:14833/1‑149 (MQ=255)

                                                                                                         aTGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGaa                                           >  8:122794/1‑149 (MQ=255)

                                                                                                          tGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAAt                                          <  5:113412/149‑1 (MQ=255)

                                                                                                            gCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATaa                                        >  4:307035/1‑149 (MQ=255)

                                                                                                             cAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAAc                                       <  1:394425/149‑1 (MQ=255)

                                                                                                                tGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTcc                                    <  5:165595/149‑1 (MQ=255)

                                                                                                                tGACGCCAGGAGCCGAACAATTATTGCCCTTTTTACAGCGTTACCGCTTCGAAACGCTAGAAAAACTGGCAGTTTTAGGCTGATATGGTTGAATGTTGAGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTcc                                    <  8:156386/149‑1 (MQ=38)

                                                                                                                agACGCCAGGAGCGTAACAATTATTGGCAGTTTTCCAGCGTTACGGCTTCGAACCGCTCGAAAAACTGGGAGTTATAGTCTGATTTGGTTGAATGTTGCGCGTGCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTcc                                    <  6:92006/148‑1 (MQ=25)

                                                                                                                  aCGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGGTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCt                                  <  3:236864/149‑1 (MQ=255)

                                                                                                                      cAGGAGCTGAACAATTATTGCCAGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATaa                              >  2:291375/1‑149 (MQ=255)

                                                                                                                       aGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAAt                             >  5:396540/1‑149 (MQ=255)

                                                                                                                          aGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGCCCGGAATAACTCCCTATAATgcg                          <  8:121328/149‑1 (MQ=255)

                                                                                                                           gCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATgcgc                         <  5:149826/149‑1 (MQ=255)

                                                                                                                             tGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGcca                       >  4:176799/1‑149 (MQ=255)

                                                                                                                              gAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGcca                       >  3:230298/1‑148 (MQ=255)

                                                                                                                                  aaTTCTTGCCCGTGTTACAGCGGTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGGCAGAAAATTATTTTAAATTTCCGCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACt                  <  4:13755/149‑1 (MQ=37)

                                                                                                                                  aaTTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACt                  <  6:291586/149‑1 (MQ=255)

                                                                                                                                    ttattGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGa                >  3:230664/1‑149 (MQ=255)

                                                                                                                                       ttGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGacac             >  6:151662/1‑149 (MQ=255)

                                                                                                                                       ttGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGacac             >  5:10957/1‑149 (MQ=255)

                                                                                                                                       ttGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGGTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGaca              <  1:397885/148‑1 (MQ=255)

                                                                                                                                         gCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACg            <  4:63145/148‑1 (MQ=255)

                                                                                                                                         gCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACg            <  4:230664/148‑1 (MQ=255)

                                                                                                                                         gCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACg            <  2:176719/148‑1 (MQ=255)

                                                                                                                                         gCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTATCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACg            <  2:278743/148‑1 (MQ=255)

                                                                                                                                         gCCCGGTTTTCAGCGTTACGGCTTCGAAACGCTGGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGGTGCTCTGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCGTATAATGCGCCTCCACTGACACg            <  6:303745/148‑1 (MQ=25)

                                                                                                                                          cccGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACgg           >  1:140743/1‑148 (MQ=255)

                                                                                                                                          cccGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa          <  2:145302/149‑1 (MQ=255)

                                                                                                                                           ccGTTTTACAGCGTTACGGCTTCGAAACTCTCGAAAAACTTGCAGTTTTAGGCTGATTTGGTTGAATGTTGCCCGGTCAGAAAAATATTTTAAATATCCTCTTGGCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaa         <  1:114209/149‑1 (MQ=25)

                                                                                                                                           ccGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaa         <  1:411870/149‑1 (MQ=255)

                                                                                                                                           ccGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGAGTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTcc                                    >  4:240778/1‑122 (MQ=255)

                                                                                                                                            cGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATACTGCGCCACCACTGACACGGaa         >  5:276989/1‑148 (MQ=255)

                                                                                                                                             gTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGGTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaac        >  4:68424/1‑148 (MQ=255)

                                                                                                                                                ttACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTTACACGGaacaac     <  2:243336/148‑1 (MQ=255)

                                                                                                                                                ttACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac     <  8:394523/148‑1 (MQ=255)

                                                                                                                                                ttACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg    >  4:288389/1‑149 (MQ=255)

                                                                                                                                                  aCAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACgg   >  5:220145/1‑148 (MQ=255)

                                                                                                                                                   cAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGc  >  5:168862/1‑148 (MQ=255)

ACCCGGCTGGTGGATTCTGGCGCAGCGATTGCTCGCCGAACGGCCTGGTTGTTAGAACATGAAGCCCCGGATGCAAAATCTGCCGATGCGAATATTGCCTTTTGTATGGCAATGACGCCAGGAGCTGAACAATTATTGCCCGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGC  >  NC_000913/4166091‑4166385

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`