BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,891,504	T→A	6.2%	intergenic (+21/‑111)	tnaB → / → mdtL	tryptophan transporter of low affinity/multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,891,504	0	T	A	6.2%	771.2 / 26.8	216	intergenic (+21/‑111)	tnaB/mdtL	tryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand): ref base T (111/90); new base A (8/6); total (120/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 220 total aligned reads displayed.

TCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCG  >  NC_000913/3891249‑3891776

        tCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCCGTGCTTGCAATATAAGCTCGCAAGAAGTTTCCCACTCAGACGTTCACGGTCTCGGGCGGCAATCTTATTCCGGCTACTGTCATTCTCTTTGGTATAACCGTCATTTTGTGCTGGTTCGGCACCGTCTTTCAAGTGTTACTCACTTTTGTCTTAATCTTTCTCGAAGCCACCCTATCCCTGActtctttaccctcccgtaa                                                                                                                                                                                                                                                        >  4:132922/1‑261 (MQ=255)

                            gggATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGACTGGGGCGGCAATCTTATTCCGGCGCTTGTCATTCTCTTTGGTATAACCGTGATTTTTTGCTGGTTCGGCCACGTCTTTAACGTTTTCCCTAAATTTGGCTTAATCCTTCAAGAAGCCACCCTTTCTCTGGCTTCTTGCCTCTCAGGAAATCCCTTATGTCGAAATGGCa                                                                                                                                                                                                                                    >  4:250282/1‑274 (MQ=255)

                                                 tgCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTACATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTTCCTCTCAGGTATTCACTTATGTCCAAATGGCAACTAGCCCGATCCTCCTTcacc                                                                                                                                                                                                              >  2:173132/1‑275 (MQ=255)

                                                                       ttattCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtc                                                                                                                                                                                                                        <  4:93784/243‑1 (MQ=255)

                                                                              cGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTactg                                                                                                                                                                                  >  2:150060/1‑271 (MQ=255)

                                                                                                           agTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCAt                                                                                                                                                                                                                                                                                  <  8:760642/149‑1 (MQ=255)

                                                                                                             tttCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATCGTCATTCTCTTTGGTATAACCGTGACTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTc                                                                                                                                                                                                                                                                                <  5:1187296/149‑1 (MQ=255)

                                                                                                               tCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGAAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGc                                                                                                                                                                                                                                                                              >  8:1253752/1‑149 (MQ=255)

                                                                                                                 ccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGAATGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                                                                                                                                                                                                            >  7:831196/1‑149 (MQ=255)

                                                                                                                  cAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                                                                                                                                                                                                           <  5:862067/149‑1 (MQ=255)

                                                                                                                    aTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGct                                                                                                                                                                                                                                                                         >  7:1428015/1‑149 (MQ=255)

                                                                                                                     tCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGctt                                                                                                                                                                                                                                                                        <  5:1433792/149‑1 (MQ=255)

                                                                                                                      cAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATGCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttc                                                                                                                                                                                                                                                                       <  5:1124107/149‑1 (MQ=255)

                                                                                                                        gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttctt                                                                                                                                                                                                                                                                     >  6:849018/1‑149 (MQ=255)

                                                                                                                         aTGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                                                                                                                                                                                                    >  6:564979/1‑149 (MQ=255)

                                                                                                                            ttCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                                                                                                                                                                                                   >  7:785999/1‑147 (MQ=255)

                                                                                                                              cACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctct                                                                                                                                                                                                                                                               >  7:1013089/1‑149 (MQ=255)

                                                                                                                                   tCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTAc                                                                                                                            >  2:108005/1‑275 (MQ=255)

                                                                                                                                     tGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaaa                                                                                                                                                                                                                                                        <  6:1427512/149‑1 (MQ=255)

                                                                                                                                       gggCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                                                                                                                                                                                                      >  5:1341039/1‑149 (MQ=255)

                                                                                                                                       gggCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                                                                                                                                                                                                      <  6:493439/149‑1 (MQ=255)

                                                                                                                                           ggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCACCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACtt                                                                                                                                                                                                                                                  <  8:1379535/149‑1 (MQ=255)

                                                                                                                                             cAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCCACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCCGGACGCTTTAGCCTATCTCCAGCTTTTAGACTTGATGTTTTGCCCTGGTTTTACttttctccgc                                                                                                                  >  4:135885/1‑268 (MQ=255)

                                                                                                                                              aaTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                                                                                                                                                                                               >  8:684673/1‑149 (MQ=255)

                                                                                                                                              aaTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                                                                                                                                                                                               <  7:230385/149‑1 (MQ=255)

                                                                                                                                                  ttATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGa                                                                                                             <  4:108005/275‑1 (MQ=255)

                                                                                                                                                  ttATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGGCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCa                                                                                                                                                                                                                                           >  7:1408684/1‑149 (MQ=255)

                                                                                                                                                   tATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaa                                                                                                                                                                                                                                          >  8:613986/1‑149 (MQ=255)

                                                                                                                                                    aTTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaaa                                                                                                                                                                                                                                         <  6:200425/149‑1 (MQ=255)

                                                                                                                                                        cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACt                                                                                                                                                                                                                                                   <  6:1235101/135‑1 (MQ=255)

                                                                                                                                                        cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACt                                                                                                                                                                                                                                                   >  5:1235101/1‑135 (MQ=255)

                                                                                                                                                        cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                                                                                                                                                                                                     >  5:1377451/1‑149 (MQ=255)

                                                                                                                                                         ggCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCa                                                                                                                                                                                                                                    >  5:435332/1‑149 (MQ=255)

                                                                                                                                                                   aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                                                                                                                                                                                                          >  5:237246/1‑149 (MQ=255)

                                                                                                                                                                       tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACGTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcct                                                                                                                                                                                                                      >  6:30744/1‑149 (MQ=255)

                                                                                                                                                                       tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcct                                                                                                                                                                                                                      >  7:375906/1‑149 (MQ=255)

                                                                                                                                                                       tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcct                                                                                                                                                                                                                      >  7:204462/1‑149 (MQ=255)

                                                                                                                                                                       tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcct                                                                                                                                                                                                                      >  8:209539/1‑149 (MQ=255)

                                                                                                                                                                         tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                                                                                                                                                                                                    >  6:1109434/1‑149 (MQ=255)

                                                                                                                                                                            ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGa                                                                                                             <  2:317112/249‑1 (MQ=255)

                                                                                                                                                                             gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                                                                                                                                                                                                <  6:1377451/149‑1 (MQ=255)

                                                                                                                                                                               attaCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCg                                                                                              <  4:3346/258‑1 (MQ=255)

                                                                                                                                                                                   ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCCAATGGCAACTCGCCTGATCCTCTTTCACCACGTATGCTTTGCGTCCCCTTACTATCAGGACGCTTTAGCCCTTGTCCCGCTTTTTGATTTGTAGTTTTTCCCTGGGTTTTCTTTATc                                                                                                                      >  2:290261/1‑233 (MQ=255)

                                                                                                                                                                                    cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTa                                                                                                                                                                                                         >  8:148459/1‑149 (MQ=255)

                                                                                                                                                                                            tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                                                                                                                                                 <  5:1004551/149‑1 (MQ=255)

                                                                                                                                                                                             gtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaa                                                                                                                                                                                                                                                                  >  7:1160130/1‑81 (MQ=37)

                                                                                                                                                                                             gtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaa                                                                                                                                                                                                                                                                  <  8:1160130/83‑3 (MQ=37)

                                                                                                                                                                                                cTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                                                                                                                                                                                             >  7:544258/1‑149 (MQ=255)

                                                                                                                                                                                                      cGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGc                                                                                                                                                                                                                              <  7:770686/110‑1 (MQ=255)

                                                                                                                                                                                                       ggCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACt                                                                                                                                                                                      <  7:621313/149‑1 (MQ=255)

                                                                                                                                                                                                         cAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTAt                                                                                                                                                                                    <  5:1422439/149‑1 (MQ=255)

                                                                                                                                                                                                          aaCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTa                                                                                                                                                                                                         <  8:652422/127‑1 (MQ=255)

                                                                                                                                                                                                          aaCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc                                                                                                                                                                                   >  6:397723/1‑149 (MQ=255)

                                                                                                                                                                                                          aaCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc                                                                                                                                                                                   >  6:968873/1‑149 (MQ=255)

                                                                                                                                                                                                           aCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCa                                                                                                                                                                                  <  5:155768/149‑1 (MQ=255)

                                                                                                                                                                                                           aCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCa                                                                                                                                                                                  <  5:251237/149‑1 (MQ=255)

                                                                                                                                                                                                            cGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAg                                                                                                                                                                                 <  8:1162567/149‑1 (MQ=255)

                                                                                                                                                                                                             gTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAgg                                                                                                                                                                                >  8:742192/1‑149 (MQ=255)

                                                                                                                                                                                                              tCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGa                                                                                                                                                                               <  7:1275562/149‑1 (MQ=255)

                                                                                                                                                                                                               cTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGAc                                                                                                                                                                              >  8:1341964/1‑149 (MQ=255)

                                                                                                                                                                                                                   aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTAACTCGTTGGTTTACAGCGCATCGCCGCCGATCTCAATGCCATCGCAGCGCAGtt                                             >  1:1926/1‑274 (MQ=255)

                                                                                                                                                                                                                      gtgtTACGTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAg                                                                                                                                                                       <  7:1253752/149‑1 (MQ=255)

                                                                                                                                                                                                                      gtgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAg                                                                                                                                                                       <  8:398802/149‑1 (MQ=255)

                                                                                                                                                                                                                         ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGccc                                                                                                                                                                    <  8:544258/149‑1 (MQ=255)

                                                                                                                                                                                                                         ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCttttt                                                                                                                                                      >  4:198976/1‑163 (MQ=255)

                                                                                                                                                                                                                         ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGa                                                                                                         >  4:68107/1‑208 (MQ=255)

                                                                                                                                                                                                                            ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATg                                                                                                                                                                 >  7:8812/1‑149 (MQ=255)

                                                                                                                                                                                                                             cTAAATTTGGCTAAATCCTTCAAGAAGCCGGCCATTCGCGGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGt                                  <  2:132922/275‑1 (MQ=255)

                                                                                                                                                                                                                             cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGt                                                                                                                                                                >  6:742910/1‑149 (MQ=255)

                                                                                                                                                                                                                             cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATg                                                                                                     <  2:68107/208‑1 (MQ=255)

                                                                                                                                                                                                                              tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTc                                                                                                                                                               >  7:859357/1‑149 (MQ=255)

                                                                                                                                                                                                                              tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTc                                                                                                                                                               >  5:669303/1‑149 (MQ=255)

                                                                                                                                                                                                                                aaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTccc                                                                                                                                                             <  7:289695/149‑1 (MQ=255)

                                                                                                                                                                                                                                      gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCtttt                                                                                                                                                       <  5:894100/149‑1 (MQ=255)

                                                                                                                                                                                                                                       cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCttttt                                                                                                                                                      <  6:602838/149‑1 (MQ=255)

                                                                                                                                                                                                                                        tAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttctt                                                                                                                                                                                                                                                                     <  5:680997/37‑1 (MQ=255)

                                                                                                                                                                                                                                        tAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttctt                                                                                                                                                                                                                                                                     >  6:680997/1‑37 (MQ=255)

                                                                                                                                                                                                                                          aaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAt                                                                                                                                                   <  6:1341039/149‑1 (MQ=255)

                                                                                                                                                                                                                                           aTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAtt                                                                                                                                                  <  5:814071/149‑1 (MQ=255)

                                                                                                                                                                                                                                           aTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAtt                                                                                                                                                  <  5:327243/149‑1 (MQ=255)

                                                                                                                                                                                                                                            tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAttt                                                                                                                                                 <  7:724992/149‑1 (MQ=255)

                                                                                                                                                                                                                                            tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAttt                                                                                                                                                 <  6:832165/149‑1 (MQ=255)

                                                                                                                                                                                                                                              cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCACTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGt                                                                                                                                               >  6:1041649/1‑149 (MQ=255)

                                                                                                                                                                                                                                               ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAGATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTa                                                                                                                                              >  8:790472/1‑149 (MQ=255)

                                                                                                                                                                                                                                               ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCtttt                                                                                                                                                       <  7:1184550/140‑1 (MQ=255)

                                                                                                                                                                                                                                               ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACGATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTa                                                                                                                                              <  6:1188303/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  aagaagCCAGCGATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtt                                                                                                                                           <  7:947320/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  aagaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtt                                                                                                                                           >  8:415301/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    gaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt                                                                                                                                         >  6:389214/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt                                                                                                                                         <  8:311337/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt                                                                                                                                         <  8:783611/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      agccagcgAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGGTCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGAAGCTTTAg                                                                                                                                                                       <  6:794208/117‑1 (MQ=255)

                                                                                                                                                                                                                                                        ccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGccc                                                                                                                                     <  6:343252/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCt                                                                                                                                    >  5:44492/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTg                                                                                                                                   >  8:561420/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           gccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGt                                                                                                                                                                                               <  3:26271/88‑1 (MQ=255)

                                                                                                                                                                                                                                                            ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGt                                                                                                                                 <  5:396917/149‑1 (MQ=255)

                                                                                                                                                                                                                                                            ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCAC‑TATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGtt                                                                                                                                <  8:1013089/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             caTTGGCTGGTTTCTTGCCTCTCTGGAAATCACTTATGCCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCg  <  4:301444/275‑1 (MQ=255)

                                                                                                                                                                                                                                                               ttCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGtttt                                                                                                                              >  8:896247/1‑149 (MQ=255)

TCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCG  >  NC_000913/3891249‑3891776

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