BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,412,779	T→G	19.1%	V131V (GTT→GTG)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,779	0	T	G	19.1%	1115.7 / 126.2	484	V131V (GTT→GTG)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base T (86/298); new base G (98/2); total (184/300)
Fisher's exact test for biased strand distribution p-value = 0.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 499 total aligned reads displayed.

ATAATGCCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTC  >  NC_000913/4412631‑4412928

ataatGCCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGttt                                                                                                                                                       <  1:4899279/150‑1 (MQ=255)

  aatGCCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCt                                                                                                                                                     <  2:6996251/150‑1 (MQ=255)

      ccGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGAt                                                                                                                                                 <  1:2304214/150‑1 (MQ=255)

      ccGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTATTGAt                                                                                                                                                 >  2:335766/1‑150 (MQ=255)

      ccGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTATTGAt                                                                                                                                                 >  1:3254495/1‑150 (MQ=255)

      ccGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGAt                                                                                                                                                 >  1:197894/1‑150 (MQ=255)

      ccGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGAt                                                                                                                                                 >  1:9132989/1‑150 (MQ=255)

        gACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGGGGAttt                                                                                                                                               >  1:4618485/1‑150 (MQ=255)

        gACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGAttt                                                                                                                                               <  1:3560030/150‑1 (MQ=255)

         aCATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGAtttt                                                                                                                                              <  2:6847810/150‑1 (MQ=255)

          cATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGAttttt                                                                                                                                             <  1:144705/150‑1 (MQ=255)

          cATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGAttttt                                                                                                                                             <  1:4273011/150‑1 (MQ=255)

            tttGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGt                                                                                                                                           >  2:5497429/1‑150 (MQ=255)

              tGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGTGGAttttttttt                                                                                                                                         >  2:113276/1‑149 (MQ=255)

               gcgcgcGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTaa                                                                                                                                        <  2:5931520/150‑1 (MQ=255)

                     gTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAATCCggg                                                                                                                                  >  1:6455713/1‑150 (MQ=255)

                     gTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTTGATTTTTTTTAAGGCggg                                                                                                                                  >  2:8366469/1‑150 (MQ=255)

                        gCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTAAGCCGGGTTa                                                                                                                               >  2:6423537/1‑150 (MQ=255)

                          tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGGTTTTTGTTAAGCCgggggatg                                                                                                                             >  2:4577953/1‑145 (MQ=255)

                          tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTTTTAAGCCGGGTTAtt                                                                                                                             >  1:3976934/1‑150 (MQ=255)

                             aacaacGCCAAAACGGCATTTTACGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGt                                                                                                                          <  1:5517791/150‑1 (MQ=255)

                              acaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTa                                                                                                                         <  2:9066455/150‑1 (MQ=255)

                              acaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTa                                                                                                                         <  2:8306303/150‑1 (MQ=255)

                                aacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTaaa                                                                                                                       <  2:734512/150‑1 (MQ=255)

                                   gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGAt                                                                                                                    <  2:3679678/150‑1 (MQ=255)

                                    ccAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAATACGGGGTATTCGTAAatatg                                                                                                                   >  2:2860238/1‑149 (MQ=255)

                                      aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATcgc                                                                                                                 <  1:6175963/150‑1 (MQ=255)

                                      aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATcgc                                                                                                                 <  1:8876959/150‑1 (MQ=255)

                                      aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATcgc                                                                                                                 <  1:591946/150‑1 (MQ=255)

                                       aanCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGGGGATTTTTGTTAAGCCCGGGTAGTCGTCAAAATcgcg                                                                                                                >  1:7144099/1‑150 (MQ=255)

                                         aCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGAt                                                                                                              <  1:8353970/150‑1 (MQ=255)

                                           ggCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGa                                                                                                            <  1:1146696/150‑1 (MQ=255)

                                           ggCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGa                                                                                                            <  2:7833922/150‑1 (MQ=255)

                                           ggCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATAGCGAtta                                                                                                            >  2:7185077/1‑148 (MQ=255)

                                            gCATTTTTTGCATATGTTACGAAGAATATGAATCAGGTGTTCTGTCAGTATAAGTACGAAAATTGCTACTACGACTATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATCCGGGTTTTTCGTAAAGATCGCGTTGAg                                                                                                           <  1:5685459/150‑1 (MQ=255)

                                            gCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAAGCCGGGTTATTAGTAAAAATAGCGATGAg                                                                                                           >  2:5037987/1‑150 (MQ=255)

                                             cATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGATTTTTTTTTAGACCGGTGTTTGGGAAAGATCGCGATGAG‑c                                                                                                         >  1:2730844/1‑150 (MQ=255)

                                              aTTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑cc                                                                                                        <  2:6117047/150‑1 (MQ=255)

                                               ttttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTTGTTTTTTTTTAAGACGGGTTATTGGTAAAGAAAGAGATGAG‑CCg                                                                                                       >  1:1010515/1‑150 (MQ=255)

                                               ttttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTTATCAGGGTTATTAGTTAAGATCGCGATGAG‑ccc                                                                                                       >  2:8486127/1‑149 (MQ=255)

                                               ttttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTTTTTATCCGGGGTGTTGGGAAAGATAGCGCTGAggac                                                                                                        >  2:2546759/1‑147 (MQ=255)

                                                 ttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTGGTGGGTTTTTTTTTATCCGGGTTATTAGGGAAGATAGAGAGGAG‑CCgtt                                                                                                     >  1:2848887/1‑150 (MQ=255)

                                                  tATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCgttt                                                                                                    <  1:1352201/150‑1 (MQ=255)

                                                    tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCgtttgt                                                                                                  <  1:5248170/150‑1 (MQ=255)

                                                    tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCgtttgt                                                                                                  <  1:7649186/150‑1 (MQ=255)

                                                    tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCgtttgt                                                                                                  <  2:1595484/150‑1 (MQ=255)

                                                      ccGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGGGGATTTTTTTTTATCCCGGGTATTCGTAAAAATAGCGCTGTG‑CCgtttgttt                                                                                                >  2:4295624/1‑150 (MQ=255)

                                                         aTGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATAGCGATGAG‑CCGGTTGTTtaaa                                                                                             >  2:2421850/1‑148 (MQ=255)

                                                             tCCGAAGAATATGACACACCAGAACTGACAGTCTAAGTACGAAAATTGCTACTACCACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCAGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAg                                                                                         <  2:8088756/150‑1 (MQ=255)

                                                               cGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGCTTTTTGTTTATCCGGGGTATTCGTAAAGATAGCGATTAT‑CCGTTTGTTTTTaaaatcg                                                                                       >  2:3526467/1‑145 (MQ=255)

                                                                 aagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAgcggc                                                                                     <  2:9345216/150‑1 (MQ=255)

                                                                   gaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGCGGGATTTTTTTTTAGCCCGGGTATTCGGGAAGATCGCGCTGGG‑CCCGTTGTTTATAAAAGcggcgg                                                                                   >  1:3059276/1‑150 (MQ=255)

                                                                    aaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGGGGGTTTTTTTTTATCCCGGGGTTTAGTTAAGATAGCGATGTGCCCGTTTGTTTTTAA‑AAcggcggg                                                                                  >  1:9688/1‑149 (MQ=255)

                                                                          aCTCAGCAGTACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCGCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCctct                                                                            <  2:8277956/150‑1 (MQ=255)

                                                                              agcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTg                                                                        <  1:5083934/150‑1 (MQ=255)

                                                                              agcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTTGATTTTTTTTTATGCGGGGTTTTGGTTAAGATAGCGATGGG‑CCGGTTGTTTTTAAAAACCGCGGGTTCGCTTTTg                                                                        >  1:8102408/1‑150 (MQ=255)

                                                                                 agATCTGTCAGTCTAATTACTATTATTGCTACTACTACAATGTTTAACAGAGTTTGATCCCGGTTGTTTCTTGATTTTTGTTAAGCCGGGTTATTCGTTAAGACCGCGATGAG‑CCGTTTGTTTATAACAGCTGCGGTTCCTCTTTTGCtt                                                                     <  1:3438439/150‑1 (MQ=255)

                                                                                 agAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGAtt                                                                     >  2:5615896/1‑150 (MQ=255)

                                                                                 agAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTATGGATTTTTGTTAATCCGGGTTATTAGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTATTTTGCtt                                                                     >  2:5691649/1‑150 (MQ=255)

                                                                                   aaCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGATTTTTGTTAAGCCCGGTTATTAGTAAAGATCGCGATTAT‑CCGTTTGTTTATAACAGCGGCGGTTACACTTTTGCGTcg                                                                   >  2:4920278/1‑150 (MQ=255)

                                                                                   aaCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGATTTTTGTTAAGCCGGGGTATTCGTAAAGATAGCGATGAGCCCGTTTG‑TTATAAAAGCCGCGGGTACGCTTTTGcatag                                                                   >  2:263275/1‑146 (MQ=255)

                                                                                    aCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgc                                                                  <  1:4437085/150‑1 (MQ=255)

                                                                                       gTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGcc                                                               <  1:9459434/150‑1 (MQ=255)

                                                                                        tCAGTCTAAGTACTAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGccc                                                              <  1:6747818/150‑1 (MQ=255)

                                                                                            tCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTc                                                          <  1:1631267/150‑1 (MQ=255)

                                                                                              tAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTg                                                        <  2:1778160/150‑1 (MQ=255)

                                                                                                aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgc                                                      <  1:5414326/150‑1 (MQ=255)

                                                                                                aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgc                                                      <  2:9256318/150‑1 (MQ=255)

                                                                                                aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgc                                                      <  1:210710/150‑1 (MQ=255)

                                                                                                 gTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgcg                                                     <  2:581184/150‑1 (MQ=255)

                                                                                                  tACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgcgc                                                    <  1:8092842/150‑1 (MQ=255)

                                                                                                   aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGGTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAAAAGCGGCGGTTACTCGTTTGCGTCGCGCCCAGTCGGGGgca                                                   >  2:5423523/1‑150 (MQ=255)

                                                                                                    cGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                  <  2:197894/150‑1 (MQ=255)

                                                                                                    agAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                  <  2:8551881/149‑1 (MQ=255)

                                                                                                           tGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCt                                           <  1:1778036/150‑1 (MQ=255)

                                                                                                            gCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTg                                          <  2:1272330/150‑1 (MQ=255)

                                                                                                             ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg                                         <  1:1123937/150‑1 (MQ=255)

                                                                                                             ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg                                         <  2:86631/150‑1 (MQ=255)

                                                                                                              tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                        <  2:9543531/150‑1 (MQ=255)

                                                                                                              tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                        <  1:8765946/150‑1 (MQ=255)

                                                                                                                ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                      <  1:2451291/150‑1 (MQ=255)

                                                                                                                ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                      <  1:2660648/150‑1 (MQ=255)

                                                                                                                ctacGACAACGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGTTTTTTTTTTATGCGGGTTTTTGGTTAATTTGGGGATGTT‑GCGTTTTTTTTTATAATCGGCGGTTTGTATTTTTGTTTGCGCGCAGTGTTGCGCATGACttgttgttt                                      >  2:8836194/1‑143 (MQ=255)

                                                                                                                        aTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCt                              <  2:470984/150‑1 (MQ=255)

                                                                                                                        aTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCt                              <  2:3998052/150‑1 (MQ=255)

                                                                                                                         tGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTg                             <  1:8486127/150‑1 (MQ=255)

                                                                                                                          gTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTgg                            <  1:5050407/150‑1 (MQ=255)

                                                                                                                              aaCAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCt                        <  2:1827084/150‑1 (MQ=255)

                                                                                                                                 agagTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCg                     <  2:2513076/150‑1 (MQ=255)

                                                                                                                                     ttGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTc                 <  2:226148/150‑1 (MQ=255)

                                                                                                                                     ttGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTc                 <  2:3089387/150‑1 (MQ=255)

                                                                                                                                       ggATCCCGTTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg               <  2:5182599/150‑1 (MQ=255)

                                                                                                                                       ggATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg               <  1:7740109/150‑1 (MQ=255)

                                                                                                                                       ggATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg               <  1:7417725/150‑1 (MQ=255)

                                                                                                                                         aTCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGc             <  2:1446150/150‑1 (MQ=255)

                                                                                                                                            ccGGTCGTTTATGGATTTTTGTTAATCCGGGTTATTCGTAAATATCGCGATGAT‑CCGTTTGTTTATAACAGCGGCGGTTACTCTTTTGCTTCGCGCCCAGTCTGGCGCCGCACTTCTGGTGATTACTCCTGGCGCTCCGGGTCTTGCaaa          >  2:3230139/1‑150 (MQ=255)

                                                                                                                                              ggTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaaa        <  1:5800408/150‑1 (MQ=255)

                                                                                                                                               gtagTTTCTGGATTTTTTTTTAGCCGGTTTATTCGTAAAGATCGCGATGAG‑CATTTTGTTTATACCAGCGGCGGTTCCTCTTTTGCTTCGCGCCCTGTCTGGCGCAGCTCTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaaaa       <  1:3901564/147‑1 (MQ=255)

                                                                                                                                                 cGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga     <  2:7607399/150‑1 (MQ=255)

                                                                                                                                                 cGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga     <  1:8366385/150‑1 (MQ=255)

ATAATGCCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTC  >  NC_000913/4412631‑4412928

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`