BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,412,794	A→T	9.1%	L136F (TTA→TTT)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,794	0	A	T	9.1%	1490.8 / 29.3	479	L136F (TTA→TTT)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base A (148/280); new base T (49/2); total (197/282)
Fisher's exact test for biased strand distribution p-value = 1.28e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 488 total aligned reads displayed.

GCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGG  >  NC_000913/4412646‑4412943

    gcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCg                                                                                                                                                  <  1:4488376/150‑1 (MQ=255)

     cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCgg                                                                                                                                                 <  2:3238239/150‑1 (MQ=255)

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

      gTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCggg                                                                                                                                                <  1:380729/150‑1 (MQ=255)

           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTAtt                                                                                                                                           <  2:9253395/150‑1 (MQ=255)

           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTAtt                                                                                                                                           <  2:8522069/150‑1 (MQ=255)

           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTGTGGATTTTTTTTTAGCCGGGTtttt                                                                                                                                           >  2:2625553/1‑150 (MQ=255)

            ggAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTc                                                                                                                                          <  1:2007826/150‑1 (MQ=255)

            ggAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTc                                                                                                                                          <  2:1130821/150‑1 (MQ=255)

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

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

                 aacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTTGGTTTTTTTTTATACCGGG‑TATTCggtaa                                                                                                                                     >  1:8383297/1‑146 (MQ=255)

                 aacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTCTGGATTTTTTTTTAGCCGGGGTTTTAGTaaa                                                                                                                                     >  2:2763779/1‑150 (MQ=255)

                 aacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTCGGGATTTTTGTTAAGCCGGGTTCTTCGTaaa                                                                                                                                     >  2:772829/1‑150 (MQ=255)

                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGAt                                                                                                                                  <  2:3679769/150‑1 (MQ=255)

                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTCTGGATTTTTTTTAATACGGGGTGTTCGTAAatat                                                                                                                                  >  2:2043807/1‑150 (MQ=255)

                      cAAAACTGCATTTTATGCCGATGTTCCTAAGAATATGACTCATCAGAACTGTCAGTCTAAGTACGTAAATTGCTACTACGTCAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATcg                                                                                                                                <  2:1458105/150‑1 (MQ=255)

                       aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATcgc                                                                                                                               <  1:1355512/150‑1 (MQ=255)

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

                           cGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTTGGATTTTT‑TTAATACAGGTTATTAGTAAATAAAGCGATg                                                                                                                           >  1:9295484/1‑150 (MQ=255)

                           cGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTTTTAAGCCGGGGTATTCGTAAAGATAGCGATg                                                                                                                           >  2:6975481/1‑150 (MQ=255)

                            ggCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTCTGGATTTTTGTTAATCCGGGTTTTTCGTTAAGATAGCGATGa                                                                                                                          >  1:7218662/1‑150 (MQ=255)

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

                               aTTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGcc                                                                                                                       <  2:8306322/150‑1 (MQ=255)

                                ttttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTGTTGTTTTTTTTTAAGACGGGTTATTGGTAAAGAAAGAGATGAGCCg                                                                                                                      >  1:1010515/1‑150 (MQ=255)

                                     tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtttgt                                                                                                                 <  2:4646620/150‑1 (MQ=255)

                                     tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtttgt                                                                                                                 <  1:7649347/150‑1 (MQ=255)

                                      gCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCgtttgtt                                                                                                                <  2:542583/150‑1 (MQ=255)

                                          aTGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATAGCGATGAGCCGGTTGTTtaaa                                                                                                            >  2:2421850/1‑148 (MQ=255)

                                             ttCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACa                                                                                                         <  2:4885414/150‑1 (MQ=255)

                                             ttCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACa                                                                                                         <  1:5461315/150‑1 (MQ=255)

                                              tCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAg                                                                                                        <  1:1878207/150‑1 (MQ=255)

                                                cGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAgcg                                                                                                      <  1:8835162/150‑1 (MQ=255)

                                                  aagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAgcggc                                                                                                    <  2:9345216/150‑1 (MQ=255)

                                                  aagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTCTGGATTTTTTTTAAGCCGGGGGATTCGGAAAAATCCCCCAGAGCCCGTTGTTGATAAAAgccgc                                                                                                    >  1:5923575/1‑150 (MQ=255)

                                                           aCTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCctct                                                                                           <  2:8278206/150‑1 (MQ=255)

                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTTAGCCGGGTTATTCGTAAAGATAGCGATGTGCCGTTTTTTTATAACAGCGGCGGTTCCTCttt                                                                                         >  2:8795885/1‑150 (MQ=255)

                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTAGTAAAGATAGCGAAGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCttt                                                                                         >  2:6380649/1‑150 (MQ=255)

                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTGTGGGTTTTTTTTTATCCGGGGTATTAGGAAATATGGCGATGTGGCGGTTGTTTTTAAAAGCGGCGGTTGctaatt                                                                                         >  1:9605807/1‑146 (MQ=255)

                                                              caggagAACTGTCAGTCTAAGTCCGAAAAGTGCTCCTACGACCATGTCGACCAGAGTTGGATCCCGGTCGT‑GTCTGGATTTTTGTTTAGCGGGGTTTTTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCtttt                                                                                        <  1:5816496/150‑1 (MQ=255)

                                                              cagcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTGTGGATTTTTGTTTAACCGGGGTATTAGTAAAAATAGCGCTTATCCGTTTGTTTTTAAAAGCGGCGGGTACTCgttt                                                                                        >  2:1748247/1‑150 (MQ=255)

                                                              cagcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCtttt                                                                                        <  2:7476976/150‑1 (MQ=255)

                                                                        gTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGcc                                                                              <  1:3256873/150‑1 (MQ=255)

                                                                         tCAGTCTAAGTACTAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGccc                                                                             <  1:6747818/150‑1 (MQ=255)

                                                                         tCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGccc                                                                             <  1:5095673/150‑1 (MQ=255)

                                                                         tCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTCTGGTTTTTTTTTTATCCGGGGGGTTATGAAAGATCGCGATGAGCCGGGTGTTTTTTACAAAGGCGGGTGCTCATTTGGTTCGCGccc                                                                             >  1:9534934/1‑150 (MQ=255)

                                                                          cAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCa                                                                            <  2:2536933/150‑1 (MQ=255)

                                                                            gTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGCTTTTTTTTAAGCCGGGGTATTCGTAAAGATCGCGATGATCCGGTTTTTTATAAAAGCGGCGGTTCCTCATTTGCTTCGCGCCCAGt                                                                          >  2:4066650/1‑150 (MQ=255)

                                                                               tAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTg                                                                       <  1:3445408/150‑1 (MQ=255)

                                                                               tAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTGTGGTTTTTTTTTTATCCGGGGTATTAGTAAAGATAGCGCTGAGCCGGTTTTTTTTAAAAACGGCGGGTGCACATTTGGTTGGCGCCCCCTCTg                                                                       >  1:5278589/1‑150 (MQ=255)

                                                                                 aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgc                                                                     <  1:439016/150‑1 (MQ=255)

                                                                                  gTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgcg                                                                    <  2:581440/150‑1 (MQ=255)

                                                                                   tACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgcgc                                                                   <  1:8092876/150‑1 (MQ=255)

                                                                                    aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTTTTAATCCCGGTTATTCGTAAAGATCGCGATTAGCCGGTTGTTTATAACAGCGGCGGGTCCTCATTTTCTTCGCGCCCACTCTGGCgca                                                                  >  2:7928964/1‑150 (MQ=255)

                                                                                    aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                                  <  1:1729839/150‑1 (MQ=255)

                                                                                     cGAAAATTGCTACTACGACACTGTTGAACAGAGGTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                                 <  2:7144099/150‑1 (MQ=255)

                                                                                     cGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                                 <  2:197894/150‑1 (MQ=255)

                                                                                     agAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                                 <  2:8551881/149‑1 (MQ=255)

                                                                                       aaaaTTGCTACGACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcagca                                                               <  2:6078227/150‑1 (MQ=255)

                                                                                             gCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTg                                                         <  2:1272330/150‑1 (MQ=255)

                                                                                              ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg                                                        <  1:1123937/150‑1 (MQ=255)

                                                                                              ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg                                                        <  2:86631/150‑1 (MQ=255)

                                                                                              ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg                                                        <  1:2687870/150‑1 (MQ=255)

                                                                                               tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                                       <  1:6302511/150‑1 (MQ=255)

                                                                                               tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                                       <  2:9543531/150‑1 (MQ=255)

                                                                                                    cGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTa                                                  <  1:1192059/150‑1 (MQ=255)

                                                                                                     gACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGAATTTTGTTAAGCCGGGGTATTCGTAAAGATCGCCATGAGCCGTTTTTTTATAAAAGCGGCGGTTACTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCGGGTGAttta                                                 >  2:670405/1‑148 (MQ=255)

                                                                                                           gTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTgg                                           <  1:303690/150‑1 (MQ=255)

                                                                                                           gTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTgg                                           <  1:5050407/150‑1 (MQ=255)

                                                                                                               aaCAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCt                                       <  2:5517440/150‑1 (MQ=255)

                                                                                                                 cAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGGTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCAGGTGATTACTCCAGGCGCTcc                                     >  2:5648854/1‑150 (MQ=255)

                                                                                                                  agagTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCg                                    <  2:2513076/150‑1 (MQ=255)

                                                                                                                  agagTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCg                                    <  1:9036910/150‑1 (MQ=255)

                                                                                                                      ttGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTc                                <  2:3089387/150‑1 (MQ=255)

                                                                                                                        ggATCCCGTTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                              <  2:5182599/150‑1 (MQ=255)

                                                                                                                        ggATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                              <  2:178660/150‑1 (MQ=255)

                                                                                                                        ggATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                              <  2:3400266/150‑1 (MQ=255)

                                                                                                                        ggATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                              <  1:8697752/150‑1 (MQ=255)

                                                                                                                        ggATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                              <  1:284671/150‑1 (MQ=255)

                                                                                                                              cGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaa                        <  1:2158630/150‑1 (MQ=255)

                                                                                                                                  cGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga                    <  1:8366350/150‑1 (MQ=255)

                                                                                                                                  cGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga                    <  2:7270525/150‑1 (MQ=255)

                                                                                                                                   cggttCTGGATTTGTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGGTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga                    <  2:5029784/147‑1 (MQ=255)

                                                                                                                                    t‑ttCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGt                  <  1:2546759/150‑1 (MQ=255)

                                                                                                                                       tCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCt                <  1:7764615/150‑1 (MQ=255)

                                                                                                                                         tGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTa              <  2:8383297/150‑1 (MQ=255)

                                                                                                                                            aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc           >  2:8975816/1‑150 (MQ=255)

                                                                                                                                            aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc           >  1:2430238/1‑150 (MQ=255)

                                                                                                                                            aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc           >  1:6817480/1‑150 (MQ=255)

                                                                                                                                             tttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActct          >  2:4294033/1‑150 (MQ=255)

                                                                                                                                              ttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTATTTTGCTTAGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGGTCCGGTTATTGCAAAAAAGAGTGTTACTCtt         >  1:4665766/1‑150 (MQ=255)

                                                                                                                                                ttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCg       <  2:4441830/150‑1 (MQ=255)

                                                                                                                                                ttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGC‑AAAAAGAGTCTTACTCTTCgg      >  1:7425636/1‑150 (MQ=255)

                                                                                                                                                 tGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCgg      <  1:4169332/150‑1 (MQ=255)

                                                                                                                                                  gTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCggg     <  2:5740851/150‑1 (MQ=255)

                                                                                                                                                   ttAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggc    >  1:397330/1‑150 (MQ=255)

                                                                                                                                                   ttAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggc    >  2:8138037/1‑150 (MQ=255)

                                                                                                                                                   ttAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTAGCGCCCAGTCTGGCGCAGCACATCTGGTGATTAATCCTGGCGCTCCGGTTCTGGCAAAAAAGATTCTTAATCTTCGggc    >  1:8126482/1‑150 (MQ=255)

                                                                                                                                                    tAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcg   >  1:1568419/1‑150 (MQ=255)

                                                                                                                                                     aaGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcgg  <  1:3171522/150‑1 (MQ=255)

GCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGG  >  NC_000913/4412646‑4412943

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: `‑`