BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,509	0	G	C	50.0%	201.2 / inf	196	L151L (CTG→CTC)	nohB	bacteriophage DNA packaging protein
Reads supporting (aligned to +/- strand): ref base G (47/51); new base C (45/53); total (92/104)
Fisher's exact test for biased strand distribution p-value = 8.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 206 total aligned reads displayed.

TCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACA  >  NC_000913/580418‑580604

                                atgccccTGTCGGTGCAGAGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGcc                                                             >  1:8389948/4‑96 (MQ=0)

                                   cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCtgcg                                                          >  2:14306642/4‑96 (MQ=0)

                                   cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCtgcg                                                          >  1:5898286/4‑96 (MQ=0)

                                                       cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCggg                                      >  2:9622432/5‑96 (MQ=0)

                                                        ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCgggg                                     >  2:12659087/4‑96 (MQ=0)

                                                        ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCgggg                                     <  2:10490969/93‑1 (MQ=0)

                                                        ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCgggc                                     <  1:11271892/93‑2 (MQ=0)

                                                             ggcaCTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCt                                <  1:5522607/93‑1 (MQ=1)

                                                                         ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat                    >  1:6178562/1‑96 (MQ=31)

                                                                         ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat                    <  1:9253283/96‑1 (MQ=31)

                                                                         ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAatatat                    <  2:4501099/96‑1 (MQ=31)

                                                                         ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                    <  2:11846428/96‑1 (MQ=11)

                                                                         ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                    <  1:3021734/96‑1 (MQ=11)

                                                                          cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                   >  1:290336/1‑96 (MQ=32)

                                                                          cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                   <  2:13952342/96‑1 (MQ=32)

                                                                          cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                   <  1:5300820/96‑1 (MQ=32)

                                                                          cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATc                   >  2:1214579/1‑96 (MQ=32)

                                                                          cGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCCGAGTGAATATATc                   <  1:13898476/96‑1 (MQ=17)

                                                                          cGACATATTGATTTCCTCATGAAGGTCATCATTAAAGCCATGATCAAAGCAGCTGCGCTGGCTGCAATAATCCCGGGGTTGCTGAGTGAatatatn                   <  1:1205098/96‑2 (MQ=2)

                                                                          cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGTAATAATACCGGGGTTGCTGAGTGAATATATc                   >  2:4632032/1‑96 (MQ=11)

                                                                          cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                   <  2:2625773/96‑1 (MQ=11)

                                                                          cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                   <  1:5598898/96‑1 (MQ=11)

                                                                          cGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATc                   >  2:5642696/1‑96 (MQ=11)

                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                  >  1:4776672/1‑96 (MQ=33)

                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                  >  2:10617778/1‑96 (MQ=33)

                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                  <  2:8970517/96‑1 (MQ=33)

                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                  <  1:8909584/96‑1 (MQ=33)

                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCg                  >  1:8774641/1‑96 (MQ=33)

                                                                           gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGATTATATCg                  >  1:4326504/1‑96 (MQ=16)

                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                  <  2:5870438/96‑1 (MQ=11)

                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                  <  1:14334549/96‑1 (MQ=11)

                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                  <  2:9774323/96‑1 (MQ=11)

                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                  >  2:1685750/1‑96 (MQ=11)

                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCg                  <  1:8626059/96‑1 (MQ=11)

                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAAGAATACCGGGGGTGCTGAGTGAATATATCg                  >  1:10411990/1‑96 (MQ=2)

                                                                           gACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATCAAATAATACCGGGGTTGCTGAGTGAATATATCg                  <  2:9829544/96‑1 (MQ=2)

                                                                            aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                 <  2:1741192/96‑1 (MQ=11)

                                                                            aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                 <  2:12801911/96‑1 (MQ=11)

                                                                            aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                 >  2:1225152/1‑96 (MQ=11)

                                                                            aCATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGa                 <  1:11794217/96‑1 (MQ=11)

                                                                             cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                >  2:13520927/1‑96 (MQ=35)

                                                                             cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                >  2:706005/1‑96 (MQ=35)

                                                                             cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                <  1:1748369/96‑1 (MQ=35)

                                                                             cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                <  2:9404469/96‑1 (MQ=35)

                                                                             cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                <  2:13045783/96‑1 (MQ=35)

                                                                             cATATTGATTTCGTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                >  1:7373227/1‑96 (MQ=11)

                                                                             cATATTGATTTCCTCAAGAAGGATATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                <  2:13298159/96‑1 (MQ=11)

                                                                             cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                >  1:9847076/1‑96 (MQ=2)

                                                                              atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc               <  2:3976798/93‑1 (MQ=11)

                                                                              atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc               <  1:7095645/93‑1 (MQ=11)

                                                                              atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc               >  2:4952824/4‑96 (MQ=11)

                                                                              atattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAAc               >  2:5459848/4‑96 (MQ=11)

                                                                              aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCTCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc               <  2:12865006/96‑1 (MQ=37)

                                                                              aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc               >  1:496869/1‑96 (MQ=37)

                                                                              aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAAc               <  2:8758841/96‑1 (MQ=37)

                                                                               tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACa              <  2:12345892/94‑1 (MQ=11)

                                                                               tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa              <  2:8949719/96‑1 (MQ=38)

                                                                               tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa              <  1:12314381/96‑1 (MQ=38)

                                                                               tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa              <  1:5560989/96‑1 (MQ=38)

                                                                                gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg             >  1:11034744/1‑96 (MQ=39)

                                                                                gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg             >  2:8079739/1‑96 (MQ=39)

                                                                                gTTGATTTCATGAAACGGGAAATCATCAACGCCATCAACAAAGCAGCCGAGCTGGATGAACTGATCCCGGGGTTGCTGAGTGAATATATCGAACAg             <  1:4111865/96‑1 (MQ=14)

                                                                                attGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAg             <  2:12411089/95‑1 (MQ=14)

                                                                                 ttGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGt            <  2:3389137/96‑1 (MQ=40)

                                                                                 ttGATTTCCTCAAGAAGGCCATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt            <  1:13409371/96‑1 (MQ=14)

                                                                                 ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt            <  2:8685417/96‑1 (MQ=14)

                                                                                 ttGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGt            >  1:1402386/1‑96 (MQ=14)

                                                                                  tGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTc           >  2:5652033/1‑96 (MQ=255)

                                                                                  tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc           >  1:1941138/1‑96 (MQ=14)

                                                                                   gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa          >  1:11734283/1‑96 (MQ=14)

                                                                                   gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa          >  1:12500414/1‑96 (MQ=14)

                                                                                   gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa          <  1:6627474/96‑1 (MQ=14)

                                                                                   gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa          >  2:7095203/1‑96 (MQ=14)

                                                                                   gATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa          >  2:7307854/1‑96 (MQ=14)

                                                                                    aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg         <  1:11512283/96‑1 (MQ=255)

                                                                                    aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg         >  2:12296472/1‑96 (MQ=255)

                                                                                    aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg         >  1:14081801/1‑96 (MQ=255)

                                                                                    aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTAAg         >  2:6824741/1‑96 (MQ=255)

                                                                                    aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAACGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg         <  2:12224460/96‑1 (MQ=14)

                                                                                    aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg         >  2:9539552/1‑96 (MQ=14)

                                                                                    aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg         <  2:10052007/96‑1 (MQ=14)

                                                                                    aTTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAg         <  2:13158339/96‑1 (MQ=14)

                                                                                     tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg        >  1:3593740/1‑96 (MQ=255)

                                                                                     tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg        >  1:931506/1‑96 (MQ=255)

                                                                                     tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg        <  2:6673310/96‑1 (MQ=255)

                                                                                     tttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAgg        <  2:6078989/96‑1 (MQ=255)

                                                                                      ttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGt       <  1:13159386/96‑1 (MQ=255)

                                                                                      ttCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGATGCTGGGAGAATATATCGAACAGTCAGGt       >  1:6208402/1‑96 (MQ=255)

                                                                                       tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt      <  1:6101574/96‑1 (MQ=255)

                                                                                       tCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGtt      >  2:2008055/1‑96 (MQ=255)

                                                                                        ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa     >  1:8196976/1‑96 (MQ=255)

                                                                                        ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAACAGTCAGGTTa     >  2:1884343/1‑96 (MQ=255)

                                                                                         cTGAAACGGGATGTCATCAAAGCCATGAACAAAGCAGCCGCGCTGTATGAAATGATATCGGGGTTGCTGAGTGAATATATCGAAAAGTCAGGTTaa    >  1:7476453/1‑96 (MQ=14)

                                                                                         cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa    >  2:1132493/1‑96 (MQ=255)

                                                                                          tGAAACGGGATATCATCAAAGCCCTGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc   <  1:10352604/96‑1 (MQ=255)

                                                                                          tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc   <  2:8259191/96‑1 (MQ=255)

                                                                                          tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc   <  2:2200517/96‑1 (MQ=255)

                                                                                          tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAc   >  1:8983927/1‑96 (MQ=255)

                                                                                           gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa  >  2:2431651/1‑96 (MQ=255)

                                                                                           gAAACGGGATATCATCAAAGCCACGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACa  <  1:10091696/96‑1 (MQ=255)

TCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACA  >  NC_000913/580418‑580604

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