BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,573,189	T→C	6.9%	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	6.9%	319.4 / 11.1	102	pseudogene (776/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand): ref base T (47/48); new base C (4/3); total (51/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

Only 100 of 109 total aligned reads displayed.

TTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCTT  >  NC_000913/4573041‑4573336

ttAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                                                                     <  2:701358/149‑1 (MQ=255)

ttAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                                                                     >  1:138936/1‑149 (MQ=255)

ttAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                                                                     <  2:18309/149‑1 (MQ=255)

ttAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                                                                     >  3:559627/1‑149 (MQ=255)

ttAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                                                                     <  3:24484/149‑1 (MQ=255)

   gTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                                  <  1:687701/149‑1 (MQ=255)

   gTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                                  <  1:242936/149‑1 (MQ=255)

   gTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                                  <  1:534878/149‑1 (MQ=255)

   gTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                                  >  3:42093/1‑149 (MQ=255)

   gTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTAGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                                  >  1:453962/1‑149 (MQ=255)

    tacttactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTagag                                                                                                                                                 <  2:272932/149‑1 (MQ=255)

       ttactCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTATATGATGAAACAGGTACGCATTTTTTTTAATGGAGTGCTATGGGCAGTCCCTGTAGAGGc                                                                                                                                               >  2:702411/1‑149 (MQ=255)

         actCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGATTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                                                                            <  2:567960/149‑1 (MQ=255)

          ctCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAg                                                                                                                                           <  2:32134/149‑1 (MQ=255)

           tCACAACTTAATAAATCACAAAGAACCAGTACAAAAGCGTGAGAAGGTTAACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATTATGAAACACGAACGCA‑TTTTTTGAAGGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGc                                                                                                                                          >  2:493306/1‑149 (MQ=255)

           tCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGc                                                                                                                                          >  3:239766/1‑149 (MQ=255)

           tCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGc                                                                                                                                          >  2:544776/1‑149 (MQ=255)

           tCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgc                                                                                                                                          >  1:70924/1‑145 (MQ=255)

            cacaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgca                                                                                                                                         <  1:46526/149‑6 (MQ=255)

            cacaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTACAGAGGCtgcgca                                                                                                                                         <  3:192264/149‑6 (MQ=255)

             acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCt                                                                                                                                        >  1:192543/1‑149 (MQ=255)

             acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGCTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCt                                                                                                                                        >  4:336174/1‑149 (MQ=255)

               aaCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                                                                      >  1:367545/1‑149 (MQ=255)

               aaCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                                                                      >  2:34407/1‑149 (MQ=255)

               aaCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCg                                                                                                                                      >  1:568070/1‑149 (MQ=255)

                aCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGTCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTGGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGa                                                                                                                                     >  4:572719/1‑149 (MQ=255)

                 ccTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGAc                                                                                                                                    >  3:662587/1‑149 (MQ=255)

                   taataaATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTa                                                                                                                                  >  4:324049/1‑149 (MQ=255)

                   taataaATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTa                                                                                                                                  >  2:120213/1‑149 (MQ=255)

                   taataaATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTa                                                                                                                                  >  3:180674/1‑149 (MQ=255)

                      taaATCACAAAGAACCAGTACAACAGCTTGAGAACGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGACGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaa                                                                                                                               >  2:479838/1‑149 (MQ=255)

                              aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                       >  2:186035/1‑149 (MQ=255)

                              aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                       >  3:81556/1‑149 (MQ=255)

                              aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                       <  4:662587/149‑1 (MQ=255)

                                  aaCCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGt                                                                                                                   >  4:208065/1‑149 (MQ=255)

                                    ccAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAc                                                                                                                 <  4:532108/149‑1 (MQ=255)

                                    ccAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAc                                                                                                                 <  1:328196/149‑1 (MQ=255)

                                           aacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAaga                                                                                                                     >  2:119860/1‑138 (MQ=255)

                                           aacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAaga                                                                                                                     <  1:119860/138‑1 (MQ=255)

                                              aGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                                                                       >  3:529406/1‑149 (MQ=255)

                                              aGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCaa                                                                                                       >  1:242120/1‑149 (MQ=255)

                                                cTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt                                                                                                     >  1:426680/1‑149 (MQ=255)

                                                cTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt                                                                                                     <  2:419487/149‑1 (MQ=255)

                                                  tGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCt                                                                                                   <  1:105333/149‑1 (MQ=255)

                                                     gaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt                                                                                                <  2:511546/149‑1 (MQ=255)

                                                     gaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAtt                                                                                                <  4:126010/149‑1 (MQ=255)

                                                         gTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCATCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTc                                                                                            <  1:693643/149‑1 (MQ=255)

                                                          ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCa                                                                                           <  3:351286/149‑1 (MQ=255)

                                                          ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCa                                                                                           <  4:191886/149‑1 (MQ=255)

                                                              cccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaa                                                                                       >  2:444242/1‑149 (MQ=255)

                                                               ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTa                                                                                                                  <  3:600164/121‑1 (MQ=255)

                                                               ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTa                                                                                                                  >  4:600164/1‑121 (MQ=255)

                                                               ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                                                                      >  2:396399/1‑149 (MQ=255)

                                                               ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                                                                      <  2:453962/149‑1 (MQ=255)

                                                               ccTACCTGGCGAAGCGAGATTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                                                                      <  2:296764/149‑1 (MQ=255)

                                                                 tACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACCGGTACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGATAAGc                                                                                    >  4:371495/1‑149 (MQ=255)

                                                                    cTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATc                                                                                 <  4:488596/149‑1 (MQ=255)

                                                                     tGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATcc                                                                                <  3:18455/149‑1 (MQ=255)

                                                                      ggCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCa                                                                               >  1:19195/1‑149 (MQ=255)

                                                                        cGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGa                                                                             <  2:169089/149‑1 (MQ=255)

                                                                         gAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGAt                                                                            >  1:106911/1‑149 (MQ=255)

                                                                           aGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGa                                                                          <  2:192543/149‑1 (MQ=255)

                                                                           aGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGa                                                                          <  3:600950/149‑1 (MQ=255)

                                                                               aGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                       >  4:604830/1‑100 (MQ=255)

                                                                               aGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                       <  3:604830/100‑1 (MQ=255)

                                                                                    ccAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGt                                                                 <  4:42093/149‑1 (MQ=255)

                                                                                      aaTACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGa                                                               >  3:106582/1‑149 (MQ=255)

                                                                                      aaTACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGa                                                               >  2:542737/1‑149 (MQ=255)

                                                                                        tACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCGAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAATAGCATCCAGATGAGCTTAGGGTGATa                                                             <  4:582792/149‑1 (MQ=255)

                                                                                             ttAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcgc                                                                                                               >  4:456064/1‑94 (MQ=255)

                                                                                                       aaaCAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGAc                                              <  3:41644/149‑1 (MQ=255)

                                                                                                       aaaCAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGAc                                              >  3:131314/1‑149 (MQ=255)

                                                                                                        aaCAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACg                                             >  1:605866/1‑149 (MQ=255)

                                                                                                         aCAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGa                                            >  2:11275/1‑149 (MQ=255)

                                                                                                              aaCGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGttt                                       <  3:6741/149‑1 (MQ=255)

                                                                                                               aCGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttc                                      <  4:180674/149‑1 (MQ=255)

                                                                                                                 gCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcct                                    <  2:614891/149‑1 (MQ=255)

                                                                                                                 gCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcct                                    >  1:430274/1‑149 (MQ=255)

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

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

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

                                                                                                                           gAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTCGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTaa                           >  1:289792/1‑149 (MQ=255)

                                                                                                                           gAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTaa                           >  2:285064/1‑149 (MQ=255)

                                                                                                                             aTGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACa                         >  4:485739/1‑149 (MQ=255)

                                                                                                                              tGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACaa                        <  4:554680/149‑1 (MQ=255)

                                                                                                                                gACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGc                      <  2:106911/149‑1 (MQ=255)

                                                                                                                                 aCTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCg                     <  2:452188/149‑1 (MQ=255)

                                                                                                                                       aTGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcc                                     >  3:351104/1‑127 (MQ=255)

                                                                                                                                       aTGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcc                                     <  4:389014/127‑1 (MQ=255)

                                                                                                                                       aTGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcc                                     <  4:351104/127‑1 (MQ=255)

                                                                                                                                       aTGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcc                                     >  3:389014/1‑127 (MQ=255)

                                                                                                                                         gcttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCa             >  3:116558/5‑149 (MQ=255)

                                                                                                                                          cttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTAtctc                                                      >  3:174301/4‑107 (MQ=255)

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

                                                                                                                                              gCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTAc        <  3:517105/149‑1 (MQ=255)

                                                                                                                                                ctctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTAcgc      >  4:318771/1‑149 (MQ=255)

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

                                                                                                                                                    gTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCtt  <  1:257761/149‑1 (MQ=255)

                                                                                                                                                     tAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCaa            <  4:481928/138‑1 (MQ=255)

                                                                                                                                                     tAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCaa            >  3:481928/1‑138 (MQ=255)

TTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCA‑TTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCGCTT  >  NC_000913/4573041‑4573336

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