BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,587,224	A→T	7.4%	Q94H (CAA→CAT)	yhhA →	DUF2756 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,587,224	0	A	T	7.4%	421.1 / 18.4	121	Q94H (CAA→CAT)	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 = 8.98e-01

Only 100 of 123 total aligned reads displayed.

CAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACT  >  NC_000913/3587078‑3587371

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

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

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

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

      aTCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgctgc                                                                                                                                              <  2:557294/148‑1 (MQ=255)

       tCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGccc                                                                                                                                            >  3:816074/1‑149 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                            aTCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTa                                                                                                                       <  3:491432/149‑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)

                                      gAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAAc                                                                                                              >  2:104841/1‑148 (MQ=255)

                                       aaaaCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCg                                                                                                            >  4:371049/1‑149 (MQ=255)

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

                                        aaaCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCg                                                                                                            >  1:471482/1‑148 (MQ=255)

                                        aaaCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGt                                                                                                           >  4:16146/1‑149 (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: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)

                                                   cAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATa                                                                                                 <  2:902908/148‑1 (MQ=255)

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

                                                     gTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGtt                                                                                              <  3:30127/149‑1 (MQ=255)

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

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

                                                        gCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCg                                                                                            <  1:560982/148‑1 (MQ=255)

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

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

                                                              gcagcaTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCATTTGCTGccgg                                                                                                                                           <  3:595468/95‑3 (MQ=255)

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

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

                                                               cagcaTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGtt                                                                                          >  3:685577/1‑143 (MQ=255)

                                                               cagcaTTTAGAAAACCAGATAAACAACAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCATTTGCTGccg                                                                                                                                            >  2:709414/1‑92 (MQ=255)

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

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

                                                                    tttAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCagcagc                                                                               <  4:291994/149‑1 (MQ=255)

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

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

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

                                                                            aCCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATAtgct                                                                        >  2:218672/1‑148 (MQ=255)

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

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

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

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

                                                                                      cAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagagg                                                              >  1:386053/1‑148 (MQ=255)

                                                                                      cAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagagg                                                              <  1:61466/148‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                        gtTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTg                                            <  1:147546/148‑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)

                                                                                                                       gggggAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCAcgc                            <  3:168679/149‑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)

                                                                                                                                  aTCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGAc                 <  3:434408/149‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                                                                           ggCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTg        >  3:232285/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)

CAGCAAATCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACT  >  NC_000913/3587078‑3587371

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