BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,573,189	T→C	5.0%	pseudogene (776/1503 nt)	yjiT →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,189	0	T	C	5.0%	466.6 / 10.1	141	pseudogene (776/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand): ref base T (66/67); new base C (3/4); total (70/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

Only 100 of 143 total aligned reads displayed.

GTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCTTT  >  NC_000913/4573044‑4573337

gTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                                   <  4:217165/149‑1 (MQ=255)

 tacttactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTagag                                                                                                                                                  >  2:409675/1‑149 (MQ=255)

 tacttactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTagag                                                                                                                                                  <  1:761582/149‑1 (MQ=255)

    ttagtCATAACCTAATAAGTCACAATGAACCAGTACAACTGCTTGAGAAGTGTCACCAAACCTGGCGTAGCGAGTTTCCAATACCGTGAGCTGATGAAACAGGAACGCATTTTTTGAGTTGACTGCTATGCGCAGCCTCTGAATAGGCt                                                                                                                                               >  2:192618/1‑149 (MQ=255)

      actCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                                                                             <  1:632683/149‑1 (MQ=255)

      actCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgc                                                                                                                                             <  3:299622/149‑3 (MQ=255)

        tCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGc                                                                                                                                           >  2:636281/1‑149 (MQ=255)

          acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTCGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCt                                                                                                                                         >  4:737635/1‑149 (MQ=255)

          acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCt                                                                                                                                         >  1:354703/1‑149 (MQ=255)

          acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCt                                                                                                                                         >  4:560717/1‑149 (MQ=255)

             aCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGa                                                                                                                                      <  2:511340/149‑1 (MQ=255)

                taataaATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTa                                                                                                                                   <  3:693112/149‑1 (MQ=255)

                     aaTCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                                   <  3:625972/128‑1 (MQ=255)

                     aaTCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                                   >  4:625972/1‑128 (MQ=255)

                     aaTCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaa                                                                                                                              >  1:352848/1‑149 (MQ=255)

                        cacaAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaaaaa                                                                                                                           >  3:655105/1‑149 (MQ=255)

                          caAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATa                                                                                                                         >  1:381171/1‑149 (MQ=255)

                          caAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATa                                                                                                                         >  3:492765/1‑149 (MQ=255)

                          caAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATa                                                                                                                         >  1:89881/1‑149 (MQ=255)

                           aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                        >  4:516895/1‑149 (MQ=255)

                           aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                        <  2:198717/149‑1 (MQ=255)

                           aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                        >  2:592552/1‑149 (MQ=255)

                           aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                        >  2:402885/1‑149 (MQ=255)

                                aCCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTa                                                                                                                   >  2:326793/1‑149 (MQ=255)

                                 ccAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAc                                                                                                                  <  2:381171/149‑1 (MQ=255)

                                  cAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcg                                                                                                                 >  3:506682/1‑149 (MQ=255)

                                    gTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcgcg                                                                                                               <  1:84539/149‑1 (MQ=255)

                                    gTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcgcg                                                                                                               >  2:413987/1‑149 (MQ=255)

                                       caacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg                                                                                                            <  2:461164/149‑1 (MQ=255)

                                        aacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTAc                                                                                                                                  >  2:710686/1‑126 (MQ=255)

                                        aacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                                                                           >  1:749926/1‑149 (MQ=255)

                                        aacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGGTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                                                                           >  4:751798/1‑149 (MQ=255)

                                         acaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                                                                             <  3:569595/114‑1 (MQ=255)

                                         acaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                                                                             >  4:569595/1‑114 (MQ=255)

                                         acaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGcc                                                                                                          >  1:345120/1‑149 (MQ=255)

                                           aGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                                                                        <  1:531226/149‑1 (MQ=255)

                                           aGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                                                                        <  2:105789/149‑1 (MQ=255)

                                             cTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt                                                                                                      <  3:60422/149‑1 (MQ=255)

                                             cTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt                                                                                                      <  3:80611/149‑1 (MQ=255)

                                                 agaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAt                                                                                                  <  3:229234/149‑1 (MQ=255)

                                                 agaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAt                                                                                                  <  3:78861/149‑1 (MQ=255)

                                                  gaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt                                                                                                 <  3:577630/149‑1 (MQ=255)

                                                  gaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt                                                                                                 <  1:453784/149‑1 (MQ=255)

                                                     ggTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGt                                                                                              <  4:643782/149‑1 (MQ=255)

                                                      gTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTc                                                                                             <  2:125229/149‑1 (MQ=255)

                                                            ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                                                                       <  1:402885/149‑1 (MQ=255)

                                                            ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                                                                       >  2:311605/1‑149 (MQ=255)

                                                              tACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGc                                                                                     <  2:89881/149‑1 (MQ=255)

                                                                 cTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATc                                                                                  <  1:602432/149‑1 (MQ=255)

                                                                   ggCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCa                                                                                <  2:330452/149‑1 (MQ=255)

                                                                      gAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGAt                                                                             >  1:443985/1‑149 (MQ=255)

                                                                      gAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGAt                                                                             <  2:227264/149‑1 (MQ=255)

                                                                        aGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGa                                                                           >  4:155836/1‑149 (MQ=255)

                                                                            aGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGACAAGCATCCAGATGAGCtt                                                                       >  4:387175/1‑149 (MQ=255)

                                                                             gTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTa                                                                      >  1:164289/1‑149 (MQ=255)

                                                                              tttCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAg                                                                     >  2:351476/1‑149 (MQ=255)

                                                                                tCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAggg                                                                   <  2:212334/149‑1 (MQ=255)

                                                                                 ccAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGt                                                                  <  1:425832/149‑1 (MQ=255)

                                                                                  cAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTa                                                                      <  2:353605/144‑1 (MQ=255)

                                                                                  cAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTg                                                                 >  1:358998/1‑149 (MQ=255)

                                                                                   aaTACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAGATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGa                                                                >  1:552977/1‑149 (MQ=255)

                                                                                   aaTACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGa                                                                >  1:641420/1‑149 (MQ=255)

                                                                                    aTACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGAt                                                               <  2:763453/149‑1 (MQ=255)

                                                                                        cGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCAt                                                                                   >  1:728031/1‑125 (MQ=255)

                                                                                        cGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCAt                                                                                   <  2:728031/125‑1 (MQ=255)

                                                                                         gTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTAt                                                          >  3:341336/1‑149 (MQ=255)

                                                                                         gTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTAt                                                          <  3:547964/149‑1 (MQ=255)

                                                                                                gatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTcc                                                   <  2:228900/149‑1 (MQ=255)

                                                                                                gatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTcc                                                   >  2:706873/1‑149 (MQ=255)

                                                                                                   gaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTGGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGa                                                <  3:476571/149‑1 (MQ=255)

                                                                                                       cAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGaa                                            <  1:312881/149‑1 (MQ=255)

                                                                                                        aGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAgg                                                                    <  4:155841/124‑1 (MQ=255)

                                                                                                         ggAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAgg                                                                    >  4:725214/1‑123 (MQ=255)

                                                                                                         ggAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAgg                                                                    <  3:725214/123‑1 (MQ=255)

                                                                                                         ggAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGt                                          <  2:357493/149‑1 (MQ=255)

                                                                                                         ggAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGt                                          <  1:401067/149‑1 (MQ=255)

                                                                                                            aCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttc                                       <  3:300876/149‑1 (MQ=255)

                                                                                                               cATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcctt                                    >  3:18417/1‑149 (MQ=255)

                                                                                                                aTTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCAt                                                                                   >  3:729053/1‑101 (MQ=255)

                                                                                                                aTTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCAt                                                                                   <  4:729053/101‑1 (MQ=255)

                                                                                                                aTTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCt                                                                        <  4:270531/112‑1 (MQ=255)

                                                                                                                aTTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCt                                                                        >  3:270531/1‑112 (MQ=255)

                                                                                                                 ttttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcctttc                                  <  4:343600/149‑1 (MQ=255)

                                                                                                                   ttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTttcctttcct                                >  1:488153/1‑149 (MQ=255)

                                                                                                                   ttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTttcctttcct                                <  4:181226/149‑1 (MQ=255)

                                                                                                                        aaTGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAAc                           >  4:207213/1‑149 (MQ=255)

                                                                                                                          tGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACaa                         >  4:499274/1‑149 (MQ=255)

                                                                                                                            gACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGc                       >  1:579351/1‑149 (MQ=255)

                                                                                                                                      cgcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCaa             <  2:488153/149‑1 (MQ=255)

                                                                                                                                       ttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAgg                                                                    <  4:385887/91‑1 (MQ=255)

                                                                                                                                       ttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAgg                                                                    >  3:385887/3‑93 (MQ=255)

                                                                                                                                       ttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAAc            >  3:596736/3‑149 (MQ=255)

                                                                                                                                       gcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGTGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAAc            >  3:464122/1‑149 (MQ=255)

                                                                                                                                        taGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACt           >  4:348301/2‑149 (MQ=255)

                                                                                                                                          gCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGa                                                                           >  4:302479/1‑83 (MQ=255)

                                                                                                                                             tctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTAcgcg      >  1:601214/1‑149 (MQ=255)

                                                                                                                                              ctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACACATTGCCAATTCTATTGGTCAGAAAAGCATCCAGGTGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTAAACTACgcgc     <  4:492765/149‑1 (MQ=255)

                                                                                                                                                gTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTAtct                                                        <  2:300285/96‑1 (MQ=255)

                                                                                                                                                gTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTAtct                                                        >  1:300285/1‑96 (MQ=255)

                                                                                                                                                 tAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCttt  <  4:341336/149‑1 (MQ=255)

GTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCTTT  >  NC_000913/4573044‑4573337

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