BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,054,888	A→T	15.1%	V186E (GTA→GAA)	ygfB ←	UPF0149 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,054,888	0	A	T	15.1%	466.6 / 16.2	176	V186E (GTA→GAA)	ygfB	UPF0149 family protein
Reads supporting (aligned to +/- strand): ref base A (52/95); new base T (27/2); total (79/97)
Fisher's exact test for biased strand distribution p-value = 7.19e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037

tACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTa                                                                                                                                                       <  2:1285067/150‑1 (MQ=255)

tACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTa                                                                                                                                                       <  1:5915/150‑1 (MQ=255)

 aCTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAc                                                                                                                                                      <  2:263722/150‑1 (MQ=255)

   ttCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACtt                                                                                                                                                    >  2:149691/1‑150 (MQ=255)

    tCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTc                                                                                                                                                   <  1:223543/150‑1 (MQ=255)

     cTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCt                                                                                                                                                  <  1:3817126/150‑1 (MQ=255)

      tGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTGTTTTGTTTTTAGTGTAGAGTCGGGTGTTGTTCGTCTg                                                                                                                                                 >  2:2542100/1‑150 (MQ=255)

        gTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGgc                                                                                                                                               <  1:1702474/150‑1 (MQ=255)

        gTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGgc                                                                                                                                               >  1:609286/1‑150 (MQ=255)

           cagcaAAAATCAGCTCGGAGCTGCCGGTTTGCATTTGCTCCACCATGGCCTGACGTCGACGCTGAAACTCTTGCCGTTATATCTCAATCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGATTCGGTTTTTGTACTTCTGGcgcg                                                                                                                                            <  2:4183/150‑1 (MQ=255)

           cagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTGCGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTCTTCTGGggcg                                                                                                                                            >  1:1657479/1‑150 (MQ=255)

            agcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCGTATGGCGcgg                                                                                                                                           >  1:1520803/1‑150 (MQ=255)

             gcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggt                                                                                                                                          <  1:1151108/150‑1 (MQ=255)

             gcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggt                                                                                                                                          <  1:2607239/150‑1 (MQ=255)

              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTGGTCTGGCGcgggg                                                                                                                                         >  1:1415699/1‑148 (MQ=255)

                  aaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTATTGTGGCGCGGGCGGtt                                                                                                                                     >  1:1758588/1‑150 (MQ=255)

                   aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTATGGCGCGGTCGGGTg                                                                                                                                    >  2:1614938/1‑150 (MQ=255)

                   aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGCGGTCGGTTg                                                                                                                                    >  1:3080800/1‑150 (MQ=255)

                     cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                  <  2:1179182/150‑1 (MQ=255)

                     cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                  <  1:4013023/150‑1 (MQ=255)

                     cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                  <  2:269872/150‑1 (MQ=255)

                     cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                  <  2:3244328/150‑1 (MQ=255)

                         gcggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTGTTGTTCTTCTTGCGCGGGCGGGTGGGGCTg                                                                                                                              >  2:766446/1‑150 (MQ=255)

                          cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  2:1401747/150‑1 (MQ=255)

                            gcgcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTATGGCGCGGTCGGTTGCGGATgtgt                                                                                                                           >  2:4079210/1‑150 (MQ=255)

                               cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTaaa                                                                                                                        <  1:3770601/150‑1 (MQ=255)

                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  1:976688/150‑1 (MQ=255)

                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  2:2209970/150‑1 (MQ=255)

                                tGCCGGGGTGCATGTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  2:3093039/150‑1 (MQ=255)

                                  ccGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGgt                                                                                                                     <  2:1609348/150‑1 (MQ=255)

                                  ccGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGgt                                                                                                                     <  1:1584982/150‑1 (MQ=255)

                                  ccGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGgt                                                                                                                     <  2:253917/150‑1 (MQ=255)

                                     gctgatttttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTc                                                                                                                  >  2:1729585/8‑150 (MQ=255)

                                      gTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTAGTCTGGCGCGGTCGGGTTGGGGTTAGTAAAGgtgtag                                                                                                                 >  1:1923168/1‑148 (MQ=255)

                                      gTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCg                                                                                                                 <  1:1721225/150‑1 (MQ=255)

                                      gTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTAAGGCGCGGTCGGTTGCGGCTGAGTAAAAGTGTCg                                                                                                                 >  1:905971/1‑150 (MQ=255)

                                        tGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTg                                                                                                               <  2:1082439/150‑1 (MQ=255)

                                            tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACATCAGGCGCGGTCGGTTGCGGATGTGTAAATGTGTCGTGGCAt                                                                                                           >  2:2552553/1‑150 (MQ=255)

                                             ttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATa                                                                                                          <  2:1772231/150‑1 (MQ=255)

                                              tGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATaa                                                                                                         <  2:242299/150‑1 (MQ=255)

                                                 tCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGGTGCGGCTGAGTCAAGGGGTCGTGGCATAAAAg                                                                                                      >  2:3968960/1‑150 (MQ=255)

                                                 tCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTCTTCTGGCGCGGGCGGGTGCGGGTTTGTAAAGGTGTGGTGGCCTAAAAg                                                                                                      >  2:745787/1‑150 (MQ=255)

                                                  ccaccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTCTGGCGCGGTCGGGTGCGGGTGAGTAAATGTGTCGTGGCATAAAAgc                                                                                                     >  2:2925368/1‑150 (MQ=255)

                                                  ccaccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTATTGCGCGGTCGGTTGCGGGTGAGTAAAGGTTTCGTGGCATAACAgc                                                                                                     >  2:2934534/1‑150 (MQ=255)

                                                      caGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTACTTCTGGCGCGGTCGGTTGCGGGTGAGTAAAGGTGTCGTGGCATAACAGCGCgg                                                                                                 >  1:1737722/1‑150 (MQ=255)

                                                        gggCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTACGTCTTGCGCGGTCGGGTGCGGATGTGTAAAGGTGTGGTGGCATAACAGCCCCGCa                                                                                               >  1:3482262/1‑150 (MQ=255)

                                                           ccTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGGTCTTCTGGCGCGGTCGGTTGCGGATGAGGTAAGGGGGCGTGGGGTAACAGCGCGGCAaca                                                                                            >  1:4147325/1‑150 (MQ=255)

                                                            cTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTTGCGCGGTCGGTTGCGGGTGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAacac                                                                                           >  2:562591/1‑150 (MQ=255)

                                                             tGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTCTTGGGCGGGCGGGTGCGGCAGAGGAAAGGTGTCGTGGCCTAACCGCCCGGCAacacc                                                                                          >  2:2483753/1‑149 (MQ=255)

                                                               aCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTGTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTAGTGGCCTAAAAGCGCGGCAAAAAGaa                                                                                        >  2:4089891/1‑150 (MQ=255)

                                                                cGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAAc                                                                                       <  1:2253118/150‑1 (MQ=255)

                                                                   cGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTTATTATGGCGCGGGCGGGTGCGGGGGTGTAAAGGTGTCGTGGGAGAAAATCCCGGCAACACCAAccaa                                                                                    >  1:3533226/1‑147 (MQ=255)

                                                                   cGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTCTTCTGGCGCGGTGGGGTGGGGGTGGTTAAAAGTGTGGTGGCGTAAAAGCGCGGCCacacccacgga                                                                                    >  1:1122901/1‑143 (MQ=255)

                                                                    gACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGGACTTCTGGCGCGGGCGGGTGCGGGTGAGTACAGGTGTCGTGGCATAACAGCGCCGCAACCCGAACGGAt                                                                                   >  2:3672257/1‑150 (MQ=255)

                                                                    gACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTTATTCTGGCGCGGGCGGGTGGGGGTGTGTGGATGTGTCGTGGCGTAAAACAGCGGCAACACCAACGGAt                                                                                   >  1:900524/1‑150 (MQ=255)

                                                                      cGCTGAAACTCGTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTc                                                                                 <  2:2387226/150‑1 (MQ=255)

                                                                        cTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTTCTTATTGCGCGGGCGGTTTCGGCTGGGTAAAGGTTTAGTGGCATAAAAGCGCGGCAACACGAACCTATTCga                                                                               >  1:913006/1‑150 (MQ=255)

                                                                         tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTATTATGTCGCGGGCGGTTGCGGGTGTGTTAAGGTGTGGGGGCATAAAAGCGCGGCAAAACCAACGTATTAgat                                                                              >  1:1033410/1‑150 (MQ=255)

                                                                         tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgat                                                                              >  2:2446725/1‑150 (MQ=255)

                                                                          gAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTTCTTATGGCGCGGTCGGTTGCGGGTGAGTAAAGGGTTGGTGGCATAAAAACCCGGCAACACCAACGAATTAgatg                                                                             >  1:58152/1‑150 (MQ=255)

                                                                           aaaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatga                                                                            <  1:1378684/150‑1 (MQ=255)

                                                                                ctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATctct                                                                       <  1:2162723/150‑1 (MQ=255)

                                                                                 ttGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTGTTGCGCGGTCGGTTTCGGATGGGTAAATGTTTGGTGGCATAACAGCCCGGCAACACGAACGGAATCGAGGATATCtt                                                                      >  2:3223377/1‑150 (MQ=255)

                                                                                 ttGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTGTGGCGCGGTGGGTTGCGGATGAGTTAAGGTGTGGTGGCATAAAAGCGCCGCAACACCAACGAATTAGATGATCACtt                                                                      >  1:2671004/1‑150 (MQ=255)

                                                                                 ttGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTTGGGCGGTCGGGTGGGGATGAGTAAAGGTTTGGTGGCATAAAAGCGCGGCCACACCAACCTAATAgatgatatatt                                                                      >  1:3085562/1‑145 (MQ=255)

                                                                                      gggATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGc                                                                 <  2:1962375/150‑1 (MQ=255)

                                                                                        gATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGa                                                               <  1:2630706/150‑1 (MQ=255)

                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTATTAAGGCGCGGTCGGGTTCGGGTGAGTAAAGGTGTCGTGGCATAAAAGCGCGGCAACACCAACGTATTAGCAGATCTAATCAAGCGACATTTc                                                         >  2:1361540/1‑150 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       <  2:894630/150‑1 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       <  1:4079220/150‑1 (MQ=255)

                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGCGGTCGGGTGCGGATGAGTAAAGGTGTCGGGGCATAACAGCGCGGCAACACGAACGTATTAGATGATATCTTCAAGCGCCATTTCaa                                                       >  2:1697605/1‑150 (MQ=255)

                                                                                                 ctcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAg                                                      <  1:2565545/150‑1 (MQ=255)

                                                                                                  tcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGc                                                     <  2:664234/150‑1 (MQ=255)

                                                                                                   caTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTGTGGCGCGGTGGGGTTGGGGTGGGTTAAGGTGTGGGGGGATAACAGCGCCGCCACACCAAACGATTCGATGATATCTTAAATAGACAATTaaatca                                                    >  2:1092913/1‑147 (MQ=255)

                                                                                                       accctctCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACGTTTCAAGCTCTTc                                                <  2:2833512/147‑1 (MQ=255)

                                                                                                             tcCTTACGTTTTTTGTTTTTAGTGTAGAGAACGTTTTTGTACGTCTGGCGAGCTAGGTTGCGGCTGAGTAAAGGGGTCGTGGCATAAAAGCGCGGCAACACGGACGGATTCGGTGATCTCGTCAAGCGGCATTTCAAGCTCTTACGGAtc                                          >  2:1359235/1‑150 (MQ=255)

                                                                                                                      tttttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTGTTCCTGAtcttcgtcttc                                 <  2:1779287/150‑1 (MQ=255)

                                                                                                                        tttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATcttcgtcttcgt                               <  1:24884/150‑1 (MQ=255)

                                                                                                                         ttGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCttcgtcttcgtc                              <  1:313604/150‑1 (MQ=255)

                                                                                                                              tttAGTGTAGAGTCGGTTTTTTTACTTAAGGCGCGGTAGGTTTCGGATGAGTAAATGTGTAGTGGCAAAAAACCCCGGCAACACCAACCGATTCCATGATCTCGTAAACCCCCAATTCAAGCGATTACAGATATTCGTCGTACTCGTAAc                         >  2:874284/1‑150 (MQ=255)

                                                                                                                              tttAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAAc                         <  1:4068216/150‑1 (MQ=255)

                                                                                                                                tAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAAccc                       <  1:1968791/150‑1 (MQ=255)

                                                                                                                                tAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAAccc                       <  1:3008863/150‑1 (MQ=255)

                                                                                                                                tAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAAccc                       <  1:3594566/150‑1 (MQ=255)

                                                                                                                                 aGTGTAGAGTCGGTTTTTGTTCTTGTGGCGCGGTCGGTTGCGGTTGTGTATAGGTGTCGTGGCATAACAGCGCGGCAACACTATCGTTTTCTATGATCTCTTGAAGCGTCATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCa                      <  1:1694632/150‑1 (MQ=255)

                                                                                                                                 aGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACCCGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCa                      <  2:481042/150‑1 (MQ=255)

                                                                                                                                 aGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCa                      <  2:3616122/150‑1 (MQ=255)

                                                                                                                                  gtgtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAg                     <  1:3806408/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTGTACTTCTGTCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  1:3754785/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  1:3284957/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  2:411536/150‑1 (MQ=255)

                                                                                                                                    gtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  1:2925368/150‑1 (MQ=255)

                                                                                                                                        agTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTgcgc               >  2:3266239/1‑150 (MQ=255)

                                                                                                                                        agTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTgcgc               <  2:672876/150‑1 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGTGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATATCTTCAAGCGACATTTCAAGCTCTTCCAGATCTTCGTCTTCGTCGTAACCCAGTTTCGCaa             >  1:677447/1‑150 (MQ=255)

                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCaa             >  2:1249327/1‑150 (MQ=255)

                                                                                                                                           cGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAAt            <  2:3968704/150‑1 (MQ=255)

                                                                                                                                             gTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGt          >  1:1567562/1‑150 (MQ=255)

                                                                                                                                                tttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTAc       >  1:3371674/1‑150 (MQ=255)

                                                                                                                                                 ttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACg      <  1:601420/150‑1 (MQ=255)

TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037

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