BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,469,455	A→G	9.7%	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	9.7%	299.4 / 20.8	113	T611A (ACC→GCC) ‡	mgtA	magnesium transporter
Reads supporting (aligned to +/- strand): ref base A (55/47); new base G (4/7); total (59/54)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

Only 100 of 113 total aligned reads displayed.

GAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACG  >  NC_000913/4469309‑4469602

gAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCa                                                                                                                                                  <  2:200571/150‑1 (MQ=255)

 aaTTACCATAAAAATCCTTACTGGCGAAAGTGAGTTAGTTGCAGCGAAAGTCTGACATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCAc                                                                                                                                                 >  1:574681/1‑150 (MQ=255)

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

       cGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGAAGCGCAAGTGTGCCATGAAGTGGGGCTGGACGAGGGAGAGGTGGTCATTGGTAGTGAGATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACGACGCTGtt                                                                                                                                           >  1:1230843/1‑150 (MQ=255)

        gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgc                                                                                                                                          <  1:470295/150‑3 (MQ=255)

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

          aaaaaTCCTTACTTGCGACCTTGAGTTAGTTGCAGCTAAAGTGTGCCATGAAGGGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGGGATATTGAAACGCTATCTGACGACTAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgcaa                                                                                                                                        <  1:1073623/150‑5 (MQ=255)

          aaaaaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgcaa                                                                                                                                        <  2:128360/150‑5 (MQ=255)

          aaaaaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGc                                                                                                                                        <  1:999765/150‑1 (MQ=255)

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

           aaaaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGcc                                                                                                                                       <  1:428264/150‑1 (MQ=255)

              aTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGc                                                                                                                                    >  1:1442232/1‑150 (MQ=255)

                    aCTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACAACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACg                                                                                                                              >  1:353476/1‑150 (MQ=255)

                        gCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTACCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGa                                                                                                                          >  2:127625/1‑150 (MQ=255)

                           aCAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGc                                                                                                                       >  1:30113/1‑150 (MQ=255)

                              gTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGCCGCCGATGCaga                                                                                                                    >  1:834554/1‑148 (MQ=255)

                               tGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGCCTGGACGCGGGAGAGGTGGTCATTGGTAGAGATATAGAAACGCTATCTGACGACGAACTGGCAAATCTTCCACAGCGAACCACGCTGTTTGCCCTCCTGACGCCGATGCATaa                                                                                                                   >  2:1120971/1‑150 (MQ=255)

                                gAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATaaa                                                                                                                  >  1:225152/1‑150 (MQ=255)

                                 agttagttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAg                                                                                                                 <  2:1442232/150‑1 (MQ=255)

                                     agttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACg                                                                                                             <  1:1244092/150‑1 (MQ=255)

                                     agttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACg                                                                                                             >  1:612056/1‑150 (MQ=255)

                                     agttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACg                                                                                                             <  1:573625/150‑1 (MQ=255)

                                      gttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGc                                                                                                            >  2:1302176/1‑150 (MQ=255)

                                      gttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGc                                                                                                            >  1:1278360/1‑150 (MQ=255)

                                      gttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGc                                                                                                            >  1:553961/1‑150 (MQ=255)

                                        ggCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCAt                                                                                                          >  1:434333/2‑150 (MQ=255)

                                         gCAGCGAAAGTGTGCCATTAAGGGGGGCGGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCAtt                                                                                                         <  2:699219/150‑1 (MQ=255)

                                          cAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGATGTGGTCATTGGTAGAGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtg                                                                                                        >  2:30628/1‑150 (MQ=255)

                                          cAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtg                                                                                                        >  2:1452411/1‑150 (MQ=255)

                                           aGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtgt                                                                                                       <  2:823181/150‑1 (MQ=255)

                                              gAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGAc                                                                                                    >  1:1473323/1‑150 (MQ=255)

                                              gAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGAc                                                                                                    >  1:1256567/1‑150 (MQ=255)

                                               aaaGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGAcc                                                                                                   <  2:225152/150‑1 (MQ=255)

                                                aaGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCt                                                                                                  <  1:1420395/150‑1 (MQ=255)

                                                  gtgtGCCATGAAGGGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTa                                                                                                <  1:619038/150‑1 (MQ=255)

                                                   tgtgCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTAc                                                                                               >  2:514418/1‑150 (MQ=255)

                                                   tgtgCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTAc                                                                                               <  2:553961/150‑1 (MQ=255)

                                                   tgtgCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTAc                                                                                               <  2:615526/150‑1 (MQ=255)

                                                    gtgCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACt                                                                                              >  2:358575/1‑150 (MQ=255)

                                                      gCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACtaa                                                                                            <  1:421765/150‑3 (MQ=255)

                                                      gCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGa                                                                                            >  2:852894/1‑150 (MQ=255)

                                                      gCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGa                                                                                            >  2:400558/1‑150 (MQ=255)

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

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

                                                            aaGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGcgcg                                                                                      >  1:698101/1‑150 (MQ=255)

                                                             aGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGa                                                                                     >  1:687519/1‑150 (MQ=255)

                                                                    cTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGATGGGCAt                                                                              >  1:171446/1‑150 (MQ=255)

                                                                     gggACGCGGGAGAGGTGGTCATTGGTAGTGATTTTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCGGTTTGCCCGCCTGACGCCGATGCTTAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCAtg                                                                             <  1:1120971/149‑1 (MQ=255)

                                                                        ccgcgGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTgg                                                                          <  2:782066/149‑1 (MQ=255)

                                                                         cgcgGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGt                                                                         <  1:108757/150‑1 (MQ=255)

                                                                            gggAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTgg                                                                      >  1:1512278/1‑150 (MQ=255)

                                                                               agagGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACCAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCtt                                                                   >  1:153850/1‑150 (MQ=255)

                                                                               agagGTGGTCATTGGTAGTGATATTGAAACGCTATCTCACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCtt                                                                   <  2:1163313/150‑1 (MQ=255)

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

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

                                                                                   gtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATg                                                               >  2:1070484/1‑150 (MQ=255)

                                                                                     ggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATggg                                                             >  2:1477173/1‑150 (MQ=255)

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

                                                                                      gtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGTGCGAAGGGCATGTGGTTGGCTTTATGGGc                                                            >  2:30831/1‑150 (MQ=255)

                                                                                      gtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGc                                                            <  1:37836/150‑1 (MQ=255)

                                                                                        cATTGGTAGTGATATTGAAACGCTATCTGAGGACTAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAATGAAAGCATTGTGACCTTACTGAAGCGAGTAGGGCATGTGGTTGGCTTTATGGGCGa                                                          <  1:500029/150‑1 (MQ=255)

                                                                                          ttGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATg                                                        >  1:1037583/1‑150 (MQ=255)

                                                                                                gTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTa                                                  >  2:1458868/1‑150 (MQ=255)

                                                                                                    tatTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatga                                              >  2:572525/1‑150 (MQ=255)

                                                                                                      ttGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatgatg                                            <  1:1392910/150‑1 (MQ=255)

                                                                                                      ttGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatgatg                                            <  1:520756/150‑1 (MQ=255)

                                                                                                        gAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCTTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGCACGCATTGAGACCTTACTGACGCGCGAAGGGCATGTGGGTCGCTTTATGGTCGATGCTATTAATGATgcg                                          >  1:314384/1‑150 (MQ=255)

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

                                                                                                          aaCGCTATCTGACGACGAACTGGCAAATCTTGCATAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGACCGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcc                                        >  2:262966/1‑150 (MQ=255)

                                                                                                          aaCGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcc                                        <  1:400558/150‑1 (MQ=255)

                                                                                                          aaCGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcc                                        <  1:1452411/150‑1 (MQ=255)

                                                                                                          aaCGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcc                                        >  1:665389/1‑150 (MQ=255)

                                                                                                          aaCGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcc                                        >  2:744591/1‑150 (MQ=255)

                                                                                                               tATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCt                                   >  1:1293787/1‑150 (MQ=255)

                                                                                                                    gacgacGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTAcgc                              >  1:438539/1‑150 (MQ=255)

                                                                                                                       gacgaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCt                           >  2:833617/1‑150 (MQ=255)

                                                                                                                         cgaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGc                         <  2:153850/150‑1 (MQ=255)

                                                                                                                           aaCTGTCAAATCTTGCACAGCGTACGACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCgg                       <  2:817953/150‑1 (MQ=255)

                                                                                                                           aaCTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTTTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCgg                       <  2:1278360/150‑1 (MQ=255)

                                                                                                                           aaCTGGCAAATCTTGCACAGCGTACCACGCCGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCgg                       <  2:416758/150‑1 (MQ=255)

                                                                                                                            aCTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGa                      >  1:372356/1‑150 (MQ=255)

                                                                                                                            aCTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGa                      <  1:1302176/150‑1 (MQ=255)

                                                                                                                              tGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGata                    <  2:470539/150‑1 (MQ=255)

                                                                                                                                  aaaTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCgg                <  1:852894/150‑1 (MQ=255)

                                                                                                                                    cgggcaaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCa              >  1:1472741/8‑150 (MQ=255)

                                                                                                                                    aTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCa              <  2:438539/150‑1 (MQ=255)

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

                                                                                                                                    aTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGGTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCa              >  2:644363/1‑150 (MQ=255)

                                                                                                                                      ggcaaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCAtt            >  2:1548564/6‑150 (MQ=255)

                                                                                                                                       ttGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTAGGCGCTGCGGATATCGGCAttt           >  1:1044080/1‑150 (MQ=255)

                                                                                                                                       gcaaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGTGCCCGCCTTACGCGCTGCGGAGATCGGCAttt           >  1:309413/5‑150 (MQ=255)

                                                                                                                                        caaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTc          >  1:439094/4‑150 (MQ=255)

                                                                                                                                            caGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAATGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGTCATTTCtgtg      >  1:111788/1‑150 (MQ=255)

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

                                                                                                                                             aGCGTACCACGCTGTTTGTCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTgg     <  1:992629/150‑1 (MQ=255)

                                                                                                                                               cGTACCACGCTGTTTGCCCGCCTTACTCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGAc   <  2:574681/150‑1 (MQ=255)

                                                                                                                                               cGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGAc   >  1:764959/1‑150 (MQ=255)

                                                                                                                                               cGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGAc   <  2:1655892/150‑1 (MQ=255)

                                                                                                                                                gTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGTGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACg  <  1:30831/150‑1 (MQ=255)

                                                                                                                                                gTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACg  >  2:623755/1‑150 (MQ=255)

GAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCTGTGGACG  >  NC_000913/4469309‑4469602

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