BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,469,455	A→G	8.1%	T611A (ACC→GCC) ‡	mgtA →	magnesium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,469,455	0	A	G	8.1%	311.5 / 14.5	111	T611A (ACC→GCC) ‡	mgtA	magnesium transporter
Reads supporting (aligned to +/- strand): ref base A (50/52); new base G (4/5); total (54/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 112 total aligned reads displayed.

TGGAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACGGC  >  NC_000913/4469307‑4469604

       aCCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTg                                                                                                                                               >  1:725840/1‑150 (MQ=255)

        ccGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTTCCATGCAGTGGGGCTGGACGCGGGAGAGGTGGTCAATGGTAGTGATATAGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGt                                                                                                                                              >  1:1651622/1‑150 (MQ=255)

         cGTAAAAATCCTTCCTGGCGACCGTGAGTTATTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGTGAGATGGTGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtg                                                                                                                                             <  2:297144/150‑2 (MQ=255)

          gTAATAATCCTTACTGGCGACAGTGAGTTAGTTGCTGCGAATGTGTGCCATGAAGTGTGGCTGGACGCGGGTGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAATTCTTTCACAGCGTACCACGCTGttt                                                                                                                                            <  2:1651809/150‑1 (MQ=255)

          gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGttt                                                                                                                                            >  1:44624/1‑150 (MQ=255)

          gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGttt                                                                                                                                            >  1:1344634/1‑150 (MQ=255)

           tAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgca                                                                                                                                           <  1:372773/150‑4 (MQ=255)

           tAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgca                                                                                                                                           <  1:1085088/150‑4 (MQ=255)

           tAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTg                                                                                                                                           <  2:1495718/150‑1 (MQ=255)

            aaaaaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgcaa                                                                                                                                          <  1:1682472/150‑5 (MQ=255)

             aaaaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgcaag                                                                                                                                         <  2:1268116/150‑6 (MQ=255)

                 tCCTTACTGGCGACAGTGGGGTAGTTGCAGCGAAAGTGGGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGcc                                                                                                                                     <  1:617761/150‑1 (MQ=255)

                 tCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGcc                                                                                                                                     <  1:1689270/150‑1 (MQ=255)

                  ccTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCt                                                                                                                                    <  2:949234/150‑1 (MQ=255)

                   cTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTg                                                                                                                                   <  2:364119/150‑1 (MQ=255)

                   cTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTg                                                                                                                                   <  2:1584517/150‑1 (MQ=255)

                    ttACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGa                                                                                                                                  >  1:932397/1‑150 (MQ=255)

                     tACTGGCGACAGTGTGTTAGTTGCAGCGAAAGTGTGCCATGTAGTGGGGCTGGACGCGTTAGAGGTGGTCATTGGTAGTGATTTTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGAc                                                                                                                                 <  2:408678/150‑1 (MQ=255)

                      aCTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACg                                                                                                                                <  1:1534511/150‑1 (MQ=255)

                       cTGTCGACAGTGAGTTTGAGGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGTGAGTGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAAGCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGc                                                                                                                               <  1:257870/150‑1 (MQ=255)

                        tGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGcc                                                                                                                              <  2:1312206/150‑1 (MQ=255)

                        tGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGcc                                                                                                                              <  2:1003165/150‑1 (MQ=255)

                         ggCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGccgccg                                                                                                                             >  1:1159652/1‑150 (MQ=255)

                         ggCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCg                                                                                                                             >  1:1501935/1‑150 (MQ=255)

                            gACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACTAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCAGATg                                                                                                                          >  2:1220223/1‑150 (MQ=255)

                               aGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGCACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCAt                                                                                                                       <  2:919538/150‑1 (MQ=255)

                                gTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATa                                                                                                                      >  1:587006/1‑150 (MQ=255)

                                gTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATa                                                                                                                      >  2:1671540/1‑150 (MQ=255)

                                gTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATa                                                                                                                      >  1:1272320/1‑150 (MQ=255)

                                 tGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATaa                                                                                                                     <  2:1117582/150‑1 (MQ=255)

                                          tGCAGCGAACGTGTGACATGAAGGGCGGCTGGACGCGGGAGGGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCAt                                                                                                            <  1:1485449/150‑1 (MQ=255)

                                            cAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGCGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtg                                                                                                          <  2:379040/150‑1 (MQ=255)

                                             aGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtgt                                                                                                         >  1:1189987/1‑150 (MQ=255)

                                                gAAAGTGTGACATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCAATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGAc                                                                                                      <  2:328633/150‑1 (MQ=255)

                                                 aaaGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACTCCGATGCATAAAGAACGCATTGTGAcc                                                                                                     >  2:585724/1‑150 (MQ=255)

                                                 aaaGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGAcc                                                                                                     >  2:955849/1‑150 (MQ=255)

                                                   actgtgCCATGAATTGGGGCTGGACTCGGCAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCtt                                                                                                   <  2:846614/148‑1 (MQ=255)

                                                      gtgCCATGAAGTGGGGCTGGACGAGGGAGAGTTTGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACt                                                                                                <  2:1608475/150‑1 (MQ=255)

                                                        gCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACTCATTGTGACCTTACTGa                                                                                              >  2:1198441/1‑150 (MQ=255)

                                                          cATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAg                                                                                            <  1:968822/150‑1 (MQ=255)

                                                          cATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAg                                                                                            <  2:1195266/150‑1 (MQ=255)

                                                           aTGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGCCGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAgc                                                                                           >  2:851252/1‑150 (MQ=255)

                                                           aTGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAgc                                                                                           <  2:892309/150‑1 (MQ=255)

                                                             gAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGCCCTTACTGCAgcgc                                                                                         >  1:1694142/1‑150 (MQ=255)

                                                                 tGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAg                                                                                     <  1:240/150‑1 (MQ=255)

                                                                   gggCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAggg                                                                                   <  1:2460/150‑1 (MQ=255)

                                                                   gggCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAggg                                                                                   >  1:797663/1‑150 (MQ=255)

                                                                            gcgGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATATTGCACAGCGTACCACGCTGTTTGCCCGGCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGtt                                                                          <  2:1248365/150‑1 (MQ=255)

                                                                                  gagGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCCCGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCttt                                                                    <  2:1501935/150‑1 (MQ=255)

                                                                                  gagGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCttt                                                                    <  2:788358/150‑1 (MQ=255)

                                                                                  gagGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCttt                                                                    <  1:1181122/150‑1 (MQ=255)

                                                                                   agGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTa                                                                   >  2:980397/1‑150 (MQ=255)

                                                                                    ggtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTAt                                                                  >  1:103368/1‑150 (MQ=255)

                                                                                     gtggtCATTGGTAGTGATATTTAAACTCTATCTGACGAAGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATg                                                                 <  2:1629847/150‑1 (MQ=255)

                                                                                      tggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATgg                                                                >  2:1348838/1‑150 (MQ=255)

                                                                                       ggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATggg                                                               <  1:1618257/150‑1 (MQ=255)

                                                                                       ggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATggg                                                               <  2:1272320/150‑1 (MQ=255)

                                                                                       gggctatGGAAGTGAAATTGAAAAGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATggg                                                               <  1:1328933/144‑1 (MQ=255)

                                                                                        gtCATTGGTAGTTATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGc                                                              <  2:464261/150‑1 (MQ=255)

                                                                                        gtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGACCGCATTGTGACCTTACTGAAGCGCAAAGGGCATGTGGTTGGCTTTATGGGc                                                              >  1:1588556/1‑150 (MQ=255)

                                                                                        gtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCTTTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGc                                                              <  2:598269/150‑1 (MQ=255)

                                                                                        gtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGc                                                              >  1:1164970/1‑150 (MQ=255)

                                                                                         tCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCg                                                             >  2:774788/1‑150 (MQ=255)

                                                                                          cATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGgaga                                                            >  2:384237/1‑150 (MQ=255)

                                                                                           aTTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGAt                                                           >  2:1260539/1‑150 (MQ=255)

                                                                                            ttGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATg                                                          <  2:747420/150‑1 (MQ=255)

                                                                                             tGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATgg                                                         >  1:670925/1‑150 (MQ=255)

                                                                                                   tGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTaa                                                   >  1:229383/1‑150 (MQ=255)

                                                                                                    gATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAat                                                  >  2:1131599/1‑150 (MQ=255)

                                                                                                      tatTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATg                                                          <  2:1609402/140‑1 (MQ=255)

                                                                                                      tatTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatga                                                <  1:1726555/150‑1 (MQ=255)

                                                                                                       atTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatgat                                               >  1:557657/1‑150 (MQ=255)

                                                                                                       atTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatgat                                               >  2:49102/1‑150 (MQ=255)

                                                                                                        ttGAAACGCAATCTGACGACGAACTGGCAAATCTTGCACACCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatgatg                                              <  2:241866/150‑1 (MQ=255)

                                                                                                          gAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATgcg                                            <  2:805660/150‑1 (MQ=255)

                                                                                                              cGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGACGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcccg                                        <  1:391015/150‑1 (MQ=255)

                                                                                                              cGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcccg                                        <  2:1344634/150‑1 (MQ=255)

                                                                                                                cTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcccgcc                                      >  1:96117/1‑150 (MQ=255)

                                                                                                                 ggaatGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCt                                     >  1:1036958/5‑150 (MQ=255)

                                                                                                                     tGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTAcg                                 >  2:533577/1‑150 (MQ=255)

                                                                                                                        cgacgaACTGGCAAATCTAGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACgcgc                              >  1:1284191/1‑150 (MQ=255)

                                                                                                                        cgacgaACTCGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTGGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACgcgc                              <  1:996201/150‑1 (MQ=255)

                                                                                                                            gcagtGGCAAATCTTTCACAGCGGACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACTCATTTTTACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCg                          <  2:1557525/146‑1 (MQ=255)

                                                                                                                                tGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGata                      >  1:1652346/1‑150 (MQ=255)

                                                                                                                                 ggCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGatat                     >  1:1183817/1‑150 (MQ=255)

                                                                                                                                   cAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCg                   >  2:1577224/1‑150 (MQ=255)

                                                                                                                                      aTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCa                >  1:1730533/1‑150 (MQ=255)

                                                                                                                                      aTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCa                >  1:539084/1‑150 (MQ=255)

                                                                                                                                         gcaaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCAttt             >  2:564343/5‑150 (MQ=255)

                                                                                                                                          caaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCAtttg            >  1:755226/4‑149 (MQ=255)

                                                                                                                                          caaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTc            >  2:1375585/4‑150 (MQ=255)

                                                                                                                                           aaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCt           >  2:1410612/3‑150 (MQ=255)

                                                                                                                                            cacaGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCtg          >  1:811477/1‑150 (MQ=255)

                                                                                                                                            cacaGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTCACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCtg          <  2:887443/150‑1 (MQ=255)

                                                                                                                                             acaGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCtgt         >  1:1645300/1‑150 (MQ=255)

                                                                                                                                              caGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCtgtg        <  2:1354630/150‑1 (MQ=255)

                                                                                                                                              caGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACAGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATACGCCCGCCTTACGCGCTGCGGATATCGGCATTTCtgtg        >  1:1379866/1‑150 (MQ=255)

                                                                                                                                               aGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGACCGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTgg       <  1:853991/150‑1 (MQ=255)

                                                                                                                                               aGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTgg       <  1:1671540/150‑1 (MQ=255)

                                                                                                                                                    aCCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACGgc  <  1:533577/150‑1 (MQ=255)

TGGAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACGGC  >  NC_000913/4469307‑4469604

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