BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,412,796	G→T	11.0%	S137I (AGC→ATC)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,796	0	G	T	11.0%	1496.7 / 27.0	495	S137I (AGC→ATC)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base G (135/288); new base T (65/2); total (205/290)
Fisher's exact test for biased strand distribution p-value = 3.57e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 500 total aligned reads displayed.

GCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCT  >  NC_000913/4412648‑4412945

gcgcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGc                                                                                                                                                       <  1:8442845/150‑1 (MQ=255)

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

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

       gCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTGTGGGTTTTTTTTTAGACGGGTTa                                                                                                                                                >  1:2727230/1‑150 (MQ=255)

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

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

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

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

               aacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTaaa                                                                                                                                        <  1:4909985/150‑1 (MQ=255)

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

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

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

                     aaaaCCGCATTTTAAGCCGATGTTCCGAAGAATATGACTAAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTTGATCCCGGTCGT‑TTATTGATTTTTTTTAAGCCGGGTTATTAGTTAAGATagc                                                                                                                                  >  2:2418781/1‑150 (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)

                            cATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑c                                                                                                                          <  1:2259747/150‑1 (MQ=255)

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

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

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

                                ttATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTGTGGGTTTTTTTTTATCCGGGGTTTTGGGAAATATAGCGATGAG‑CCggt                                                                                                                      >  1:2847892/1‑148 (MQ=255)

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

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

                                   tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCgtttgt                                                                                                                   <  2:1460123/150‑1 (MQ=255)

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

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

                                       gATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTtat                                                                                                               <  1:2040242/150‑1 (MQ=255)

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

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

                                            tCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAg                                                                                                          <  1:1737237/150‑1 (MQ=255)

                                            tCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAg                                                                                                          <  1:1878264/150‑1 (MQ=255)

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

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

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

                                                     tatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTc                                                                                                 <  2:3113633/150‑1 (MQ=255)

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

                                                           tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTGTTGATTTTTTTTAATCCGGGGTATTAGTAAAAATAGCGATGAT‑CCGGTTTTTTATAAAAGCGGCGGTTACTCatt                                                                                           >  2:1979848/1‑150 (MQ=255)

                                                           tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCCAAGAG‑CCGTTTGTTTATAACACCGGCGGTTCCTCttt                                                                                           >  2:8667177/1‑150 (MQ=255)

                                                           tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTATGGAATTTTTTTAATCCGGGGTATTAGTAAAGAAAGCGATGTG‑CCGTTTTTTTATAAAACACGCGGTTACTCatt                                                                                           >  2:551218/1‑150 (MQ=255)

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

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

                                                            cagcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCtttt                                                                                          <  1:122136/150‑1 (MQ=255)

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

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

                                                              gcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGGTGAACAGAGTTGGATCCCGGTCGG‑TTCGGGATGTTTGTTAAGCCGGGGTATTCGTAAAGATCGCGATGAG‑CCGGTTGTTTATAACAGCGGCGGTGGCTCTTGTGc                                                                                        >  2:7263449/1‑150 (MQ=255)

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

                                                                    cTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgcg                                                                                  <  2:5626585/150‑1 (MQ=255)

                                                                    cTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgcg                                                                                  <  1:139755/150‑1 (MQ=255)

                                                                     tGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTATTGATTTTTTTTTAACCGGGGTTTTGGGAAAGATAGCGCGGAT‑GCGGTTGTTTTTAAAACCGGCGGGTGCTCATTTTGTTAgcgc                                                                                 >  2:3413866/1‑150 (MQ=255)

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

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

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

                                                                              aaGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTGTGGGTTTTTTTTTATGCGGGGTTTTGTTAAAGATGGCGGTTGG‑CCGGGTGTTTATAAAAGCGGGGGTTGCTCGTTTTGTTGGCGCCCCGTCAgg                                                                        >  2:1754711/1‑150 (MQ=255)

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

                                                                               aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgc                                                                       <  2:1314387/150‑1 (MQ=255)

                                                                               aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgc                                                                       <  1:3266771/150‑1 (MQ=255)

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

                                                                                 taagAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGgcgc                                                                     <  1:159188/147‑1 (MQ=255)

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

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

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

                                                                                   cGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                                   <  2:8552014/150‑1 (MQ=255)

                                                                                            ctactaCGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTgg                                                          <  1:1123941/150‑1 (MQ=255)

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

                                                                                             tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                                         <  1:6697314/150‑1 (MQ=255)

                                                                                               ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                                       <  1:2450954/150‑1 (MQ=255)

                                                                                               ctacGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                                       <  2:6087530/150‑1 (MQ=255)

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

                                                                                                tacGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTGTGGATTTTTTTTAAGCCCGGGTATTCGTAAATATCGCGATGAG‑CCGGTTTTTTATAAAAACGGCGGTTACACGTTTTCATAGCGCCCAGTATGGCGCCGCACATATGGTGag                                                      >  2:469035/1‑149 (MQ=255)

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

                                                                                                   gACAATGTTGAACAGAGTTGGATCCCGGTCGT‑GTCTGGATTTTTTTTAAGCCGGGGTATTAGGAAAGATAGCGCTGAG‑CCGTTTTTTTATAAAAACGGCGGTTCCTCTTTTTCTTCGCCCCCCCGCTGGCGCCGCAAAGCAGGTGATTaa                                                   >  2:1810865/1‑149 (MQ=255)

                                                                                                        tGTTGAACAGAGTTGGATCCCGGTCGT‑GTGTGGATTTTTTTTAATGCGGGTTTTTGGTAAATATCGCGTTTTG‑CCCTTTTTTTTTAAAAGAGGCGGTTTCAGTTTTTGGTTGCGCCCCCGCGGGCGCGGCACATGTGGTGttttaacagg                                              >  2:3297136/1‑143 (MQ=255)

                                                                                                              aCAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTc                                        <  1:2550535/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                      ggATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                                <  1:3057424/150‑1 (MQ=255)

                                                                                                                           ccGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaa                           >  2:6412608/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                         gATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTAct              >  1:6453361/1‑150 (MQ=255)

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

                                                                                                                                          aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc             >  2:8975869/1‑150 (MQ=255)

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

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

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

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

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

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

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

                                                                                                                                                   aaGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcgg    <  1:3171600/150‑1 (MQ=255)

                                                                                                                                                   aaGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcgg    <  2:4515119/150‑1 (MQ=255)

                                                                                                                                                    aGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcggc   >  2:9075229/1‑150 (MQ=255)

                                                                                                                                                     gCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTGCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTGCTGGCGCTCCGGTTGTGGCAAAAAAGAGTCTTTATCTTCGGGCGGCt  >  1:6273559/1‑150 (MQ=255)

GCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGT‑TTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAG‑CCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCT  >  NC_000913/4412648‑4412945

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