BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	18,885	A→G	G57G (GGA→GGG)	nhaR →	NhaR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	18,885	0	A	G	99.7%	979.2 / ‑5.1	314	G57G (GGA→GGG)	nhaR	NhaR
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (152/161); total (153/161)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 321 total aligned reads displayed.

GGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGT  >  NC_000913/18795‑18974

ggAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCg                                                                                         <  2:3823685/93‑1 (MQ=255)

ggAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCg                                                                                         >  2:49147/1‑93 (MQ=255)

 gAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGt                                                                                        >  2:2717756/1‑93 (MQ=255)

 gAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGt                                                                                        >  1:7310172/1‑93 (MQ=255)

 gAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGt                                                                                        <  1:6983432/93‑1 (MQ=255)

   ggCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTgg                                                                                      >  2:593849/1‑93 (MQ=255)

      gcTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGtct                                                                                   >  2:1350455/1‑93 (MQ=255)

      gcTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGtct                                                                                   >  2:5282553/1‑93 (MQ=255)

      gcTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGtct                                                                                   >  1:7691119/1‑93 (MQ=255)

      gcTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGtcc                                                                                   >  1:6054030/1‑92 (MQ=255)

       cTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGtctc                                                                                  <  1:6502532/93‑1 (MQ=255)

          tatttaACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGac                                                                               >  1:2682517/1‑92 (MQ=255)

            tttaACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAAcc                                                                             >  2:1608090/1‑93 (MQ=255)

            tttaACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAAcc                                                                             >  2:3243492/1‑93 (MQ=255)

                aCTCCACAAACCATTACAGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGc                                                                         <  1:7021566/93‑1 (MQ=255)

                 cTCCACAATCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCg                                                                        <  2:1383142/93‑1 (MQ=255)

                 cTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCg                                                                        <  2:7257574/93‑1 (MQ=255)

                  tCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGa                                                                       >  2:6741505/1‑93 (MQ=255)

                  accACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGa                                                                       >  2:4794631/2‑93 (MQ=255)

                   ccACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAg                                                                      <  1:5387130/93‑1 (MQ=255)

                    cacaAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGc                                                                     >  2:2482001/1‑93 (MQ=255)

                    cacaAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGc                                                                     <  2:469677/93‑1 (MQ=255)

                     acaAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCt                                                                    >  1:100472/1‑93 (MQ=255)

                     acaAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCt                                                                    >  2:6048574/1‑93 (MQ=255)

                       aaaCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTgg                                                                  <  1:5178383/93‑1 (MQ=255)

                         aCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGga                                                                <  2:6370732/93‑1 (MQ=255)

                         aCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGga                                                                <  1:4548425/93‑1 (MQ=255)

                         aCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGga                                                                <  2:3319928/93‑1 (MQ=255)

                         aCCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGga                                                                <  1:2352821/93‑1 (MQ=255)

                           cATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGgaga                                                              >  2:242289/1‑93 (MQ=255)

                            aTTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGaa                                                             <  1:7097462/93‑1 (MQ=255)

                            aTTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGaa                                                             >  2:1448286/1‑93 (MQ=255)

                            aTTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGaa                                                             >  1:1343854/1‑93 (MQ=255)

                            aTTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGaa                                                             <  2:1843211/93‑1 (MQ=255)

                            aTTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGaa                                                             <  2:1874263/93‑1 (MQ=255)

                             ttACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAAc                                                            <  1:953232/93‑1 (MQ=255)

                                 cagACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGt                                                        >  2:4102708/3‑93 (MQ=255)

                                 cGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGt                                                        <  2:29400/93‑1 (MQ=255)

                                  ggACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTc                                                       >  1:3935471/1‑93 (MQ=255)

                                      aGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTAtg                                                   >  1:5041173/1‑92 (MQ=255)

                                       gATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCg                                                  >  1:1369288/1‑93 (MQ=255)

                                          tCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTa                                               <  2:5599983/93‑1 (MQ=255)

                                           cGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTAt                                              <  1:6690494/93‑1 (MQ=255)

                                            gAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATg                                             <  2:3545046/93‑1 (MQ=255)

                                             aGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGc                                            <  1:4302762/93‑1 (MQ=255)

                                              gcgcTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGcc                                           <  2:5461198/93‑1 (MQ=255)

                                                gcTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGa                                         <  1:4770338/93‑1 (MQ=255)

                                                 cTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGAt                                        <  1:150294/93‑1 (MQ=255)

                                                 cTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGAt                                        <  2:368061/93‑1 (MQ=255)

                                                  tGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATa                                       <  2:1718714/93‑1 (MQ=255)

                                                  tGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATa                                       >  1:2976255/1‑93 (MQ=255)

                                                   ggAAGAGCGCCTGCAAGGCAAATTTTTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATaa                                      <  1:6174384/93‑1 (MQ=255)

                                                   ggAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATaa                                      <  2:6184730/93‑1 (MQ=255)

                                                     aaGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATaaaa                                    <  1:7462766/93‑1 (MQ=255)

                                                      agagCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAAt                                   <  2:901299/93‑1 (MQ=255)

                                                      agagCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAAt                                   <  2:1188922/93‑1 (MQ=255)

                                                      agagCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAAt                                   <  1:3566998/93‑1 (MQ=255)

                                                       gagCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATg                                  >  1:4470012/1‑93 (MQ=255)

                                                       gagCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATg                                  >  1:4969176/1‑93 (MQ=255)

                                                       gagCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATg                                  <  2:6943345/93‑1 (MQ=255)

                                                        agCGCCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGt                                 <  1:4777807/93‑1 (MQ=255)

                                                         gcgcCTGCAAGGCAAATTATTCAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGtt                                <  1:560812/93‑1 (MQ=255)

                                                          cgcCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTc                               >  1:2540641/1‑93 (MQ=255)

                                                           gcCTGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCa                              >  1:1734413/1‑93 (MQ=255)

                                                              tGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCt                           >  1:5252929/1‑93 (MQ=255)

                                                              tGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCt                           >  1:1852193/1‑93 (MQ=255)

                                                              tGCAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCt                           >  2:6460280/1‑93 (MQ=255)

                                                                gaaGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGACCTGGTCTATCGCTATGCCGATAAAATGTTCACCTTa                         <  1:3401719/92‑1 (MQ=255)

                                                                cAAGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTa                         <  1:6632754/93‑1 (MQ=255)

                                                                 aaGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTaa                        >  2:5706606/1‑93 (MQ=255)

                                                                 aaGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTaa                        >  2:3934940/1‑93 (MQ=255)

                                                                  aGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAg                       >  2:4382921/1‑93 (MQ=255)

                                                                  aGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAg                       >  2:4530783/1‑93 (MQ=255)

                                                                  aGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAg                       <  1:4457458/93‑1 (MQ=255)

                                                                  aGGCAAATTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAg                       >  2:3442941/1‑93 (MQ=255)

                                                                     cAAATTATTTAAACGCAAGGGGAGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCa                    >  1:4214662/1‑93 (MQ=255)

                                                                      aaaTTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAg                   >  2:5781777/1‑93 (MQ=255)

                                                                       aaTTATTTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAgg                  >  1:282304/1‑93 (MQ=255)

                                                                       aaTTATTTAAACGCAAGGGGCGTGGTCTCGAAACCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAgg                  >  2:7385810/1‑93 (MQ=255)

                                                                         ttattTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGaa                <  1:3164763/93‑1 (MQ=255)

                                                                         ttattTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGaa                >  1:2925891/1‑93 (MQ=255)

                                                                         ttattTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGACTATCGCTATGCGGATAAAATGTTCACTTTAAGCCAGGaa                >  1:2189755/1‑93 (MQ=255)

                                                                          tattTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGaaa               <  2:4138586/93‑1 (MQ=255)

                                                                          tattTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGaaa               <  1:6336586/93‑1 (MQ=255)

                                                                           attTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAAt              >  1:3891026/1‑93 (MQ=255)

                                                                           attTAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAAt              >  1:208696/1‑93 (MQ=255)

                                                                            tttAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATg             >  2:7886068/1‑93 (MQ=255)

                                                                            tttAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATg             >  2:768404/1‑93 (MQ=255)

                                                                            tttAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATg             >  1:720132/1‑93 (MQ=255)

                                                                            tttAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATg             >  1:2846217/1‑93 (MQ=255)

                                                                            tttAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATg             >  1:3715130/1‑93 (MQ=255)

                                                                            tttAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATg             >  1:5518042/1‑93 (MQ=255)

                                                                              tAAACGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCt           <  2:6793426/93‑1 (MQ=255)

                                                                                aaCGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTgg         <  2:5193613/93‑1 (MQ=255)

                                                                                 aCGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGa        <  1:5163911/93‑1 (MQ=255)

                                                                                 aCGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGa        <  1:5944402/93‑1 (MQ=255)

                                                                                 aCGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGa        <  1:6251585/93‑1 (MQ=255)

                                                                                  cGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGat       >  2:2187840/1‑93 (MQ=255)

                                                                                  cGCAAGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGat       <  1:1173289/93‑1 (MQ=255)

                                                                                      aGGGGCGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATtg   >  1:4645982/1‑93 (MQ=255)

GGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGT  >  NC_000913/18795‑18974

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