BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,265,675	C→A	5.1%	N454K (AAC→AAA)	dnaB →	replicative DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,265,675	0	C	A	5.1%	381.0 / 5.2	112	N454K (AAC→AAA)	dnaB	replicative DNA helicase
Reads supporting (aligned to +/- strand): ref base C (57/48); new base A (3/4); total (60/52)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 114 total aligned reads displayed.

TCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGAC  >  NC_000913/4265526‑4265824

 cTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACg                                                                                                                                                      <  1:677145/150‑1 (MQ=255)

   ggCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGt                                                                                                                                                    <  2:865740/150‑1 (MQ=255)

   ggCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGt                                                                                                                                                    >  2:293625/1‑150 (MQ=255)

    gCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTc                                                                                                                                                   <  1:686041/150‑1 (MQ=255)

         aTCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATgg                                                                                                                                              >  2:1212025/1‑150 (MQ=255)

         aTCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggcg                                                                                                                                              <  1:1195605/150‑5 (MQ=255)

            gAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcg                                                                                                                                           >  2:785297/1‑150 (MQ=255)

            gAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggcgtac                                                                                                                                           <  2:1363120/150‑8 (MQ=255)

            gAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggcgtac                                                                                                                                           <  2:1368649/150‑8 (MQ=255)

                 ggATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCtt                                                                                                                                      >  1:221511/1‑150 (MQ=255)

                   aTGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCg                                                                                                                                    <  2:250740/150‑1 (MQ=255)

                     gCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAAAGGTCGCGCTTCGAc                                                                                                                                  <  1:12210/150‑1 (MQ=255)

                      cGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                                                                                 >  2:452463/1‑150 (MQ=255)

                      cGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                                                                                 >  1:871594/1‑150 (MQ=255)

                       ggACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGAATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCTGGACGGTACGCATGACCTTTAAAGTTCAATGGTCGCGCTTCGACaa                                                                                                                                >  2:576889/1‑150 (MQ=255)

                            tGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATg                                                                                                                           >  1:438363/1‑150 (MQ=255)

                                cATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCggg                                                                                                                       >  1:1119897/1‑150 (MQ=255)

                                  tGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGc                                                                                                                     >  2:346594/1‑150 (MQ=255)

                                  tGTTTATATAACGTGATGAGGTGTATCACGAAAACATTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGCCCTTTAAACGTCAATGGTCGCGCTTCGACAACTATGCtggcc                                                                                                                     >  2:1334119/1‑148 (MQ=255)

                                   gTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGcc                                                                                                                    >  2:35187/1‑150 (MQ=255)

                                     ttATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc                                                                                                                  <  1:785297/150‑1 (MQ=255)

                                     ttATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc                                                                                                                  <  1:1524655/150‑1 (MQ=255)

                                       atctatcGTGATGAGGTTTATCACGAAAACATTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGTTACGCCTTACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAg                                                                                                                <  2:174079/150‑1 (MQ=255)

                                         ctatcGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcgcg                                                                                                                                         <  1:285374/123‑1 (MQ=255)

                                         ctatcGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcgcg                                                                                                                                         >  2:285374/1‑123 (MQ=255)

                                         ctatcGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTa                                                                                                              >  2:481392/1‑150 (MQ=255)

                                           atcGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacg                                                                                                            >  1:174136/1‑150 (MQ=255)

                                                gatgaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgacgac                                                                                                       >  2:632592/1‑150 (MQ=255)

                                                   gaGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaa                                                                                                    <  2:1066209/150‑1 (MQ=255)

                                                    aGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaat                                                                                                   <  1:342337/150‑1 (MQ=255)

                                                     ggTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGATACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaata                                                                                                  >  1:109111/1‑150 (MQ=255)

                                                       tgtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataat                                                                                                >  1:316661/1‑150 (MQ=255)

                                                        gtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataata                                                                                               >  1:1273595/1‑150 (MQ=255)

                                                        gtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataata                                                                                               >  2:616717/1‑150 (MQ=255)

                                                        gtATCACGAAAACAGTCATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataata                                                                                               >  2:249563/1‑150 (MQ=255)

                                                         tATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGTGACGGTACGCCTGACCTTTAACGGTAAATGGTCGCGCTTCTACAACTATGCGGGGCCGCAGTACGACGACGaataataa                                                                                              <  2:70652/150‑1 (MQ=255)

                                                           tCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAatt                                                                                            >  2:332892/1‑150 (MQ=255)

                                                            cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAatta                                                                                           >  2:855641/1‑150 (MQ=255)

                                                            cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAatta                                                                                           >  2:327822/1‑150 (MQ=255)

                                                            cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAatta                                                                                           >  2:1503308/1‑150 (MQ=255)

                                                            cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAatta                                                                                           >  2:1438251/1‑150 (MQ=255)

                                                            cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAatta                                                                                           >  1:155074/1‑150 (MQ=255)

                                                                aaaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTAtttt                                                                                       >  2:693266/1‑150 (MQ=255)

                                                                 aaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGGCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTa                                                                                      <  2:1285564/150‑1 (MQ=255)

                                                                    cAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGa                                                                                   >  2:595482/1‑150 (MQ=255)

                                                                     aGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGaa                                                                                  <  1:53594/150‑1 (MQ=255)

                                                                      gTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAt                                                                                 >  2:468921/1‑150 (MQ=255)

                                                                      gTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAt                                                                                 >  1:132971/1‑150 (MQ=255)

                                                                      gTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAt                                                                                 >  1:1386812/1‑150 (MQ=255)

                                                                      gTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAt                                                                                 >  2:1458069/1‑150 (MQ=255)

                                                                       tGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAtt                                                                                <  1:1350428/150‑1 (MQ=255)

                                                                            tAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATTGTCTCTCTTCGACAACTATGCGGGGCCGCAGTACTACGACGAATAATAATTATTTTATGAATTAGGTa                                                                           <  2:810219/150‑1 (MQ=255)

                                                                            tAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTa                                                                           >  2:806201/1‑150 (MQ=255)

                                                                            tAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTa                                                                           >  2:494830/1‑150 (MQ=255)

                                                                             aaaaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTaa                                                                          >  2:712332/1‑150 (MQ=255)

                                                                              aaaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAAt                                                                         >  1:965275/1‑150 (MQ=255)

                                                                               aaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAAtt                                                                        <  2:221511/150‑1 (MQ=255)

                                                                                    aTCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACgg                                                                                                                                                     >  2:531223/1‑68 (MQ=255)

                                                                                    aTCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACgg                                                                                                                                                     <  1:531223/68‑1 (MQ=255)

                                                                                      cgcgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaa                                                                 <  2:770610/150‑1 (MQ=255)

                                                                                      cgcgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaa                                                                 >  2:783961/1‑150 (MQ=255)

                                                                                       gcgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaaa                                                                >  2:21856/1‑150 (MQ=255)

                                                                                        cgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAac                                                               <  1:35187/150‑1 (MQ=255)

                                                                                        cgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAac                                                               <  2:316661/150‑1 (MQ=255)

                                                                                        cgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAac                                                               <  1:522406/150‑1 (MQ=255)

                                                                                        cgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAac                                                               <  1:346594/150‑1 (MQ=255)

                                                                                         ggAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAaca                                                              >  2:787517/1‑150 (MQ=255)

                                                                                          gAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAacac                                                             >  2:1128280/1‑150 (MQ=255)

                                                                                                attatCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCa                                                       >  1:1092371/1‑150 (MQ=255)

                                                                                                      ggTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAac                                                 >  1:951712/1‑150 (MQ=255)

                                                                                                      ggTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCACAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAac                                                 <  2:890529/150‑1 (MQ=255)

                                                                                                       gTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAaca                                                <  1:1212025/150‑1 (MQ=255)

                                                                                                        tAAACAACGTAACGGCCCAATCGGGACGGTAAGGCTGTCCTTTAACGGTCAATGGTCGCGGTTCGACAACTTTGCGGGGCCGGAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGGAAACACTTATCAAGGAacac                                               <  2:871594/150‑1 (MQ=255)

                                                                                                         aaaCAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAaca                                                              >  2:533809/1‑134 (MQ=255)

                                                                                                         aaaCAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAAcaca                                              >  2:204130/1‑150 (MQ=255)

                                                                                                           acaacGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACaaa                                            >  2:862038/1‑150 (MQ=255)

                                                                                                              acGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGc                                         <  1:481392/150‑1 (MQ=255)

                                                                                                              acGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGc                                         >  1:451834/1‑150 (MQ=255)

                                                                                                                 tAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAg                                      <  1:452463/150‑1 (MQ=255)

                                                                                                                  aaCGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGc                                     <  1:595482/150‑1 (MQ=255)

                                                                                                                      gCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCa                                 >  2:422486/1‑150 (MQ=255)

                                                                                                                        ccAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAAc                               <  1:787517/150‑1 (MQ=255)

                                                                                                                           aTCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtgt                            >  1:1355191/1‑150 (MQ=255)

                                                                                                                            tCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACtgtt                           <  1:158447/150‑1 (MQ=255)

                                                                                                                               ggACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTtgtg                        >  1:128461/1‑150 (MQ=255)

                                                                                                                               ggACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTtgtg                        <  1:1323960/150‑1 (MQ=255)

                                                                                                                                gACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGa                       <  1:1320840/150‑1 (MQ=255)

                                                                                                                                 aCGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGAt                      <  2:486689/150‑1 (MQ=255)

                                                                                                                                    gTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGAtt                     <  2:1001579/148‑1 (MQ=255)

                                                                                                                                    gTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGAtt                     >  1:1001579/1‑148 (MQ=255)

                                                                                                                                    gTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTaa                   <  1:1108710/150‑1 (MQ=255)

                                                                                                                                     tACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAAc                  >  1:301865/1‑150 (MQ=255)

                                                                                                                                      accattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAAcc                 >  2:597292/6‑150 (MQ=255)

                                                                                                                                           tGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACcgccgc            >  2:1462441/1‑150 (MQ=255)

                                                                                                                                             aCCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCgcgc          <  1:293625/150‑1 (MQ=255)

                                                                                                                                               cTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGctc        <  2:451834/150‑1 (MQ=255)

                                                                                                                                                   aaCGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCg                                    <  1:145425/118‑1 (MQ=255)

                                                                                                                                                   aaCGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCg                                    >  2:145425/1‑118 (MQ=255)

                                                                                                                                                   aaCGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCg    <  2:174136/150‑1 (MQ=255)

                                                                                                                                                     cGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGac  <  1:468921/150‑1 (MQ=255)

TCTGGCTCTATCGAGCAGGATGCGGACTTGATCATGTTTATCTATCGTGATGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGGAACACAAATGCAAGCGGCAACTGTTGTGATTAACCGCCGCGCTCTGCGAC  >  NC_000913/4265526‑4265824

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