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.3%	881.5 / ‑3.6	276	G412G (GGT→GGC)	nfrB	NfrB
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (138/136); total (139/137)
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.41e-01

Only 100 of 282 total aligned reads displayed.

GAAAAGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATA  >  NC_000913/591076‑591255

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

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

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

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

   aaGACGAAGGAGGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCaaa                                                                                      <  2:7555977/93‑1 (MQ=255)

    aGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAg                                                                                     >  1:4069086/1‑93 (MQ=255)

    aGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAg                                                                                     <  1:7764026/93‑1 (MQ=255)

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

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

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

     gACGCAGGACGGAAAGCAGCCCCGGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGt                                                                                    <  1:1189493/93‑1 (MQ=255)

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

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

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

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

            gACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCca                                                                             >  1:1264739/1‑93 (MQ=255)

            gACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCca                                                                             >  2:2996037/1‑93 (MQ=255)

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

              cGGAAAGCAGCCCCTTCGTCAGGCCGTAGTAGCCAGTAACGAAAATCCCCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCcaca                                                                           <  1:1130654/93‑1 (MQ=255)

              cGGAAAGCAGCCCCTGCATCAGCCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCTCTATGCGGTTAGCCATCAAGCCAAAGTTTAGCcaca                                                                           <  1:4310396/93‑1 (MQ=255)

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

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

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

                    ccgtccccTGCGTGAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAggg                                                                     <  2:446997/89‑1 (MQ=255)

                     cagcCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGt                                                                    >  1:4755858/1‑93 (MQ=255)

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

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

                        ccccTGTGACAGGCCGTAGTAGCCAGTAACGAAAAGCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCAt                                                                 <  2:6419623/93‑1 (MQ=255)

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

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

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

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

                          ccTGCGTCAGCCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCTAAAGTTTAGCCACAGCAGGGTCATTa                                                               <  2:3163401/93‑1 (MQ=255)

                           cTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTaa                                                              <  1:614379/93‑1 (MQ=255)

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

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

                               gTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCCCAGCAGGGCCATTAACCAt                                                          >  1:5682431/1‑93 (MQ=255)

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

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

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

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

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

                                   ggCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATgcgc                                                      >  2:7786088/1‑93 (MQ=255)

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

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

                                      cttagtagCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGc                                                   <  2:1256981/91‑1 (MQ=255)

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

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

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

                                      cGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGc                                                   >  1:3270416/1‑93 (MQ=255)

                                        tagtagCCAGTAACGAAAAACCCCCGCTGCACGATGCGGTTAACCATCAAGACAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCg                                                 <  1:2448953/93‑1 (MQ=255)

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

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

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

                                               cAGTAACGAAAATCACCCGCTGCACGAGGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGaaaa                                          <  1:683622/93‑1 (MQ=255)

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

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

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

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

                                                   aaCGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATACGCCGCCGCTGAAAATAGa                                      >  1:6901913/1‑93 (MQ=255)

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

                                                      gAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAg                                   <  1:1571382/93‑1 (MQ=255)

                                                      gAAAATCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAg                                   <  1:2090037/93‑1 (MQ=255)

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

                                                       aaaaTCACCCGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAgg                                  <  2:2214303/93‑1 (MQ=255)

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

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

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

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

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

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

                                                              cccGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGcc                           <  2:1712319/93‑1 (MQ=255)

                                                              cccGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGcc                           <  1:1665134/93‑1 (MQ=255)

                                                                cGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGACAg                         <  1:7861754/93‑1 (MQ=255)

                                                                cGCTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGAGTCATTATCCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAg                         <  1:3697401/93‑1 (MQ=255)

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

                                                                  cTGCACGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGc                       >  2:2562401/1‑93 (MQ=255)

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

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

                                                                    gCACGTTGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCAt                     <  2:1954521/93‑1 (MQ=255)

                                                                     gaCGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATc                    <  2:5598143/92‑1 (MQ=255)

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

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

                                                                      aCGATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCg                   <  1:5364923/93‑1 (MQ=255)

                                                                        gATGCGGTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCggg                 <  2:5798837/93‑1 (MQ=255)

                                                                         aTGCGGTTAACCATCAAGCCAAAGTTTAGCCATAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGc                >  2:4030743/1‑93 (MQ=255)

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

                                                                             ggTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAGTAGAAAGGAAATGCCAGGCATCGGGCcaca            >  2:6606878/1‑93 (MQ=255)

                                                                             ggTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCcaca            <  2:3836673/93‑1 (MQ=255)

                                                                             ggTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCcaca            <  2:3469171/93‑1 (MQ=255)

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

                                                                              gTTAACCATCAAGCCAAAGTTTAGCCACAGCAGGTTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACaa           <  1:6013615/93‑1 (MQ=255)

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

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

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

                                                                                   ccATCAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACttt      >  2:189718/1‑93 (MQ=255)

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

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

                                                                                       cAAGCCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATa  >  1:6583318/1‑93 (MQ=255)

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

GAAAAGACGCAGGACGGAAAGCAGCCCCTGCGTCAGGCCGTAGTAGCCAGTAACGAAAATCACCCGCTGCACGATGCGGTTAACCATCAAACCAAAGTTTAGCCACAGCAGGGTCATTAACCATGCGCTGCCGCTGAAAATAGAAAGGAAATGCCAGGCATCGGGCCACAAACTTTCATA  >  NC_000913/591076‑591255

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: `‑`