BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	216,444	C→T	S187N (AGC→AAC)	yaeF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	216,444	0	C	T	99.4%	1226.2 / ‑4.8	346	S187N (AGC→AAC)	yaeF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (170/174); total (172/174)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 356 total aligned reads displayed.

TTTGGCAAAGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGCTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCG  >  NC_000913/216353‑216535

  tGGCAAAGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGa                                                                                          <  1:555955/93‑1 (MQ=255)

  tGGCAAAGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGa                                                                                          <  2:7880222/93‑1 (MQ=255)

  tGGCAAAGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGa                                                                                          <  1:357649/93‑1 (MQ=255)

  tGGCAAAGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGa                                                                                          <  2:3917121/93‑1 (MQ=255)

       aaGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCAc                                                                                     <  2:6856103/93‑1 (MQ=255)

        aGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACt                                                                                    <  1:2142130/93‑1 (MQ=255)

          gCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTTACTGCGCTTTCGCCAGGCCGTTGACGCACTGt                                                                                  <  1:3956895/93‑1 (MQ=255)

          gCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGt                                                                                  <  2:7669310/93‑1 (MQ=255)

            aTCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTg                                                                                >  2:2938142/1‑93 (MQ=255)

             tCCGTGACAAATTCCGAACAAAACCACGATTTTTTGTCTCCTTCGCCTACACTGCTTAACTGCGCTTTCGCCAGGCCGTGGACGCACTGTTgg                                                                               >  1:6993847/1‑93 (MQ=255)

             tCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTgg                                                                               >  1:5056540/1‑93 (MQ=255)

             tCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTgg                                                                               >  1:6781187/1‑93 (MQ=255)

             tCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTgg                                                                               >  2:1178878/1‑93 (MQ=255)

              ccGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGc                                                                              >  2:2659061/1‑93 (MQ=255)

              ccGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGc                                                                              >  2:6122564/1‑93 (MQ=255)

              ccGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGc                                                                              <  2:1711772/93‑1 (MQ=255)

              ccGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGc                                                                              >  1:7897952/1‑93 (MQ=255)

                gTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGa                                                                            >  2:2402122/1‑93 (MQ=255)

                  gACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGaaa                                                                          >  2:3545996/1‑93 (MQ=255)

                   aCAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGaaaa                                                                         >  1:124141/1‑93 (MQ=255)

                   aCAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGaaaa                                                                         >  1:1144209/1‑93 (MQ=255)

                    cAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAAt                                                                        <  2:6250467/93‑1 (MQ=255)

                     aaaTTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATc                                                                       >  1:2870677/1‑93 (MQ=255)

                     aaaTTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATc                                                                       <  2:2110818/93‑1 (MQ=255)

                      aaTTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATcc                                                                      <  1:3614173/93‑1 (MQ=255)

                      aaTTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATcc                                                                      <  2:701240/93‑1 (MQ=255)

                       aTTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCt                                                                     >  2:6390717/1‑93 (MQ=255)

                       aTTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCt                                                                     >  2:946109/1‑93 (MQ=255)

                       aTTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCt                                                                     >  1:2254701/1‑93 (MQ=255)

                        ttCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTc                                                                    <  2:3122206/93‑1 (MQ=255)

                        ttCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTc                                                                    >  1:6517892/1‑93 (MQ=255)

                          ccGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCgg                                                                  >  2:6207182/1‑93 (MQ=255)

                           cGAACAAAGCCACGACTTTTTGTCTCCTTCGCCCACGCTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGga                                                                 <  1:1964558/93‑1 (MQ=255)

                            gAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGgag                                                                <  1:7809546/93‑1 (MQ=255)

                              aCAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGaa                                                              >  2:7562097/1‑93 (MQ=255)

                               cAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAAc                                                             <  2:231711/93‑1 (MQ=255)

                               cAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAAc                                                             <  1:7227381/93‑1 (MQ=255)

                                aaaaCCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACg                                                            >  2:4429898/1‑93 (MQ=255)

                                 aaaCCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACgg                                                           >  1:4210867/1‑93 (MQ=255)

                                 aaaCCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACgg                                                           >  1:5420615/1‑93 (MQ=255)

                                  aaCCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGa                                                          <  2:667169/93‑1 (MQ=255)

                                    ccACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGAtt                                                        >  1:5935682/1‑93 (MQ=255)

                                     cACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTc                                                       <  2:1552641/93‑1 (MQ=255)

                                      aCGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCa                                                      >  2:5962729/1‑93 (MQ=255)

                                       cGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAg                                                     >  2:3799632/1‑93 (MQ=255)

                                       cGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAg                                                     <  1:3674800/93‑1 (MQ=255)

                                        gACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGt                                                    <  1:4191750/93‑1 (MQ=255)

                                         aCTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTg                                                   <  1:2838394/93‑1 (MQ=255)

                                          cTTTTTGTCTCCTTCGCCCACACTGCTTATCTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGa                                                  <  2:893702/93‑1 (MQ=255)

                                          cTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGa                                                  <  1:3094993/93‑1 (MQ=255)

                                          cTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGa                                                  <  1:744585/93‑1 (MQ=255)

                                           tttttGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTAGGAGACAGGATTCAGTGAg                                                 >  1:6142312/1‑93 (MQ=255)

                                              ttGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGACAATCCTCGGAGAACGGATTCAGTGAGcac                                              >  1:1738633/1‑93 (MQ=255)

                                                gTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACAt                                            <  1:5181126/93‑1 (MQ=255)

                                                 tctcCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATc                                           <  1:7747998/93‑1 (MQ=255)

                                                 tctcCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATc                                           >  2:5840230/1‑93 (MQ=255)

                                                   tcCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTg                                         >  1:4016011/1‑93 (MQ=255)

                                                   tcCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTg                                         <  2:4939295/93‑1 (MQ=255)

                                                    ccTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTgg                                        >  1:873462/1‑93 (MQ=255)

                                                        cGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAg                                    >  1:347317/1‑93 (MQ=255)

                                                        cGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAg                                    <  1:7416894/93‑1 (MQ=255)

                                                         gCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGt                                   <  2:7142418/93‑1 (MQ=255)

                                                         gCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTCGCGAGt                                   >  1:4100344/1‑93 (MQ=255)

                                                          cccACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTc                                  <  2:7428028/93‑1 (MQ=255)

                                                           ccACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCa                                 <  2:2031626/93‑1 (MQ=255)

                                                           ccACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCa                                 <  2:3449709/93‑1 (MQ=255)

                                                            cacaCTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAc                                >  2:4792103/1‑93 (MQ=255)

                                                            cacaCTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAc                                >  2:6559739/1‑93 (MQ=255)

                                                            cacaCTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAc                                >  2:7259418/1‑93 (MQ=255)

                                                            cacaCTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAc                                >  1:4395858/1‑93 (MQ=255)

                                                             acacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAcc                               <  1:3172962/93‑1 (MQ=255)

                                                             acacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAcc                               >  1:3252868/1‑93 (MQ=255)

                                                             acacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAcc                               >  1:661042/1‑93 (MQ=255)

                                                             acacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAcc                               >  1:6353328/1‑93 (MQ=255)

                                                              cacTGCTTAACTGCGCTTTTGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCa                              >  1:962683/1‑93 (MQ=255)

                                                              cacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCa                              <  2:1738009/93‑1 (MQ=255)

                                                               acTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCAt                             <  1:1094148/93‑1 (MQ=255)

                                                               acTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCAt                             <  1:4840162/93‑1 (MQ=255)

                                                               acTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCAt                             <  2:3731511/93‑1 (MQ=255)

                                                                cTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATa                            <  2:4264617/93‑1 (MQ=255)

                                                                cTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATa                            <  2:4030118/93‑1 (MQ=255)

                                                                 tGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATaa                           <  1:4605118/93‑1 (MQ=255)

                                                                    ttAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAgg                        <  2:1854169/93‑1 (MQ=255)

                                                                    ttAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAgg                        <  2:6710026/93‑1 (MQ=255)

                                                                    ttAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAgg                        >  2:135191/1‑93 (MQ=255)

                                                                         tGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATa                   >  2:4127783/1‑93 (MQ=255)

                                                                           cgcTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATaaa                 >  2:5556423/1‑93 (MQ=255)

                                                                              tttCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTc              <  1:3294091/93‑1 (MQ=255)

                                                                              tttCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTc              >  2:4540235/1‑93 (MQ=255)

                                                                              tttCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTc              >  1:1500908/1‑93 (MQ=255)

                                                                                 cGCCAGGCCGTTGATGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGAc           <  1:2097501/93‑1 (MQ=255)

                                                                                  gCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACa          <  1:2509443/93‑1 (MQ=255)

                                                                                  gCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACa          <  1:3741513/93‑1 (MQ=255)

                                                                                  gCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACa          >  1:7573024/1‑93 (MQ=255)

                                                                                    cAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAAt        <  1:5793428/93‑1 (MQ=255)

                                                                                    cAGGCCGTTGACGCACTGGTGGCGAAAATCCTCGGAGAACGGATTCAGTTAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAAt        <  1:2715429/93‑1 (MQ=255)

                                                                                     aGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATg       >  2:2241958/1‑93 (MQ=255)

                                                                                      ggCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGc      >  1:6214289/1‑93 (MQ=255)

                                                                                      ggCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGc      <  2:1485022/93‑1 (MQ=255)

                                                                                      ggCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGc      >  2:5282356/1‑93 (MQ=255)

TTTGGCAAAGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGCTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCG  >  NC_000913/216353‑216535

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