BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,367,702	A→C	5.0%	I209S (ATC→AGC) ‡	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,367,702	0	A	C	5.0%	391.1 / 6.4	130	I209S (ATC→AGC) ‡	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand): ref base A (64/58); new base C (3/4); total (67/63)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 132 total aligned reads displayed.

AGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATCAGG  >  NC_000913/4367553‑4367851

aGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                                                                                                                                       <  2:1643441/150‑1 (MQ=255)

        aCTTCAGCTAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                                                                                                               <  2:589797/150‑1 (MQ=255)

        aCTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                                                                                                               <  1:717524/150‑1 (MQ=255)

         cTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAg                                                                                                                                              >  1:99855/1‑150 (MQ=255)

         cTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAg                                                                                                                                              >  2:1059086/1‑150 (MQ=255)

         cTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAg                                                                                                                                              <  2:1385851/150‑1 (MQ=255)

         attCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAg                                                                                                                                              <  2:1892974/149‑1 (MQ=255)

          ttCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAga                                                                                                                                             <  2:1270744/150‑2 (MQ=255)

          ttCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAga                                                                                                                                             <  1:352160/150‑2 (MQ=255)

           tCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgcg                                                                                                                                            >  2:338250/1‑150 (MQ=255)

           tCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgcg                                                                                                                                            >  1:499402/1‑150 (MQ=255)

           tCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgcg                                                                                                                                            >  1:911728/1‑150 (MQ=255)

            cAGCCAGTTATTTCACTGCCAGCGGAGAGTACTCTTTCGGAGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAgaag                                                                                                                                           <  2:940733/150‑4 (MQ=255)

               ccAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGa                                                                                                                                        <  1:607124/150‑1 (MQ=255)

                    tttttCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTcc                                                                                                                                   <  1:485170/150‑1 (MQ=255)

                    tttttCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTcc                                                                                                                                   <  2:526292/150‑1 (MQ=255)

                         cACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGAcc                                                                                                                              <  1:1024203/150‑1 (MQ=255)

                           cTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGa                                                                                                                            >  1:1056166/1‑150 (MQ=255)

                           cTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGa                                                                                                                            >  1:1064489/1‑150 (MQ=255)

                           cTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGa                                                                                                                            >  2:1116987/1‑150 (MQ=255)

                             gCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAgg                                                                                                                          >  1:1346581/1‑150 (MQ=255)

                                     agagTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCgg                                                                                                                  <  1:1545564/150‑1 (MQ=255)

                                      gagTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGt                                                                                                                 <  1:1907868/150‑1 (MQ=255)

                                        gTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTAc                                                                                                               <  2:1319710/150‑1 (MQ=255)

                                         tACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACt                                                                                                              <  2:1029792/150‑1 (MQ=255)

                                            cctTTCGGCGTGTTCAGACCAGTACCGATTGCTGGTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGCACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCg                                                                                                           <  1:96621/149‑1 (MQ=255)

                                             ctTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGt                                                                                                          >  2:449723/1‑150 (MQ=255)

                                               ttCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCAACAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTcc                                                                                                        <  2:333868/150‑1 (MQ=255)

                                                 cGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTg                                                                                                      >  2:1856828/1‑150 (MQ=255)

                                                       gtTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTg                                                                                                >  1:354163/1‑150 (MQ=255)

                                                       gtTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTg                                                                                                <  1:876069/150‑1 (MQ=255)

                                                         tCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTggg                                                                                              <  2:892643/150‑1 (MQ=255)

                                                         tCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTggg                                                                                              >  1:452077/1‑150 (MQ=255)

                                                             aCCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACg                                                                                          >  1:37537/1‑150 (MQ=255)

                                                              ccAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACgc                                                                                         >  2:216905/1‑149 (MQ=255)

                                                              ccAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGa                                                                                         >  2:645022/1‑150 (MQ=255)

                                                              ccAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGGTCCGCGGTCCGTTGGATGTGTTTCACTTCTTCTTTCAGCAGGCTGCGGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGAGTCCTGCAGCTGGTTACgg                                                                                         >  1:239110/1‑149 (MQ=255)

                                                                   aCCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCAt                                                                                    >  1:897190/1‑150 (MQ=255)

                                                                      gATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCAtttt                                                                                 >  1:1637642/1‑150 (MQ=255)

                                                                        ttGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTATCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCa                                                                               <  1:899311/150‑1 (MQ=255)

                                                                         tGCTGTTGGACCTATGTTAACTTCCTGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGGCATCGGTACTGCGTCCTGCAGCTGGGTACGAGCCATTTTCAg                                                                              <  1:243391/150‑1 (MQ=255)

                                                                         tGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCTGAATTCCTGACCCAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAg                                                                              >  2:1791119/1‑150 (MQ=255)

                                                                         tGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAg                                                                              >  1:31380/1‑150 (MQ=255)

                                                                         tGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAg                                                                              >  2:1159066/1‑150 (MQ=255)

                                                                         tGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAg                                                                              <  2:1695290/150‑1 (MQ=255)

                                                                         tGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGGTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAg                                                                              <  1:726416/150‑1 (MQ=255)

                                                                         tGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTGCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGGCATCGGTACTGCGTCCTGCAGTTGGGTACGACCCATTTTCAg                                                                              <  1:226477/150‑1 (MQ=255)

                                                                          gCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAgg                                                                             <  2:1444674/150‑1 (MQ=255)

                                                                          gCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAgg                                                                             >  2:1552872/1‑150 (MQ=255)

                                                                           cTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGa                                                                            >  1:442140/1‑150 (MQ=255)

                                                                           cTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGAGGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGa                                                                            <  1:793306/150‑1 (MQ=255)

                                                                                gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTcc                                                                       >  1:1819754/1‑150 (MQ=255)

                                                                                gCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTcc                                                                       >  1:1525824/1‑150 (MQ=255)

                                                                                 cACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCGTGCAGCTGGGTACGACCCATTTTCAGGATGTCCt                                                                      >  1:308386/1‑150 (MQ=255)

                                                                                 cACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGCCCCATTTTCAGGATGTCCt                                                                      >  2:971065/1‑150 (MQ=255)

                                                                                   ccAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTgg                                                                    <  1:1227175/150‑1 (MQ=255)

                                                                                   ccAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTgg                                                                    >  1:724313/1‑150 (MQ=255)

                                                                                     aaGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGGCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGaa                                                                  <  2:1389462/150‑1 (MQ=255)

                                                                                      aGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTTCGTCCTTCCGCTGGGTACGACCCATTTTCAGGATGTCCTGGAAt                                                                 >  2:1684943/1‑150 (MQ=255)

                                                                                      aGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAAt                                                                 >  2:1684942/1‑150 (MQ=255)

                                                                                       ggTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTTCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAttt                                                                >  1:282420/1‑150 (MQ=255)

                                                                                        gTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTc                                                               <  1:526543/150‑1 (MQ=255)

                                                                                         ttAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCg                                                              <  1:1856828/150‑1 (MQ=255)

                                                                                          tAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATCTTCAGGATGTCCTGGAATTCGa                                                             >  2:1346811/1‑150 (MQ=255)

                                                                                            aCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACa                                                           <  2:1755279/150‑1 (MQ=255)

                                                                                             nnnccAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTAACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGATCATCTGTACTGCGTACTGCAGCTGGGTACGACCCATTTTCAGGATGTCATGGAATTCGACAg                                                          >  2:98281/4‑150 (MQ=255)

                                                                                             cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAg                                                          >  2:782667/1‑150 (MQ=255)

                                                                                             cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAg                                                          <  1:140087/150‑1 (MQ=255)

                                                                                              ttCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGc                                                         <  1:980224/150‑1 (MQ=255)

                                                                                              ttCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGc                                                         <  1:1712506/150‑1 (MQ=255)

                                                                                               tCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACCGCt                                                        >  2:1331119/1‑150 (MQ=255)

                                                                                                ccAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCtt                                                       <  1:833492/150‑1 (MQ=255)

                                                                                                 cagcagcagTTCAGCGGTACGTTGGATGTTGTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCttt                                                      >  1:753451/1‑150 (MQ=255)

                                                                                                   gcagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACTCATTTTCAGGATGTCCTGGAATTCGACAGCTTTAc                                                    >  2:833670/1‑150 (MQ=255)

                                                                                                     agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGt                                                  <  1:1059086/150‑1 (MQ=255)

                                                                                                     agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGt                                                  <  1:1249763/150‑1 (MQ=255)

                                                                                                       cagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTc                                                <  1:868966/150‑1 (MQ=255)

                                                                                                         gTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCaa                                              >  2:643014/1‑150 (MQ=255)

                                                                                                         gTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCaa                                              <  1:783305/150‑1 (MQ=255)

                                                                                                           tCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAg                                            >  2:1838355/1‑150 (MQ=255)

                                                                                                            cAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCAATTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGc                                           >  2:914469/1‑150 (MQ=255)

                                                                                                                ggTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTc                                       <  2:1346581/150‑1 (MQ=255)

                                                                                                                        ggATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGtt                               >  2:76722/1‑150 (MQ=255)

                                                                                                                         gATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTg                              <  1:1730506/150‑1 (MQ=255)

                                                                                                                         gATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACAGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGAACCATTTTCAGGATGTCCTGGCATTCGACAGCTTTACGTTCAAAGACTTCACGCAGTTg                              >  1:126378/1‑150 (MQ=255)

                                                                                                                           tGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGt                            <  1:1657357/150‑1 (MQ=255)

                                                                                                                             tttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGATGGGTACGACCCATTTTCAGGATGTCCTGGTATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTa                          >  1:1141374/1‑150 (MQ=255)

                                                                                                                                 tCACGTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTACTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCg                      <  2:1258524/150‑1 (MQ=255)

                                                                                                                                  cACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGc                     >  1:1716870/1‑150 (MQ=255)

                                                                                                                                   aCTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACAGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGCCCCATTTTCAGGATGTCCTGGAATTCGACAGCTATACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCa                    >  1:1225512/1‑150 (MQ=255)

                                                                                                                                     ttcttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGCCCGAGGGTAATCGGTACTGCGTCCTGCAGCTGGGTACGCCCCATTTTCAGGATGTCCTGGAATTCGCCAGCTGTACGATCATAGCCTTCACGCAGTTGGTTAATCGCATc                  >  1:331037/1‑150 (MQ=255)

                                                                                                                                      tcttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCt                 >  2:68472/1‑150 (MQ=255)

                                                                                                                                         tcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAcc              <  2:499402/150‑1 (MQ=255)

                                                                                                                                         tcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAcc              >  1:165970/1‑150 (MQ=255)

                                                                                                                                          cttTCAGCATCCTGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAccc             >  2:161323/1‑149 (MQ=255)

                                                                                                                                          cttTCAGCATCCTGCTGAAAGCGCGGAATACCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATAGCATCTACCa             >  2:894642/1‑150 (MQ=255)

                                                                                                                                              cagcagGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACTCAGTTGGTTAATCGCATCTACCAGCtt         >  2:1761388/1‑150 (MQ=255)

                                                                                                                                                   ggATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGTAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATCa    >  2:592694/1‑150 (MQ=255)

                                                                                                                                                    gATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATCAg   >  2:559407/1‑150 (MQ=255)

                                                                                                                                                     aTGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATCAgg  <  2:387323/150‑1 (MQ=255)

AGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCAGCTTAATCAGG  >  NC_000913/4367553‑4367851

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