BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,174,504	A→T	6.7%	L221Q (CTG→CAG)	ycfT ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,174,504	0	A	T	6.7%	288.5 / 6.7	99	L221Q (CTG→CAG)	ycfT	inner membrane protein
Reads supporting (aligned to +/- strand): ref base A (53/39); new base T (2/5); total (55/44)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

TCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTA  >  NC_000913/1174355‑1174651

tCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATcagca                                                                                                                                                     >  2:208967/1‑150 (MQ=255)

   gcagGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACACCCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGAt                                                                                                                                                  >  2:910897/1‑150 (MQ=255)

        cTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATggcgg                                                                                                                                             <  1:1009410/150‑1 (MQ=255)

           gtggCGCGCATACCGAATCGTTGCTCGTATTGCTAAAACAGCTTCATGATCACCACAATCGCTACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGCCCGGTAAGCAAAGAAGCCATCAGCAGATGGCGgcgc                                                                                                                                          >  2:1003125/1‑150 (MQ=255)

           gtggAGCGCATACCGAATCGTTGCTCGTATTGATAAAAGAGCTTCATGATCACCACAATCGATACGAGCGATAACAGCAGCGAGATAGTAAACAACCAGGCACAGACCGCCAGAACGGTAAGCAAAGAAGCCATCTGCTGATggcttctt                                                                                                                                          <  2:201958/150‑6 (MQ=255)

           gtggAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCTGCTGATggcttctt                                                                                                                                          <  1:1693902/150‑6 (MQ=255)

           gtggAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGCAAGCAAAGAAGCCATCTGCTGATggcttctt                                                                                                                                          <  2:515250/150‑6 (MQ=255)

            tggAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAATGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCa                                                                                                                                         <  1:1352002/150‑1 (MQ=255)

            tggAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCTGCTGATggcttcttt                                                                                                                                         <  2:1733777/150‑7 (MQ=255)

             ggAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCTGCTGATggcttctttg                                                                                                                                        <  1:562155/150‑8 (MQ=255)

               aGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACg                                                                                                                                      <  1:1732805/150‑1 (MQ=255)

                 cgcATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACgga                                                                                                                                    <  1:1027504/150‑2 (MQ=255)

                 cgcATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGAACCGACCGCCAGAACGGTAAGCAAAGAAGCCAACAGCAGATGGCGGCGCAACGGc                                                                                                                                    <  1:302128/150‑1 (MQ=255)

                 cgcATACCGAATCGTTGCGCGTATTGATAAAACAGCTTCATGATCACCACAATCGAGACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGc                                                                                                                                    <  1:959995/150‑1 (MQ=255)

                   cATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGACATCAGCAGATGGCGGCGCAACGGCAc                                                                                                                                  >  1:438340/1‑150 (MQ=255)

                     tACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGAGACAACAGCAGCGAGATAGTAAACAACCAGGCCCCGACCGCCAGAGCGGTAAGCCAAGACGCCATCAGCAGATGGCGGCGCAACGGCACct                                                                                                                                >  2:273217/1‑150 (MQ=255)

                       ccGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACctct                                                                                                                              >  2:1535598/1‑150 (MQ=255)

                         gAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCttt                                                                                                                            <  2:607038/150‑1 (MQ=255)

                            tCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAac                                                                                                                         <  1:709830/150‑1 (MQ=255)

                            tCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAac                                                                                                                         <  2:467051/150‑1 (MQ=255)

                              gTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAacac                                                                                                                       >  1:533054/1‑150 (MQ=255)

                                 gCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAgg                                                                                                                    <  1:1127487/150‑1 (MQ=255)

                                   tCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTc                                                                                                                  <  2:1402564/150‑1 (MQ=255)

                                    cGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCTATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCa                                                                                                                 >  1:438711/1‑150 (MQ=255)

                                    cGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCATCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCa                                                                                                                 >  1:340504/1‑150 (MQ=255)

                                         tGATAAAACAGCTTCATGATCACCACAATAGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCattatt                                                                                                            >  2:946739/1‑150 (MQ=255)

                                          gATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACACCCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTAttg                                                                                                           <  1:675826/150‑1 (MQ=255)

                                            tAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTAttgtt                                                                                                         <  2:448422/150‑1 (MQ=255)

                                             aaaaCAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCCGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACAGTAAGCAAAGAAGCCATCAGCAGATGGCGCCGCAACGGCACCTCTATAACCCAGGTCATTAttgtac                                                                                                        >  1:80340/1‑148 (MQ=255)

                                              aaaCAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTgc                                                                                                       >  2:653980/1‑150 (MQ=255)

                                                aCAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTgcgc                                                                                                     <  1:402931/150‑1 (MQ=255)

                                                     ttCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAAc                                                                                                >  1:1430189/1‑150 (MQ=255)

                                                     ttCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAAc                                                                                                >  1:452876/1‑150 (MQ=255)

                                                      tCATGATCACCACAATCGATACCAGCGACAACCGCAGCGAGATAGTACACAACCAGGCACCGAACGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAAcc                                                                                               <  1:1448985/150‑1 (MQ=255)

                                                      tCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAAcc                                                                                               <  1:1126317/150‑1 (MQ=255)

                                                        aTGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAt                                                                                             >  1:1269636/1‑150 (MQ=255)

                                                         tGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCTCCACGGCACCTCTTTAACACAGGTCATTATTTTTGCGCCAAACCATg                                                                                            >  2:1064103/1‑150 (MQ=255)

                                                           aTCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATgcg                                                                                          <  2:1269636/150‑1 (MQ=255)

                                                           aTCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATgcg                                                                                          <  1:1020111/150‑1 (MQ=255)

                                                               ccacAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCACACCATGCGCCaa                                                                                      >  1:466752/1‑150 (MQ=255)

                                                                cacaATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAg                                                                                     <  1:208967/150‑1 (MQ=255)

                                                                 acaATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCCCCGACCGCCAGAACGGTAAGCAAAGCAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAgg                                                                                    >  2:1512298/1‑150 (MQ=255)

                                                                 acaATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAgg                                                                                    >  2:1706065/1‑150 (MQ=255)

                                                                  caATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGc                                                                                   <  1:461591/150‑1 (MQ=255)

                                                                  caATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGc                                                                                   >  2:748481/1‑150 (MQ=255)

                                                                     acGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGt                                                                                <  1:1579353/149‑1 (MQ=255)

                                                                           ccAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCATCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAAGCGCCAAGGCTGTAATAAg                                                                          >  2:453085/1‑150 (MQ=255)

                                                                           ccAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAg                                                                          >  1:983311/1‑150 (MQ=255)

                                                                             aGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGc                                                                        >  2:605330/1‑150 (MQ=255)

                                                                                gACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCATCGGCACCTCTTTAACACAGGTCATTATTGTTACGCCAAACCATGCGCCAAGGCTGTAATAAGGCaaa                                                                     >  2:733353/1‑150 (MQ=255)

                                                                                   aacaGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTg                                                                  >  2:1445700/1‑150 (MQ=255)

                                                                                   aacaGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTg                                                                  >  1:478054/1‑150 (MQ=255)

                                                                                    acaGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCTTTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGc                                                                 >  1:176911/1‑150 (MQ=255)

                                                                                    acaGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGc                                                                 >  2:665783/1‑150 (MQ=255)

                                                                                       gcagcGAGATAGTAAACAACCAGGCACCGACCGCCAGATCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGa                                                              <  1:287507/150‑1 (MQ=255)

                                                                                       gcagcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGa                                                              <  2:1162799/150‑1 (MQ=255)

                                                                                          gcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATca                                                           >  1:517284/1‑150 (MQ=255)

                                                                                            gagaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCTATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATcaca                                                         <  2:438711/150‑1 (MQ=255)

                                                                                            gagaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATcaca                                                         <  2:340504/150‑1 (MQ=255)

                                                                                             agaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCacac                                                        >  2:1115316/1‑150 (MQ=255)

                                                                                              gaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACt                                                       <  1:1119283/150‑1 (MQ=255)

                                                                                               aTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTg                                                      >  2:939398/1‑150 (MQ=255)

                                                                                                 aGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGtt                                                    >  1:1182286/1‑150 (MQ=255)

                                                                                                     aacaacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCAtt                                                >  1:1681360/1‑150 (MQ=255)

                                                                                                       caacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGGCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTcc                                              <  1:1216109/150‑1 (MQ=255)

                                                                                                        aacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAACCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTccc                                             >  1:1061904/1‑150 (MQ=255)

                                                                                                        aacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTccc                                             >  1:1056560/1‑150 (MQ=255)

                                                                                                           cAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAACCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccac                                          >  1:387011/1‑150 (MQ=255)

                                                                                                           cAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccac                                          >  1:1566791/1‑150 (MQ=255)

                                                                                                           cAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccac                                          >  2:620602/1‑150 (MQ=255)

                                                                                                            aGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacc                                         <  2:631097/150‑1 (MQ=255)

                                                                                                                  cgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCg                                   <  1:1344657/150‑1 (MQ=255)

                                                                                                                   gaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGt                                  >  2:1045732/1‑150 (MQ=255)

                                                                                                                     ccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTgg                                <  2:533054/150‑1 (MQ=255)

                                                                                                                      cgCCAGCTCGGTATGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCCCGGCGTggg                               >  2:616476/1‑150 (MQ=255)

                                                                                                                      cgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGTGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTggg                               <  2:1599768/150‑1 (MQ=255)

                                                                                                                      cgCCAGAACGGTAAGAAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACTCAGGTCAGTATTGTTGCGCAAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCCGATCACACTGATCATTCCCCCCCACGGCGTggg                               >  2:847639/1‑150 (MQ=255)

                                                                                                                         cAGAACGGTAAGCAATGAAGCCATCAGCAGATGGCGGCGCAACGGCACCGCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAAc                            <  1:1165213/150‑1 (MQ=255)

                                                                                                                          aGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACg                           >  1:600278/1‑150 (MQ=255)

                                                                                                                          aGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGAGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACg                           <  2:1625337/150‑1 (MQ=255)

                                                                                                                              cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAt                       >  2:1568286/1‑150 (MQ=255)

                                                                                                                              cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAt                       <  1:946739/150‑1 (MQ=255)

                                                                                                                                  aaGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAACCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAAc                   >  1:925484/1‑150 (MQ=255)

                                                                                                                                  aaGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGCTCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAAc                   >  2:450722/1‑150 (MQ=255)

                                                                                                                                  aaGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAAc                   >  1:1175699/1‑150 (MQ=255)

                                                                                                                                  aaGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAAc                   >  1:1628177/1‑150 (MQ=255)

                                                                                                                                  aaGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACAATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGTAACGAAATTaaa                   >  2:1399817/1‑149 (MQ=255)

                                                                                                                                         gccgCCATCTGCTGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCcaca            >  2:382263/4‑150 (MQ=255)

                                                                                                                                         gccgCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCcaca            >  1:1388372/4‑150 (MQ=255)

                                                                                                                                          ccgCCATCTGCTGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCAGTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCacac           >  1:1113973/3‑150 (MQ=255)

                                                                                                                                          aagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCacac           <  2:735308/150‑1 (MQ=255)

                                                                                                                                             ccATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTca        >  1:1006919/1‑150 (MQ=255)

                                                                                                                                             ccATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTca        >  1:791301/1‑150 (MQ=255)

                                                                                                                                             ccATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTca        <  1:748481/150‑1 (MQ=255)

                                                                                                                                              cATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCcagcaccag       >  2:248940/1‑143 (MQ=255)

                                                                                                                                              cATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcag       >  1:589880/1‑150 (MQ=255)

                                                                                                                                               aTCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcagc      <  1:939398/150‑1 (MQ=255)

                                                                                                                                                 cagcagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCCCCACGGCGTGGGCACGAAATTAACCGCCACACTcagcag    >  1:720126/1‑150 (MQ=255)

                                                                                                                                                 cagcagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcagcag    >  1:1651341/1‑150 (MQ=255)

                                                                                                                                                   gcagATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTa  >  1:1295449/1‑150 (MQ=255)

TCAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTA  >  NC_000913/1174355‑1174651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`