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	99.4%	1123.8 / ‑5.7	343	G412G (GGT→GGC)	nfrB	NfrB
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (166/175); total (167/176)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 361 total aligned reads displayed.

GAAAAGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATA  >  NC_000913/591076‑591255

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

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

  aaaGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCaa                                                                                       <  2:827268/93‑1 (MQ=255)

   aaGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCaaa                                                                                      >  1:6355144/1‑93 (MQ=255)

   aaGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCaaa                                                                                      <  2:289382/93‑1 (MQ=255)

   aaGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCaaa                                                                                      <  1:7484790/93‑1 (MQ=255)

      aCGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGtt                                                                                   <  1:4364471/93‑1 (MQ=255)

       cGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGttt                                                                                  >  1:5054274/1‑93 (MQ=255)

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

          aGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGc                                                                               >  2:2478551/1‑93 (MQ=255)

           gggcGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGcc                                                                              <  2:8635418/90‑1 (MQ=255)

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

             aCGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCCCCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCcac                                                                            <  1:7336461/93‑1 (MQ=255)

             aCGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCcac                                                                            >  1:5984412/1‑93 (MQ=255)

             aCGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCcac                                                                            <  1:5681890/93‑1 (MQ=255)

              cGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCcaca                                                                           <  2:5152168/93‑1 (MQ=255)

              cGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCcaca                                                                           <  1:9423984/93‑1 (MQ=255)

                gAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCAcagc                                                                         <  1:4176210/93‑1 (MQ=255)

                gAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCAcagc                                                                         >  1:6507678/1‑93 (MQ=255)

                  aaGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCAcagcag                                                                       >  2:5402482/1‑93 (MQ=255)

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

                   agcaccCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAgg                                                                      <  2:5716595/93‑1 (MQ=255)

                       gcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCa                                                                  >  1:902803/1‑93 (MQ=255)

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

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

                         cccTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACAATAAAGCCAAAGTTTAGCCAAAGCAGGGTCAta                                                                >  1:2101934/1‑92 (MQ=255)

                          ccTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTa                                                               >  1:6697365/1‑93 (MQ=255)

                           cTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTac                                                              >  1:1090639/1‑92 (MQ=255)

                            tGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAAc                                                             >  1:6268807/1‑93 (MQ=255)

                            tGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAAc                                                             >  1:9191560/1‑93 (MQ=255)

                            tGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAAc                                                             >  2:3455771/1‑93 (MQ=255)

                            tGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAAc                                                             >  2:7349576/1‑93 (MQ=255)

                            tGCGGCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAAc                                                             >  1:9935734/1‑93 (MQ=255)

                              ccccAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCa                                                           <  1:8599883/90‑1 (MQ=255)

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

                              cGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCAAAAGTTTAGCCACAGCAGGGTCATTAACCa                                                           <  2:6367/93‑1 (MQ=255)

                               gTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCAt                                                          <  2:8144712/93‑1 (MQ=255)

                                tacggCTGTATTAGCAAGTAACGAAACCCCGCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATg                                                         <  2:1947342/90‑1 (MQ=255)

                                tCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATg                                                         >  1:6959636/1‑93 (MQ=255)

                                tCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATg                                                         <  1:2496354/93‑1 (MQ=255)

                                tCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATg                                                         >  1:4970861/1‑93 (MQ=255)

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

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

                                   gcccGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATgcgc                                                      <  1:479263/91‑1 (MQ=255)

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

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

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

                                         agtagCCAGTAGCGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGc                                                <  2:4217399/93‑1 (MQ=255)

                                         agtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGc                                                >  2:895354/1‑93 (MQ=255)

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

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

                                          gtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCt                                               >  2:1073170/1‑93 (MQ=255)

                                           tagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTg                                              >  1:7786966/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                    aCGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGaa                                     >  2:3574483/1‑93 (MQ=255)

                                                    aCGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGCTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGaa                                     >  1:3756487/1‑93 (MQ=255)

                                                     cGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGaaa                                    <  2:3144062/93‑1 (MQ=255)

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

                                                       aaaaTCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAgg                                  >  1:8238984/1‑93 (MQ=255)

                                                          aTCGCCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGaaa                               <  1:8679955/93‑1 (MQ=255)

                                                           tcccccGCTGCACGATGCGGTTACCCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAAt                              <  2:4978234/90‑1 (MQ=255)

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

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

                                                                cGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAg                         <  1:9642882/93‑1 (MQ=255)

                                                                 gCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAgg                        <  2:4429139/93‑1 (MQ=255)

                                                                 gCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAgg                        <  2:8593701/93‑1 (MQ=255)

                                                                 gCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAgg                        <  1:5413352/93‑1 (MQ=255)

                                                                  cTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGc                       <  2:1465481/93‑1 (MQ=255)

                                                                   tGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCa                      <  1:3427710/93‑1 (MQ=255)

                                                                    gCAGGATGCTGTTAACCATCAAGCGAAAGTTTAGCCACAGCAGGGTCATTACCCAGGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCAt                     <  2:5647008/93‑1 (MQ=255)

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

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

                                                                     cACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATc                    <  2:7165832/93‑1 (MQ=255)

                                                                       cGATGCGGTGAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCgg                  <  1:1642337/93‑1 (MQ=255)

                                                                          tGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGcc               >  2:922749/1‑93 (MQ=255)

                                                                            cGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCcac             <  1:7060961/93‑1 (MQ=255)

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

                                                                              gTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACaa           >  1:397866/1‑93 (MQ=255)

                                                                              gTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACaa           >  1:9258178/1‑93 (MQ=255)

                                                                               ttAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACaaa          >  1:47035/1‑93 (MQ=255)

                                                                               ttAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACaaa          <  2:6761800/93‑1 (MQ=255)

                                                                                tAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAAc         >  2:5214589/1‑93 (MQ=255)

                                                                                 aaCCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACt        <  2:4135165/93‑1 (MQ=255)

                                                                                  aCCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACtt       >  2:3943364/1‑93 (MQ=255)

                                                                                   ccATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGATAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACttt      >  2:5268458/1‑93 (MQ=255)

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

                                                                                    cATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTc     >  1:8013216/1‑93 (MQ=255)

                                                                                    cATCAAGCAAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTc     <  1:2328550/93‑1 (MQ=255)

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

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

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

                                                                                      tCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCAt   <  2:3666479/93‑1 (MQ=255)

                                                                                       cAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATa  <  1:2546103/93‑1 (MQ=255)

                                                                                       cAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATa  <  2:3026193/93‑1 (MQ=255)

GAAAAGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATA  >  NC_000913/591076‑591255

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