BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,075,220	(T)_5→6	coding (256/357 nt)	yggP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,075,215	1	.	T	97.0%	577.6 / 7.5	168	coding (261/357 nt)	yggP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/4); new base T (79/84); total (80/88)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

TAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑TTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAA  >  NC_000913/3075124‑3075302

taAGTGGCATTTCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑t                                                                                         >  2:900617/1‑93 (MQ=255)

 aaGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tt                                                                                        >  2:4388442/1‑93 (MQ=255)

 aaGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tt                                                                                        <  1:2619576/93‑1 (MQ=255)

 aaGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tt                                                                                        <  1:1871843/93‑1 (MQ=255)

  aGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑ttt                                                                                       >  2:575302/1‑93 (MQ=255)

   gTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt                                                                                      >  1:1923163/1‑93 (MQ=255)

   gTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt                                                                                      >  2:4512758/1‑93 (MQ=255)

   gTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt                                                                                      >  2:1482081/1‑93 (MQ=255)

   gTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt                                                                                      >  1:3041186/1‑93 (MQ=255)

    tGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt‑t                                                                                    >  1:1670157/1‑93 (MQ=255)

    tGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt‑t                                                                                    <  1:4273736/93‑1 (MQ=255)

     ggCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt‑tt                                                                                   >  1:3696680/1‑92 (MQ=255)

      gCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt‑tta                                                                                  <  2:4031403/93‑3 (MQ=255)

        aTTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt‑ttata                                                                                <  1:2799970/93‑5 (MQ=255)

        aTTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt‑ttata                                                                                >  1:933667/1‑89 (MQ=255)

        aTTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑tttt‑ttata                                                                                >  2:98178/1‑89 (MQ=255)

          tGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATAtt                                                                               <  2:3668653/93‑1 (MQ=255)

           gccgccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATAttt                                                                              <  1:3721721/93‑1 (MQ=255)

            ccgccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTg                                                                             <  1:1301126/93‑1 (MQ=255)

             cgccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGc                                                                            <  1:2562554/93‑1 (MQ=255)

             caccAAAAGATTCGTAGCTATAGCCCGCTGAATATCCGCTTAATACCCCGATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGc                                                                            <  2:386547/91‑1 (MQ=255)

              cccAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGcc                                                                           <  2:3032146/92‑1 (MQ=255)

                cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAg                                                                         >  2:4446651/1‑93 (MQ=255)

                cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAg                                                                         >  2:4732258/1‑93 (MQ=255)

                cAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAg                                                                         >  1:397384/1‑93 (MQ=255)

                  aaaTATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTg                                                                       >  2:4179315/1‑93 (MQ=255)

                     tatTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTgaga                                                                    <  1:1246610/93‑1 (MQ=255)

                       ttCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGAtt                                                                  >  1:2298789/1‑93 (MQ=255)

                       ttCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGAtt                                                                  <  1:4749867/93‑1 (MQ=255)

                        tCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGAttt                                                                 <  1:4432663/93‑1 (MQ=255)

                         cGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGAtttt                                                                <  1:3882106/93‑1 (MQ=255)

                          gTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGAttttt                                                               <  1:881421/93‑1 (MQ=255)

                           tAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATATAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTa                                                              >  1:3132052/1‑93 (MQ=255)

                           tAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTa                                                              >  2:1160780/1‑93 (MQ=255)

                             gCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTAcc                                                            <  1:2267641/93‑1 (MQ=255)

                               tGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCa                                                          <  1:3447072/93‑1 (MQ=255)

                               tGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCa                                                          >  2:686755/1‑93 (MQ=255)

                                  aGCACGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACtt                                                       <  2:2682485/93‑1 (MQ=255)

                                   gCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTc                                                      >  2:3089224/1‑93 (MQ=255)

                                   gCCCCCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTc                                                      <  1:1603198/93‑1 (MQ=255)

                                     ccGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGa                                                    >  1:3019348/1‑93 (MQ=255)

                                      cGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGAc                                                   <  1:4261762/93‑1 (MQ=255)

                                      cGCTGAATATCAGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGCTTTTTACCCACTTCGAc                                                   >  2:622203/1‑93 (MQ=255)

                                       gCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACa                                                  >  2:1325099/1‑93 (MQ=255)

                                       gCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACa                                                  >  2:1208396/1‑93 (MQ=255)

                                         tGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAAt                                                <  2:2846364/93‑1 (MQ=255)

                                          gAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATa                                               >  2:2700478/1‑93 (MQ=255)

                                          gAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATa                                               <  2:1077435/93‑1 (MQ=255)

                                           aaTATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATaa                                              <  2:3981550/93‑1 (MQ=255)

                                            atatCCGCTTGGTACCCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAc                                             <  1:1520507/93‑1 (MQ=255)

                                            atatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAc                                             >  2:2355933/1‑93 (MQ=255)

                                             tatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcc                                            <  1:2335964/93‑1 (MQ=255)

                                             tatCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcc                                            <  1:3469664/93‑1 (MQ=255)

                                              ttCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAccc                                           <  1:3605125/92‑1 (MQ=255)

                                              atCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAccc                                           <  1:2812710/93‑1 (MQ=255)

                                               tCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcccc                                          >  2:36098/1‑93 (MQ=255)

                                               tCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcccc                                          <  1:2687073/93‑1 (MQ=255)

                                               tCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAAcccc                                          <  1:3091886/93‑1 (MQ=255)

                                                ccGCTTGGTAACCCCATCGCCGGTTGCAATACNAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCg                                         <  1:1913638/93‑1 (MQ=255)

                                                   cTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGca                                      <  1:1861573/93‑1 (MQ=255)

                                                    ttGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGcac                                     <  2:1863146/93‑1 (MQ=255)

                                                      ggTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACAt                                   <  2:2197300/93‑1 (MQ=255)

                                                       gTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACAtt                                  <  2:3376804/93‑1 (MQ=255)

                                                        tAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTc                                 >  1:1553152/1‑93 (MQ=255)

                                                        tAACCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTc                                 >  2:2497468/1‑93 (MQ=255)

                                                         aaCCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCa                                <  1:3416030/93‑1 (MQ=255)

                                                          aCCCCATCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCAt                               <  2:4167759/93‑1 (MQ=255)

                                                               aTCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCg                          >  1:3176740/1‑93 (MQ=255)

                                                               aTCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCg                          >  2:2608679/1‑93 (MQ=255)

                                                                tCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCgg                         >  2:778090/1‑93 (MQ=255)

                                                                tCGCCGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCgg                         <  1:2976465/93‑1 (MQ=255)

                                                                 cGCCGGTGGCAATACAAAACGCTGGCC‑TTTTGTATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGt                        <  1:2259273/93‑1 (MQ=255)

                                                                   ccGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTaa                      <  1:4170820/93‑1 (MQ=255)

                                                                   ccGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTaa                      >  1:1792982/1‑93 (MQ=255)

                                                                    cGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAAt                     <  2:3130032/93‑1 (MQ=255)

                                                                    cGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAAt                     >  1:1471402/1‑93 (MQ=255)

                                                                    cGGTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAAt                     <  2:364177/93‑1 (MQ=255)

                                                                     ggTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg                    <  2:255307/93‑1 (MQ=255)

                                                                     ggTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATg                    >  2:1905117/1‑93 (MQ=255)

                                                                      gTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa                   >  1:2388599/1‑93 (MQ=255)

                                                                      gTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa                   <  2:2586472/93‑1 (MQ=255)

                                                                      gTTGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGa                   <  1:1122085/93‑1 (MQ=255)

                                                                        tGCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGAcc                 <  1:1063433/93‑1 (MQ=255)

                                                                         gCAATACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGCAATGACCg                >  1:3375521/1‑93 (MQ=255)

                                                                           aaTACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGga              >  2:92212/1‑93 (MQ=255)

                                                                             tACAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGgatg            >  2:3418934/1‑93 (MQ=255)

                                                                              aCAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCACGGCACATTCATGCCCGGTAATGACCGgatga           >  1:4591816/1‑93 (MQ=255)

                                                                               cAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGgatgat          >  2:3114930/1‑93 (MQ=255)

                                                                               cAAAACGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGgatgat          <  2:2444501/93‑1 (MQ=255)

                                                                                aaaaCGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtt         <  2:4339043/93‑1 (MQ=255)

                                                                                aaaaCGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtt         >  2:4240446/1‑93 (MQ=255)

                                                                                aaaaCGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtt         >  2:1618777/1‑93 (MQ=255)

                                                                                aaaaCGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtt         >  1:3813921/1‑93 (MQ=255)

                                                                                  aaCGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtttt       <  2:3860415/93‑1 (MQ=255)

                                                                                  aaCGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGAtttt       >  1:62398/1‑93 (MQ=255)

                                                                                   aCGCTGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTc      >  2:2994698/1‑93 (MQ=255)

                                                                                       tGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTaaa  <  2:3253856/93‑1 (MQ=255)

                                                                                       tGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTaaa  >  2:4640588/1‑93 (MQ=255)

                                                                                       tGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTaaa  >  2:2236733/1‑93 (MQ=255)

                                                                                       tGGCCTTTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTaaa  <  1:4378006/93‑1 (MQ=255)

TAAGTGGCATTGCCGCCAAAATATTCGTAGCTGTAGCCCGCTGAATATCCGCTTGGTAACCCCATCGCCGGTTGCAATACAAAACGCTGGCC‑TTTT‑TATATTTGCCAGTGAGATTTTTACCCACTTCGACAATAACCCCGGCACATTCATGCCCGGTAATGACCGGATGATTTTCTAAA  >  NC_000913/3075124‑3075302

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