BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	591,166	A→G	G412G (GGT→GGC)	nfrB ←	NfrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	591,166	0	A	G	97.4%	474.0 / 1.5	152	G412G (GGT→GGC)	nfrB	NfrB
Reads supporting (aligned to +/- strand): ref base A (2/2); new base G (73/75); total (75/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 158 total aligned reads displayed.

GAAAAGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCA  >  NC_000913/591076‑591253

gAAAAGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGcc                                                                                       >  1:2874452/1‑93 (MQ=255)

  aaaGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCaa                                                                                     >  2:1800737/1‑93 (MQ=255)

     gACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGt                                                                                  >  1:3007097/1‑93 (MQ=255)

     aacacAGGACGAAAAGCAGCCCCAGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGt                                                                                  <  2:1298647/89‑1 (MQ=255)

       cGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGttt                                                                                <  2:1564375/93‑1 (MQ=255)

        gCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTa                                                                               >  2:455314/1‑93 (MQ=255)

        gCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTa                                                                               <  1:2847007/93‑1 (MQ=255)

         cAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAg                                                                              >  2:3048136/1‑93 (MQ=255)

          aGGACGGAAAGCAGACCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATAAAGCCAAAGTTTAGc                                                                             >  1:4067772/1‑93 (MQ=255)

           ggACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGcc                                                                            >  2:772573/1‑93 (MQ=255)

                 aaaGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCAcagca                                                                      >  1:1344814/1‑93 (MQ=255)

                 aaaGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCAcagca                                                                      <  1:934068/93‑1 (MQ=255)

                  aaGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCAcagcag                                                                     <  2:4290437/93‑1 (MQ=255)

                    gcagcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAggg                                                                   >  2:2681514/1‑93 (MQ=255)

                    gcagcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAggg                                                                   <  1:2160198/93‑1 (MQ=255)

                     cagcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGt                                                                  <  2:3424662/93‑1 (MQ=255)

                     cagcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGt                                                                  <  2:3424658/93‑1 (MQ=255)

                     cagcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGt                                                                  <  1:4492502/93‑1 (MQ=255)

                     cagcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGt                                                                  <  1:4376525/93‑1 (MQ=255)

                      agcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTc                                                                 >  2:3564288/1‑93 (MQ=255)

                        ccccTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCAt                                                               <  1:3207152/93‑1 (MQ=255)

                        ccccTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCAt                                                               >  1:2444203/1‑93 (MQ=255)

                         cccTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCAtt                                                              <  2:1701665/93‑1 (MQ=255)

                         cccTGCGTCAGGCAGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCAtt                                                              >  1:1814489/1‑93 (MQ=255)

                           cTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTaa                                                            >  2:4140270/1‑93 (MQ=255)

                            tGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAAc                                                           <  2:2057458/93‑1 (MQ=255)

                             gCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAAcc                                                          <  1:877054/93‑1 (MQ=255)

                              gttCAGGCTGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCa                                                         <  2:3300247/91‑1 (MQ=255)

                              cGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCa                                                         <  2:2952156/93‑1 (MQ=255)

                              cGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCa                                                         <  2:3847719/93‑1 (MQ=255)

                               gTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCAt                                                        >  2:1465360/1‑93 (MQ=255)

                               gTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCAt                                                        >  1:987842/1‑93 (MQ=255)

                               gTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCAt                                                        <  2:2330188/93‑1 (MQ=255)

                                 cAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATgc                                                      <  1:3897257/93‑1 (MQ=255)

                                 cAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATgc                                                      <  1:1965123/93‑1 (MQ=255)

                                 cAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATgc                                                      >  2:3902786/1‑93 (MQ=255)

                                 cAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGCCATTAACCATgc                                                      >  2:200903/1‑93 (MQ=255)

                                   ggCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATgcgc                                                    <  2:3155738/93‑1 (MQ=255)

                                    gCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCt                                                   >  1:2897172/1‑93 (MQ=255)

                                    gCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCt                                                   <  1:667203/93‑1 (MQ=255)

                                     ccGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCtc                                                  >  2:2210759/1‑92 (MQ=255)

                                      cGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGc                                                 <  2:3380252/93‑1 (MQ=255)

                                       gtagtaGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGcc                                                <  1:4474773/93‑1 (MQ=255)

                                       gtagtaGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGcc                                                <  2:1524869/93‑1 (MQ=255)

                                       gtagtaGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGcc                                                >  1:2640295/1‑93 (MQ=255)

                                        tagtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCg                                               >  1:1090828/1‑93 (MQ=255)

                                        tagtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCg                                               >  2:4148746/1‑93 (MQ=255)

                                         agtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGc                                              <  1:393780/93‑1 (MQ=255)

                                          gtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCt                                             <  1:2924110/93‑1 (MQ=255)

                                            agCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGa                                           >  2:1651264/1‑93 (MQ=255)

                                            agCCAGTAACGAAAATCACCCGCGGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGa                                           <  2:291293/93‑1 (MQ=255)

                                             gCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGaa                                          >  2:4684772/1‑93 (MQ=255)

                                              ccAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGaaa                                         <  1:4185765/93‑1 (MQ=255)

                                               cAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGaaaa                                        <  2:3098471/93‑1 (MQ=255)

                                                aGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAAt                                       <  1:50776/93‑1 (MQ=255)

                                                aGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAAt                                       >  1:4070800/1‑93 (MQ=255)

                                                aGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAAt                                       >  1:3906733/1‑93 (MQ=255)

                                                aGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAAt                                       >  2:3166484/1‑93 (MQ=255)

                                                 gTAACGAAAATCGCCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATa                                      <  2:826694/93‑1 (MQ=255)

                                                 gTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATa                                      <  1:1664122/93‑1 (MQ=255)

                                                 gTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATa                                      <  2:3983522/93‑1 (MQ=255)

                                                 gTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATa                                      <  1:3009585/93‑1 (MQ=255)

                                                 gTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATa                                      >  1:941727/1‑93 (MQ=255)

                                                 gTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATa                                      <  2:2746673/93‑1 (MQ=255)

                                                  tAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAg                                     <  2:1724333/93‑1 (MQ=255)

                                                  tAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAg                                     >  1:3217594/1‑93 (MQ=255)

                                                   caCGAAAACCCCCCGCTGCACGATGCGGTCACCCATCAAGCCAAAGATTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGa                                    <  1:1802747/92‑1 (MQ=255)

                                                    aCGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGaa                                   <  2:1150082/93‑1 (MQ=255)

                                                      gAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAg                                 >  1:2114775/1‑93 (MQ=255)

                                                       aaaaTCACCCGCTGCACGATGCGGTTAACCATCAAGCAAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAgg                                <  2:4346836/93‑1 (MQ=255)

                                                        aaaTCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGa                               >  2:2138134/1‑93 (MQ=255)

                                                          aTCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGaaa                             <  2:686521/93‑1 (MQ=255)

                                                          aTCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGaaa                             <  1:3846130/93‑1 (MQ=255)

                                                           tCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAAt                            <  2:3892262/93‑1 (MQ=255)

                                                               ccGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTTATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCa                        <  1:1373799/93‑1 (MQ=255)

                                                               ccGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCa                        <  1:1789453/93‑1 (MQ=255)

                                                                cGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAg                       >  2:3138891/1‑93 (MQ=255)

                                                                cGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAg                       >  2:4693185/1‑93 (MQ=255)

                                                                 gCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAgg                      >  1:1702042/1‑93 (MQ=255)

                                                                 ccTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAgg                      <  1:1516736/92‑1 (MQ=255)

                                                                   tGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCa                    >  1:3112929/1‑93 (MQ=255)

                                                                    gCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCAt                   <  2:3823283/93‑1 (MQ=255)

                                                                    gCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCAt                   >  2:2456159/1‑93 (MQ=255)

                                                                    gCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCAt                   >  1:2286473/1‑93 (MQ=255)

                                                                    gCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCAt                   <  1:4518812/93‑1 (MQ=255)

                                                                      aCGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCg                 >  1:587886/1‑93 (MQ=255)

                                                                       cGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCgg                <  2:2741945/93‑1 (MQ=255)

                                                                        gATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCggg               >  1:835121/1‑93 (MQ=255)

                                                                        gATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCggg               >  1:1528245/1‑93 (MQ=255)

                                                                        gATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCggg               <  1:4016135/93‑1 (MQ=255)

                                                                          tGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGcc             <  2:1199282/93‑1 (MQ=255)

                                                                             ggTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCcaca          >  1:1948710/1‑93 (MQ=255)

                                                                                tAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCAACaaat       >  2:4487241/1‑92 (MQ=255)

                                                                                  aCCAACAAGCCACAGTTAAGCCCAAGCAGGGTAATTAACCATGCGCTGGCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACtt     <  1:1058958/93‑1 (MQ=255)

                                                                                   ccATCAAGCCAAAGTTTAGCCACAGCAGGGTCTTTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACttt    <  1:4534184/93‑1 (MQ=255)

                                                                                   ccATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACttt    <  2:1815506/93‑1 (MQ=255)

                                                                                   ccATAAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACttt    <  1:2599002/93‑1 (MQ=255)

                                                                                     aTCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCa  >  2:4573978/1‑93 (MQ=255)

                                                                                     aTCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCa  <  1:3204569/93‑1 (MQ=255)

                                                                                     aTCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCa  <  2:222961/93‑1 (MQ=255)

GAAAAGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCA  >  NC_000913/591076‑591253

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