BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,513	0	C	A	48.2%	360.5 / inf	305	R153R (CGG→AGG) ‡	nohD	DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand): ref base C (77/81); new base A (67/80); total (144/161)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-04
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 319 total aligned reads displayed.

ATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGC  >  NC_000913/580438‑580613

               cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGa                                                              >  1:14720517/4‑101 (MQ=0)

                                  gcttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                           <  1:4220264/96‑1 (MQ=0)

                                  gcttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTg                                           <  1:12507530/96‑1 (MQ=0)

                                   cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGc                                          <  1:7755946/97‑1 (MQ=0)

                                    ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCt                                         <  1:11778593/98‑1 (MQ=1)

                                          gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGGTGCTGAGTGa                                   >  1:630828/4‑101 (MQ=1)

                                                gAAATCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                             <  1:10898400/101‑1 (MQ=11)

                                                 aaaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                            >  1:9965410/1‑101 (MQ=31)

                                                 aaaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                            >  2:6159687/1‑101 (MQ=31)

                                                 aaaaCCGACATATTGTTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                            >  1:13383260/1‑101 (MQ=11)

                                                 aaaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                            <  2:5780982/101‑1 (MQ=11)

                                                  nnnccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                           <  2:9181481/98‑1 (MQ=17)

                                                  nnnccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                           <  2:8431523/98‑1 (MQ=11)

                                                  aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                           <  1:5657184/101‑1 (MQ=32)

                                                  aaaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCg                           >  1:3000928/1‑101 (MQ=17)

                                                  aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                           <  1:6061366/101‑1 (MQ=11)

                                                  aaaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                           <  1:10068588/101‑1 (MQ=11)

                                                   aaCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa                          >  2:94273/1‑101 (MQ=33)

                                                   aaCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                          >  1:5370157/1‑101 (MQ=11)

                                                    aCCGACATATTGATTTCCTCAAGAAGGATATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                         <  1:13682067/101‑1 (MQ=11)

                                                    aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATcnnn                         >  2:12264370/1‑98 (MQ=11)

                                                    aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                         <  1:722215/101‑1 (MQ=11)

                                                    aCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                         <  2:5724343/101‑1 (MQ=11)

                                                     ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc                        <  1:9027107/101‑1 (MQ=35)

                                                      nnncATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                       <  2:14780816/98‑1 (MQ=11)

                                                      cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                       <  1:9222576/101‑1 (MQ=38)

                                                      cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                       >  2:297457/1‑101 (MQ=38)

                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGannn                       >  2:11176917/1‑98 (MQ=11)

                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                       >  1:9908078/1‑101 (MQ=11)

                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                       <  1:12713906/101‑1 (MQ=11)

                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                       <  1:8108679/101‑1 (MQ=11)

                                                      cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCGGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa                       >  1:3691924/1‑101 (MQ=11)

                                                       gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaannn                      >  2:9043237/1‑98 (MQ=25)

                                                       gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                      <  2:5416161/101‑1 (MQ=38)

                                                       ctcATGTTGATTTCATGAAACGGGAAATCATCAACGCCATCAACAAAGCAGCCGAGCTGGATGAACTGATCCCGGGGTTGCTGAGTGAATATATCGAACAg                      <  2:4238117/99‑1 (MQ=14)

                                                        nnntattGATTTCCTCAAGAAGGCCATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt                     <  2:13795528/96‑1 (MQ=14)

                                                        aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGt                     <  1:3487982/101‑1 (MQ=40)

                                                        aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATAGATCGAACAGt                     >  1:13112047/1‑101 (MQ=255)

                                                        aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAgc                     >  1:5833639/1‑100 (MQ=14)

                                                         nnngTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                    <  2:10697805/98‑1 (MQ=255)

                                                         cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                    >  1:6803412/1‑101 (MQ=40)

                                                         cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                    >  1:13909430/1‑101 (MQ=40)

                                                         cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                    <  1:1744323/101‑1 (MQ=40)

                                                         cATGTTGATTTCCTGAAACGGGATATCAGCAACGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATCCCGGGGTTGCTGAGTGAATATATCGAACAGTc                    <  1:5504613/101‑1 (MQ=255)

                                                         cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                    >  1:11544639/1‑101 (MQ=14)

                                                          atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                   >  1:5118609/4‑101 (MQ=14)

                                                          atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                   <  2:6846006/98‑1 (MQ=14)

                                                          aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                   <  1:13455354/101‑1 (MQ=255)

                                                          aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa                   >  2:3894571/1‑101 (MQ=255)

                                                           tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                  <  1:13539510/99‑1 (MQ=14)

                                                           tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                  <  1:1819050/99‑1 (MQ=14)

                                                           tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                  <  2:6437919/99‑1 (MQ=14)

                                                           tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg                  <  1:10350385/99‑1 (MQ=14)

                                                           tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGtnnn                  >  2:7390025/1‑98 (MQ=255)

                                                            nnngATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATGCCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                 <  2:12650397/98‑1 (MQ=14)

                                                            gatGATTTCCTCAAACGGTATATCAACAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                 <  2:872337/99‑1 (MQ=14)

                                                            gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                 <  1:9727896/101‑1 (MQ=255)

                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                 <  1:4504266/100‑1 (MQ=14)

                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                 <  1:12409723/100‑1 (MQ=14)

                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                 <  2:6366051/100‑1 (MQ=14)

                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                 <  2:139798/100‑1 (MQ=14)

                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                 <  1:10969857/100‑1 (MQ=14)

                                                            attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg                 <  1:14337152/100‑1 (MQ=14)

                                                             ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt                <  2:4640238/101‑1 (MQ=14)

                                                             ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt                >  2:1438647/1‑101 (MQ=14)

                                                             ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt                <  1:10560909/101‑1 (MQ=14)

                                                              tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCCGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCCAACAGTCAGGtt               <  1:12603315/101‑1 (MQ=14)

                                                               gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAggnnn              >  2:12085215/1‑98 (MQ=14)

                                                               gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAggnnn              >  2:12870794/1‑98 (MQ=14)

                                                               gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa              <  1:14055473/101‑1 (MQ=17)

                                                               gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa              <  2:3387974/101‑1 (MQ=17)

                                                               gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa              >  1:7333780/1‑101 (MQ=17)

                                                                aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtnnn             >  2:14484031/1‑98 (MQ=255)

                                                                aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGTTTGCTGAGTGAATATATCGAACAGTCAGGTTaa             <  1:14464767/101‑1 (MQ=17)

                                                                 tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc            >  2:1879178/1‑101 (MQ=255)

                                                                 tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc            >  2:953008/1‑101 (MQ=255)

                                                                 tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc            <  1:8516342/101‑1 (MQ=255)

                                                                 nttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc            <  1:2267351/100‑1 (MQ=255)

                                                                  nnncTGAAACGGGATATCATCAAAGCCACGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa           <  2:10390509/98‑1 (MQ=255)

                                                                   tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCTCTAGATGAACTCATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg          <  1:13231470/101‑1 (MQ=255)

                                                                   tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg          <  2:2372361/101‑1 (MQ=255)

                                                                   tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTACTGAGTGAATATATCGAACAGTCAGGTTAACAg          <  1:11944049/101‑1 (MQ=255)

                                                                    ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg         <  1:4413040/101‑1 (MQ=255)

                                                                    ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGCGAATATATCGAACAGTCAGGGTAACAgg         >  1:10155597/1‑101 (MQ=255)

                                                                     cTGAAACGGGATGTCATCAAAGCCATGAACAAAGCAGCCGCGCTGTATGAAATGATATCGGGGTTGCTGAGTGAATATATCGAAAAGTCAGGTTAACannn        >  2:7720875/1‑98 (MQ=17)

                                                                     cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc        <  1:8842837/101‑1 (MQ=255)

                                                                     cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc        >  1:12863638/1‑101 (MQ=255)

                                                                       nnnaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg      <  2:13480561/98‑1 (MQ=255)

                                                                       gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg      <  1:4848505/101‑1 (MQ=255)

                                                                       gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTg      >  1:2504775/1‑101 (MQ=255)

                                                                        nnncGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc     <  2:11561535/98‑1 (MQ=255)

                                                                        acaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGcnnn     >  2:10907706/3‑98 (MQ=255)

                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc     <  2:4763328/101‑1 (MQ=255)

                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc     <  1:8211750/101‑1 (MQ=255)

                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc     <  1:496726/101‑1 (MQ=255)

                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc     <  1:9275113/101‑1 (MQ=255)

                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGc     <  1:3388808/101‑1 (MQ=255)

                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAACAGCCGCGCTGGATGAACTGAGACCGGGGTTGCTTAGTGAAAATATCAAACAGTCAGGATAACAGGCTgg     >  2:3126275/1‑100 (MQ=17)

                                                                         aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg    <  1:700469/101‑1 (MQ=255)

                                                                         aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCg    <  1:2072358/101‑1 (MQ=255)

ATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGC  >  NC_000913/580438‑580613

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