BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,389,966	A→G	noncoding (29/76 nt)	glyV →	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,389,966	0	A	G	99.1%	1053.2 / ‑1.6	324	noncoding (29/76 nt)	glyV	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (179/141); total (182/142)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 333 total aligned reads displayed.

ATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGA  >  NC_000913/4389876‑4390055

  aaTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCtt                                                                                       <  1:5842929/93‑1 (MQ=255)

  aaTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCtt                                                                                       >  2:1464750/1‑93 (MQ=255)

  aaTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCtt                                                                                       >  1:2413366/1‑93 (MQ=255)

   aTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTg                                                                                      >  2:4042881/1‑93 (MQ=255)

    tGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGc                                                                                     >  2:5187315/1‑93 (MQ=255)

      cgcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCa                                                                                   <  2:1633616/93‑1 (MQ=255)

      cgcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCa                                                                                   <  1:1730958/93‑1 (MQ=255)

       gcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCaa                                                                                  >  2:3429943/1‑93 (MQ=255)

          tCCCGTAACGACGCAGAAATGCGAAAATTACGACAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGACAAGGt                                                                               >  1:3875638/1‑93 (MQ=255)

          tCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAggc                                                                               >  1:223379/1‑92 (MQ=255)

          tCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGt                                                                               >  2:5359017/1‑93 (MQ=255)

          tCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGt                                                                               <  2:2480229/93‑1 (MQ=255)

           cccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTc                                                                              >  2:734446/1‑93 (MQ=255)

           cccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTc                                                                              <  1:4461098/93‑1 (MQ=255)

            ccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCg                                                                             <  2:2095215/93‑1 (MQ=255)

            ccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCg                                                                             >  2:6675427/1‑93 (MQ=255)

            ccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCg                                                                             <  1:7925724/93‑1 (MQ=255)

            ccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCg                                                                             <  1:7062397/93‑1 (MQ=255)

             cGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGTCTTGCCAAGGTggg                                                                            >  2:296720/1‑93 (MQ=255)

                aaCGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGt                                                                         >  2:4074366/1‑93 (MQ=255)

                aaCGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGt                                                                         <  1:2633845/93‑1 (MQ=255)

                aaCGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGt                                                                         >  2:6207754/1‑93 (MQ=255)

                 acgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTc                                                                        >  1:7407642/1‑93 (MQ=255)

                   gacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTcgc                                                                      >  1:3492011/1‑93 (MQ=255)

                      gCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACTCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAg                                                                   >  1:1134433/1‑93 (MQ=255)

                      gCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAg                                                                   >  1:7578623/1‑93 (MQ=255)

                       cAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGt                                                                  >  2:4401177/1‑93 (MQ=255)

                         gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                <  1:6615184/93‑1 (MQ=255)

                         gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                <  2:1068616/93‑1 (MQ=255)

                         gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                >  1:1468754/1‑93 (MQ=255)

                         gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                <  2:630425/93‑1 (MQ=255)

                          aaaTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                               >  1:4969834/1‑93 (MQ=255)

                          aaaTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTAGCGAGTTCg                                                               >  1:1381976/1‑93 (MQ=255)

                            aTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                             >  2:1770906/1‑93 (MQ=255)

                            aTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                             <  2:641043/93‑1 (MQ=255)

                              gCGAAAATTACGAAAGCAAACTTAAGTAGCAAGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc                                                           <  2:7555702/93‑1 (MQ=17)

                              gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc                                                           >  1:5685604/1‑93 (MQ=255)

                              gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc                                                           <  2:1896732/93‑1 (MQ=255)

                               cGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtct                                                          <  2:7635044/93‑1 (MQ=255)

                                gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                         >  2:5863975/1‑93 (MQ=255)

                                gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                         <  2:1246676/93‑1 (MQ=255)

                                gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                         >  1:7081341/1‑93 (MQ=255)

                                  aaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt                                                       >  1:4227931/1‑93 (MQ=255)

                                  aaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt                                                       <  1:3319666/93‑1 (MQ=255)

                                  aaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt                                                       <  2:7030051/93‑1 (MQ=255)

                                   aagtACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGTCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                                      <  2:3224642/90‑1 (MQ=255)

                                   aaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGTGTTCGAGTCTCGtt                                                      <  1:8669017/93‑1 (MQ=255)

                                   aaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                                      <  2:7950689/93‑1 (MQ=255)

                                    aTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGttt                                                     >  2:3737461/1‑93 (MQ=255)

                                       aCGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTccc                                                  <  1:969413/93‑1 (MQ=255)

                                           aaGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                              <  1:5558933/93‑1 (MQ=255)

                                           aaGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                              >  2:194660/1‑93 (MQ=255)

                                             gCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCa                                            >  2:2367626/1‑93 (MQ=255)

                                               aaaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaa                                          >  2:8255934/1‑93 (MQ=255)

                                               aaaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaa                                          >  1:3700119/1‑93 (MQ=255)

                                                aaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaaa                                         <  2:6481255/93‑1 (MQ=255)

                                                   ttAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAttt                                      <  2:3380344/93‑1 (MQ=255)

                                                    tAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTg                                     >  1:8590462/1‑93 (MQ=255)

                                                     aaGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGa                                    <  1:7708713/93‑1 (MQ=255)

                                                      agtagtACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaa                                   <  1:6877272/93‑1 (MQ=255)

                                                      agtagtACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaa                                   >  1:3741056/1‑93 (MQ=255)

                                                        tagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaaa                                 <  1:6105286/93‑1 (MQ=255)

                                                        tagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaaa                                 <  1:5069047/93‑1 (MQ=255)

                                                        tagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaaa                                 <  1:7130150/93‑1 (MQ=255)

                                                        gagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaaa                                 <  1:4813461/92‑1 (MQ=255)

                                                          gtaCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAga                               >  2:3522618/1‑92 (MQ=255)

                                                          gtaCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGt                               <  1:1921836/93‑1 (MQ=255)

                                                           taCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtg                              <  1:676084/93‑1 (MQ=255)

                                                           taCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtg                              <  2:6928540/93‑1 (MQ=255)

                                                            aCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtgc                             <  1:8416692/93‑1 (MQ=255)

                                                             cgcgGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtgct                            <  2:8289963/93‑1 (MQ=255)

                                                               cgGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaaa                                 >  2:4363172/1‑86 (MQ=14)

                                                               cgGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtgctgc                          >  2:5574502/1‑93 (MQ=255)

                                                                gggAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCa                         <  1:550560/93‑1 (MQ=255)

                                                                gggAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCa                         >  1:4704872/1‑93 (MQ=255)

                                                                gggAATAGCTCAGTGGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTGCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCa                         <  2:876334/93‑1 (MQ=255)

                                                                 ggAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaa                        <  2:6855848/93‑1 (MQ=255)

                                                                  gAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGCCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaag                       >  1:7496344/1‑92 (MQ=255)

                                                                   aaTAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAg                      >  2:8399326/1‑93 (MQ=255)

                                                                   aaTAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAACTACTGCAAAg                      >  2:1727189/1‑93 (MQ=255)

                                                                    aTAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAACTGCTGCatata                     >  2:8736923/1‑89 (MQ=11)

                                                                     tAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGca                    >  2:2705143/1‑93 (MQ=255)

                                                                     tAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGca                    <  1:8317568/93‑1 (MQ=255)

                                                                        cTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAg                 <  1:1021774/93‑1 (MQ=255)

                                                                          cAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGac               <  1:7986537/93‑1 (MQ=255)

                                                                          cAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGac               >  2:3839360/1‑93 (MQ=255)

                                                                          cAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAACTGCTGCAAAGCACAAac               >  2:7217142/1‑93 (MQ=14)

                                                                           aGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGCAAAGTGCTGCAAAGCACAGacc              >  2:3036530/1‑93 (MQ=255)

                                                                             ttGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCACAATCTGTAAAATGCTGAAAAGCACAGacacc            >  1:8047823/1‑90 (MQ=11)

                                                                               ggTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCAccc          >  1:2877025/1‑93 (MQ=255)

                                                                               ggTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCAccc          >  1:3358860/1‑93 (MQ=255)

                                                                                 tAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCaa        <  2:4201549/93‑1 (MQ=255)

                                                                                 tAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAAACCACCCaa        >  2:6757818/1‑93 (MQ=255)

                                                                                   aggCACGGCCGTGGCAAAGTCGGGGTCGCGACTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCGAGc      <  1:6324736/91‑1 (MQ=14)

                                                                                    agCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCg     <  2:5721879/93‑1 (MQ=255)

                                                                                     gCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCgg    >  1:8150676/1‑93 (MQ=255)

                                                                                     gCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGACACGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCgg    <  2:7808600/93‑1 (MQ=255)

                                                                                       aCGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAATGCACAGACCACCCAAGCGGGa  <  2:8676187/93‑1 (MQ=255)

                                                                                       aCGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa  >  2:8073367/1‑93 (MQ=255)

                                                                                       aCGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGa  >  2:4607328/1‑93 (MQ=255)

ATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGA  >  NC_000913/4389876‑4390055

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