BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,512	0	A	G	48.8%	211.3 / inf	201	K152K (AAA→AAG)	nohB	bacteriophage DNA packaging protein
Reads supporting (aligned to +/- strand): ref base A (50/53); new base G (45/53); total (95/106)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-06
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 234 total aligned reads displayed.

AGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGC  >  NC_000913/580422‑580607

agttGAAATTGCCACTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGa                                                                                            <  1:1176939/93‑1 (MQ=0)

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

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

                                                   cttcccGGAGCTGGAAACCCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCggg                                         <  1:2936405/92‑1 (MQ=0)

                                                   cttcccGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCggg                                         <  2:12302344/92‑1 (MQ=0)

                                                    ttcccGGAGCTGGAAAACCGACATATTGATTTCCTCACGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCgggg                                        <  2:12231953/93‑1 (MQ=0)

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

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

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

                                                          gagcTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGGTGCTg                                  >  2:616927/4‑96 (MQ=0)

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

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

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

                                                                     ccGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAatatat                       >  2:81197/1‑96 (MQ=2)

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

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

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

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

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

                                                                       ggcATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGCTGCTGAGTGAATATATCg                     >  1:11881950/3‑96 (MQ=11)

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

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

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

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

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

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

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

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

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

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

                                                                        aCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGa                    >  1:9597967/1‑96 (MQ=34)

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

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

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

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

                                                                         cATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCGGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaa                   >  2:6773774/1‑96 (MQ=25)

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

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

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

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

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

                                                                         cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                   <  1:5543789/96‑1 (MQ=11)

                                                                         cATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGaa                   >  1:12213835/1‑96 (MQ=11)

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

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

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

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

                                                                          aTGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAAc                  >  2:8727217/1‑96 (MQ=37)

                                                                           tattGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGATAGAATATATCGAACa                 >  2:8379632/3‑96 (MQ=11)

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

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

                                                                           tGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACa                 >  1:7150515/1‑96 (MQ=38)

                                                                            gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAg                <  1:8666682/96‑1 (MQ=39)

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

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

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

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

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

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

                                                                              tGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc              <  1:741555/96‑1 (MQ=14)

                                                                              cgATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc              >  1:8839427/2‑96 (MQ=14)

                                                                               gATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCa             >  2:9180481/1‑96 (MQ=255)

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

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

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

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

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

                                                                                aTTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATAGCGAACAGTCAg            >  1:5788516/1‑96 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                    ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATTAACTGATACCGGGGTTGCTGAGCGCATATATCGAACAGTCAGGTTa        >  2:1579745/1‑96 (MQ=255)

                                                                                    ccTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTa        <  2:259270/96‑1 (MQ=255)

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

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

                                                                                     cTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaa       <  2:2428948/96‑1 (MQ=255)

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

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

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

                                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg    >  2:1132414/1‑96 (MQ=255)

                                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAg    >  2:11029557/1‑96 (MQ=255)

                                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGGGAATATATCGAACAGTCAGGTTAACAg    >  1:3983769/1‑96 (MQ=255)

                                                                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTACTGAGTGAATATATCGAACAGTCAGGTTAACAg    <  2:11599402/96‑1 (MQ=255)

                                                                                         aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAGTATATCGAACAGTCAGGTTAACAgg   >  1:10750777/1‑96 (MQ=255)

                                                                                         aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg   >  2:4954413/1‑96 (MQ=255)

                                                                                         aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg   <  2:4278371/96‑1 (MQ=255)

                                                                                         aaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAgg   <  2:14151845/96‑1 (MQ=255)

                                                                                          aCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc  <  2:13793896/96‑1 (MQ=255)

                                                                                          aCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGc  >  2:11945497/1‑96 (MQ=255)

AGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGC  >  NC_000913/580422‑580607

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