BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,587,225	A→T	7.4%	M95L (ATG→TTG)	yhhA →	DUF2756 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,587,225	0	A	T	7.4%	424.5 / 18.6	121	M95L (ATG→TTG)	yhhA	DUF2756 family protein
Reads supporting (aligned to +/- strand): ref base A (52/60); new base T (6/3); total (58/63)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 125 total aligned reads displayed.

CAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTA  >  NC_000913/3587078‑3587372

cAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCaaa                                                                                                                                                     >  1:146999/1‑148 (MQ=255)

  gCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgc                                                                                                                                                  >  4:355415/1‑149 (MQ=255)

   cAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgc                                                                                                                                                  >  2:596317/1‑148 (MQ=255)

     aaTCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgctg                                                                                                                                                <  1:596317/148‑1 (MQ=255)

     aaTCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgctg                                                                                                                                                <  2:719309/148‑1 (MQ=255)

     aaTCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgctg                                                                                                                                                >  2:681937/1‑148 (MQ=255)

        ccAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGccc                                                                                                                                             >  2:404432/1‑148 (MQ=255)

         cAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCa                                                                                                                                            >  1:678137/1‑148 (MQ=255)

         cAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCa                                                                                                                                            >  2:729520/1‑148 (MQ=255)

         cAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAAAAGCAGCATTTAGAAAACCAGATAAGCAATTATTCGCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCa                                                                                                                                            <  2:663627/148‑1 (MQ=255)

           gCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAGATCCCGCCCGGCAGCAAATGCTGCCCAac                                                                                                                                          >  1:579931/1‑148 (MQ=255)

             aaaaaGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACAcc                                                                                                                                       >  3:158205/1‑149 (MQ=255)

                aaGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgctg                                                                                                                                                >  1:345383/1‑137 (MQ=255)

                aaGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgctg                                                                                                                                                <  2:345383/137‑1 (MQ=255)

                     aTGCTGAATCAGCAACTGAAAACGCAAACGCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAAcggcg                                                                                                                                <  1:832448/148‑1 (MQ=255)

                         tGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATg                                                                                                                           >  3:32876/1‑149 (MQ=255)

                          gAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGTCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGt                                                                                                                          <  4:747739/149‑1 (MQ=255)

                          gAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATg                                                                                                                           >  2:147546/1‑148 (MQ=255)

                            aTCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGtt                                                                                                                         <  1:96692/148‑1 (MQ=255)

                            aTCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTa                                                                                                                        <  3:491432/149‑1 (MQ=255)

                             tCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTaa                                                                                                                       >  3:62418/1‑149 (MQ=255)

                              cAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTaa                                                                                                                       <  1:404432/148‑1 (MQ=255)

                              cAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTaa                                                                                                                       <  1:387510/148‑1 (MQ=255)

                              cAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATGCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGGAATCCCGCCCGGCAGCAAATGCTGCCCAACTCCAACGGCGGGATGTTaaa                                                                                                                      <  4:119178/149‑1 (MQ=255)

                                gCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAAc                                                                                                                     <  1:420914/148‑1 (MQ=255)

                                gCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACa                                                                                                                    >  3:166006/1‑149 (MQ=255)

                                  aaCTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAg                                                                                                                   >  1:576084/1‑148 (MQ=255)

                                     tGAAAACGCCAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCATCACCAACGGCGGGATGTTAAACAGCaaa                                                                                                               >  4:778499/1‑148 (MQ=255)

                                        aaaCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAt                                                                                                                                                    <  2:190679/109‑1 (MQ=255)

                                         aaCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTa                                                                                                           >  4:769146/1‑149 (MQ=255)

                                          aCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTa                                                                                                           <  2:741800/148‑1 (MQ=255)

                                            gCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAAt                                                                                                         <  1:729520/148‑1 (MQ=255)

                                               aaCTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCg                                                                                                      >  2:560982/1‑148 (MQ=255)

                                               aaCTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCg                                                                                                      >  2:507129/1‑148 (MQ=255)

                                                aCTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGa                                                                                                    <  3:167479/149‑1 (MQ=255)

                                                  tCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCCGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGCTa                                                                                                  >  4:30897/1‑149 (MQ=255)

                                                  tCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGAt                                                                                                   <  2:324377/148‑1 (MQ=255)

                                                  tCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATa                                                                                                  >  4:30127/1‑149 (MQ=255)

                                                   cAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAg                                                                                                 <  3:351440/149‑1 (MQ=255)

                                                       tGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgctgc                                                                                                                                               >  3:604558/1‑99 (MQ=255)

                                                       tGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgctgc                                                                                                                                               <  4:604558/99‑1 (MQ=255)

                                                         cAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGtt                                                                                           <  3:208132/149‑1 (MQ=255)

                                                          aaCAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGa                                                                                                    <  4:172733/139‑1 (MQ=255)

                                                          aaCAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGa                                                                                                    >  3:172733/1‑139 (MQ=255)

                                                              gcagcaTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAAt                                                                                       >  2:509352/1‑148 (MQ=255)

                                                               cagcaTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGtt                                                                                           <  4:685577/143‑1 (MQ=255)

                                                               cagcaTTTAGAAAACCAGATAAACAACAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCATTTGCTGccg                                                                                                                                             <  1:709414/93‑2 (MQ=255)

                                                                  caTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagc                                                                                   <  2:576084/148‑1 (MQ=255)

                                                                   acttAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagca                                                                                  >  1:850094/3‑148 (MQ=255)

                                                                      tAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCat                                                                              >  4:168679/1‑149 (MQ=255)

                                                                         aaaaCCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCatat                                                                            >  1:480564/1‑148 (MQ=255)

                                                                          aaaCCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATAtg                                                                           >  1:612807/1‑148 (MQ=255)

                                                                               aGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGAGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCg                                                                     >  4:620950/1‑149 (MQ=255)

                                                                                   aaaCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggaga                                                                 <  4:61009/149‑1 (MQ=255)

                                                                                   aaaCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggaga                                                                 <  4:166006/149‑1 (MQ=255)

                                                                                   aaaCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggaga                                                                 <  3:326182/149‑1 (MQ=255)

                                                                                     aCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATAtgctgctgattca                                                                <  1:509352/148‑8 (MQ=255)

                                                                                        ataatTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTa                                                                                                           >  2:16428/1‑102 (MQ=255)

                                                                                        ataatTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTa                                                                                                           <  1:16428/102‑1 (MQ=255)

                                                                                        ataatTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagaggaga                                                            <  3:620950/149‑1 (MQ=255)

                                                                                         taatTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATg                                                                                                                           >  2:799379/1‑85 (MQ=255)

                                                                                            ttCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACg                                                        <  3:726164/149‑1 (MQ=255)

                                                                                             tctcAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATAtgctgc                                                                       <  1:667570/133‑1 (MQ=255)

                                                                                             tctcAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATAtgctgc                                                                       >  2:667570/1‑133 (MQ=255)

                                                                                                        gtTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTg                                             <  1:147546/148‑1 (MQ=255)

                                                                                                         ttGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGaa                                           <  4:158205/149‑1 (MQ=255)

                                                                                                           gCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCCGCATATGCGGCCGGAGAGGAGAAACGGCGACATGCTGAATc                                         >  4:739925/1‑149 (MQ=255)

                                                                                                           gCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAAt                                          <  2:678137/148‑1 (MQ=255)

                                                                                                            cAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCa                                        <  3:371049/149‑1 (MQ=255)

                                                                                                             aGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCa                                        <  2:612807/148‑1 (MQ=255)

                                                                                                              gTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAg                                       <  2:146999/148‑1 (MQ=255)

                                                                                                              gTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGc                                      <  3:778499/149‑1 (MQ=255)

                                                                                                                       gggggAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCAcg                              <  1:218672/148‑1 (MQ=255)

                                                                                                                        ggggAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCAcgc                             >  1:408177/1‑148 (MQ=255)

                                                                                                                           gAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCa                         >  4:382132/1‑149 (MQ=255)

                                                                                                                             gCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGc                       <  4:32876/149‑1 (MQ=255)

                                                                                                                                aaaTCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCt                     >  2:41175/1‑148 (MQ=255)

                                                                                                                                    cccgcccgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagca                                                                                  >  2:248879/1‑83 (MQ=255)

                                                                                                                                    cccgcccgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagca                                                                                  <  1:248879/83‑1 (MQ=255)

                                                                                                                                      cgcccgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTc              >  3:550220/1‑149 (MQ=255)

                                                                                                                                       gttgggCAGCATTTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTc              >  1:544951/5‑148 (MQ=255)

                                                                                                                                        ttgggCAGCATTTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagag                                                                >  4:267205/4‑97 (MQ=255)

                                                                                                                                        ttgggCAGCATTTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagag                                                                <  3:267205/94‑1 (MQ=255)

                                                                                                                                          cgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGtt          <  3:567886/149‑1 (MQ=255)

                                                                                                                                          cgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGtt          <  3:650016/149‑1 (MQ=255)

                                                                                                                                          cgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGtt          <  3:884774/149‑1 (MQ=255)

                                                                                                                                          cgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGtt          <  3:768569/149‑1 (MQ=255)

                                                                                                                                          cgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGtt          <  3:878161/149‑1 (MQ=255)

                                                                                                                                           ggCAGCATTTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCggt          >  1:161142/1‑146 (MQ=255)

                                                                                                                                           ggCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTg         >  3:232285/1‑149 (MQ=255)

                                                                                                                                           ggCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTg         <  3:16146/149‑1 (MQ=255)

                                                                                                                                           ggCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTg         >  4:376312/1‑149 (MQ=255)

                                                                                                                                            gcagcaTTTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTg         >  1:914460/1‑148 (MQ=255)

                                                                                                                                              agcaAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTc                                                                                              <  4:460576/61‑1 (MQ=255)

                                                                                                                                              agcaAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTc                                                                                              >  3:460576/1‑61 (MQ=255)

                                                                                                                                               gcaAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGaaaa     >  4:384169/1‑149 (MQ=255)

                                                                                                                                               gcaAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGaaaa     <  3:72675/149‑1 (MQ=255)

                                                                                                                                                  aaTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACt   <  1:507129/148‑1 (MQ=255)

                                                                                                                                                   aTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTa  <  2:35706/148‑1 (MQ=255)

                                                                                                                                                   aTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTa  >  1:35706/1‑148 (MQ=255)

CAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTA  >  NC_000913/3587078‑3587372

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