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	82.5%	534.4 / 105.7	251	noncoding (29/76 nt)	glyV	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base A (26/18); new base G (110/97); total (136/115)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

ATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATA  >  NC_000913/4389876‑4390058

ataATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGAcc                                                                                            <  1:3368215/93‑1 (MQ=255)

 taATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCt                                                                                           <  2:7316276/93‑1 (MQ=255)

 taATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCt                                                                                           >  1:6364670/1‑93 (MQ=255)

 taATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCt                                                                                           <  2:1766286/93‑1 (MQ=255)

 taATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCt                                                                                           <  2:2029221/93‑1 (MQ=255)

  aaTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCtt                                                                                          <  1:4618266/93‑1 (MQ=255)

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

     gcgcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGcc                                                                                       >  2:59014/1‑93 (MQ=255)

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

       gcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCaa                                                                                     <  1:3319606/93‑1 (MQ=255)

       gcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCaa                                                                                     <  1:7334254/93‑1 (MQ=255)

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

       gcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCaa                                                                                     <  2:470082/93‑1 (MQ=255)

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

        ccTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAg                                                                                    <  1:6223471/93‑1 (MQ=255)

         cTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAgg                                                                                   >  1:303997/1‑93 (MQ=255)

         cTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAgg                                                                                   >  1:5240436/1‑93 (MQ=255)

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

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

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

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

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

              gTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCggg                                                                              <  2:1342706/93‑1 (MQ=255)

              gTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCggg                                                                              >  2:580456/1‑93 (MQ=255)

              gTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCggg                                                                              >  1:5316057/1‑93 (MQ=255)

               tAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCgggg                                                                             >  2:1446951/1‑93 (MQ=255)

               tAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCgggg                                                                             <  1:1115590/93‑1 (MQ=255)

               tAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                             >  1:3935032/1‑93 (MQ=255)

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

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

                   gacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTcgc                                                                         <  2:5930916/93‑1 (MQ=255)

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

                    acgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTcgcg                                                                        >  2:225581/1‑93 (MQ=255)

                    acgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgcg                                                                        >  2:2422234/1‑93 (MQ=255)

                      gCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAg                                                                      <  1:1310345/93‑1 (MQ=255)

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

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

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

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

                        aGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGtt                                                                    <  1:2197232/93‑1 (MQ=255)

                        aGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGtt                                                                    <  1:7520102/93‑1 (MQ=255)

                        aGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGtt                                                                    <  2:5220682/93‑1 (MQ=255)

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

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

                         gAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                   <  1:354618/93‑1 (MQ=255)

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

                           aaTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                 <  1:6966701/93‑1 (MQ=255)

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

                             tGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                               <  2:4466043/93‑1 (MQ=255)

                              gCGTAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTTGCGAGTTCGAGtc                                                              <  2:4890802/93‑1 (MQ=21)

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

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

                               cGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtct                                                             >  2:1650261/1‑93 (MQ=255)

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

                                gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                            <  1:1915175/93‑1 (MQ=255)

                                gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                            <  1:4147538/93‑1 (MQ=255)

                                 aaaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCg                                                           >  1:669652/1‑93 (MQ=255)

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

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

                                     ttACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                       <  2:1871573/93‑1 (MQ=255)

                                     ttACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                       >  1:5948201/1‑93 (MQ=255)

                                          aaaGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                  >  1:2035331/1‑93 (MQ=255)

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

                                            aGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTcc                                                <  1:2485431/93‑1 (MQ=255)

                                              cAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaa                                              <  2:6712976/93‑1 (MQ=255)

                                                aaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaaa                                            >  1:7475833/1‑93 (MQ=255)

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

                                                aaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaaa                                            >  1:5771355/1‑93 (MQ=255)

                                                aaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaaa                                            <  1:3661617/93‑1 (MQ=255)

                                                 aaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAt                                           <  2:6794344/93‑1 (MQ=255)

                                                 aaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAt                                           >  2:6588844/1‑93 (MQ=255)

                                                  aTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAtt                                          >  1:38963/1‑93 (MQ=255)

                                                  aTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAtt                                          >  2:1557890/1‑93 (MQ=255)

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

                                                    tAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTg                                        <  1:6775926/93‑1 (MQ=255)

                                                    tAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTg                                        >  1:7393518/1‑93 (MQ=255)

                                                       gtagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaa                                     >  2:4629408/1‑93 (MQ=255)

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

                                                        tagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaaa                                    <  2:6348370/93‑1 (MQ=38)

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

                                                           taCGCGGGAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc                                                              >  2:5366077/1‑64 (MQ=255)

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

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

                                                                 ggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaa                           >  1:3426757/1‑93 (MQ=40)

                                                                  gAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaaa                          >  2:6125952/1‑93 (MQ=255)

                                                                  gAATAGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaaa                          <  1:4893032/93‑1 (MQ=255)

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

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

                                                                   aaTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAg                         <  2:5609899/93‑1 (MQ=255)

                                                                      aGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcac                      <  2:3672907/93‑1 (MQ=255)

                                                                      aGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcac                      <  2:3040428/93‑1 (MQ=255)

                                                                      aGCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcac                      <  2:6893642/93‑1 (MQ=255)

                                                                       gCTCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcaca                     <  1:5067246/93‑1 (MQ=255)

                                                                         tCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGa                   >  2:2246607/1‑93 (MQ=255)

                                                                         tCAGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGa                   >  1:1869189/1‑93 (MQ=255)

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

                                                                           aGTTGGTAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGacc                 <  1:1691276/93‑1 (MQ=255)

                                                                                 tAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCaa           >  2:7111495/1‑93 (MQ=255)

                                                                                 tAGAGCACGGCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCaa           >  2:2887125/1‑93 (MQ=255)

                                                                                          aCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATa  <  1:1635283/93‑1 (MQ=255)

ATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATA  >  NC_000913/4389876‑4390058

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