BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,154,920	C→T	P192S (CCG→TCG)	tmk →	thymidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,154,920	0	C	T	85.7%	660.3 / 68.4	238	P192S (CCG→TCG)	tmk	thymidylate kinase
Reads supporting (aligned to +/- strand): ref base C (16/18); new base T (97/107); total (113/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

Only 100 of 238 total aligned reads displayed.

GAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACGACCT  >  NC_000913/1154830‑1155011

gAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCg                                                                                           >  1:6431569/1‑93 (MQ=255)

gAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCg                                                                                           >  1:781911/1‑93 (MQ=255)

gAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCg                                                                                           <  1:2197636/93‑1 (MQ=255)

gAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCg                                                                                           >  2:4867083/1‑93 (MQ=255)

 aaTCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGc                                                                                          <  2:677185/93‑1 (MQ=255)

  aTCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCt                                                                                         >  2:7469381/1‑93 (MQ=255)

       tCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAgg                                                                                    <  1:126095/93‑1 (MQ=255)

         gATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGcc                                                                                  <  2:2472424/93‑1 (MQ=255)

          aTTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCg                                                                                 <  2:5543922/93‑1 (MQ=255)

           tttCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGt                                                                                >  2:7397151/1‑93 (MQ=255)

           tttCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGt                                                                                >  1:6064939/1‑93 (MQ=255)

           tttCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGt                                                                                >  2:6557700/1‑93 (MQ=255)

           tttCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGt                                                                                >  2:6362451/1‑93 (MQ=255)

             tCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTga                                                                              <  1:5825271/93‑1 (MQ=255)

             tCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTga                                                                              <  1:1725647/93‑1 (MQ=255)

             tCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTga                                                                              >  2:570048/1‑93 (MQ=255)

              cTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTgat                                                                             <  2:1362687/93‑1 (MQ=255)

              cTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTgat                                                                             >  1:4188968/1‑93 (MQ=255)

               tttAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTgatg                                                                            >  2:6031380/1‑93 (MQ=255)

                ttAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTgatgg                                                                           <  1:2048627/93‑1 (MQ=255)

                   aTCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTgatggatg                                                                        >  1:3566190/1‑93 (MQ=255)

                    tCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGc                                                                       >  2:5442395/1‑93 (MQ=255)

                     cGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCa                                                                      <  1:3066919/93‑1 (MQ=255)

                      gCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCaa                                                                     <  1:1587006/93‑1 (MQ=255)

                        aCCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATc                                                                   <  2:4617937/93‑1 (MQ=255)

                                ccGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAc                                                           <  2:5815903/93‑1 (MQ=255)

                                ccGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAc                                                           <  2:2136555/93‑1 (MQ=255)

                                ccGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAc                                                           <  2:4939439/93‑1 (MQ=255)

                                 cGCTGTCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAcc                                                          >  2:7372424/1‑93 (MQ=255)

                                 cGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAcc                                                          >  2:5932897/1‑93 (MQ=255)

                                 cGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAcc                                                          >  1:3651388/1‑93 (MQ=255)

                                 cGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTAcc                                                          >  2:556643/1‑93 (MQ=255)

                                 cGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTAcc                                                          <  1:6536945/93‑1 (MQ=255)

                                  gCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCg                                                         >  2:3499263/1‑93 (MQ=255)

                                  gCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCg                                                         >  1:1124042/1‑93 (MQ=255)

                                  gCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCg                                                         >  1:2836624/1‑93 (MQ=255)

                                    tATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTg                                                       >  1:4856896/1‑93 (MQ=255)

                                     aTCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGa                                                      >  1:7121799/1‑93 (MQ=255)

                                        tGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGAccc                                                   <  2:6740352/93‑1 (MQ=255)

                                         ggAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCAGTACCGTGACCCa                                                  <  1:5324512/93‑1 (MQ=255)

                                          gAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCAc                                                 <  2:1179218/93‑1 (MQ=255)

                                          gAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCAc                                                 <  1:3314594/93‑1 (MQ=255)

                                           aaCTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACt                                                <  2:2030225/93‑1 (MQ=255)

                                             cTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTgg                                              <  1:4357135/93‑1 (MQ=255)

                                             cTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTgg                                              <  1:6266154/93‑1 (MQ=255)

                                              tGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTggg                                             >  2:1189417/1‑93 (MQ=255)

                                               ggCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGt                                            <  2:7047070/93‑1 (MQ=255)

                                                  agcaCAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGaa                                         >  1:1000370/1‑93 (MQ=255)

                                                   gcaCAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAg                                        <  1:2815620/93‑1 (MQ=255)

                                                   gcaCAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAg                                        >  1:4429860/1‑93 (MQ=255)

                                                   gcaCAAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAg                                        <  1:4603801/93‑1 (MQ=255)

                                                    cacaAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAgg                                       <  2:5474705/93‑1 (MQ=255)

                                                     acaAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGa                                      >  1:7447444/1‑93 (MQ=255)

                                                     acaAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGa                                      <  2:2528716/93‑1 (MQ=255)

                                                     acaAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGa                                      <  2:3952360/93‑1 (MQ=255)

                                                     acaAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGa                                      <  2:3692030/93‑1 (MQ=255)

                                                      caAGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAg                                     >  2:872838/1‑93 (MQ=255)

                                                       aaGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGt                                    <  1:5172455/93‑1 (MQ=255)

                                                        aGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGtt                                   >  2:7106282/1‑93 (MQ=255)

                                                        aGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGtt                                   >  1:7017837/1‑93 (MQ=255)

                                                        aGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGtt                                   <  2:299609/93‑1 (MQ=255)

                                                        aGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGtt                                   >  2:3843157/1‑93 (MQ=255)

                                                        aGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGtt                                   >  2:3228300/1‑93 (MQ=255)

                                                        aGATAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGtt                                   <  1:5412099/93‑1 (MQ=255)

                                                         gATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTg                                  <  2:4073589/93‑1 (MQ=255)

                                                          aTAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTgg                                 <  1:6983890/93‑1 (MQ=255)

                                                           tAAAAGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGa                                <  2:2948495/93‑1 (MQ=255)

                                                            aaaaGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGAc                               <  2:534028/93‑1 (MQ=255)

                                                            aaaaGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGAc                               >  2:1538134/1‑93 (MQ=255)

                                                            aaaaGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGAc                               <  1:2403272/93‑1 (MQ=255)

                                                              aaGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGc                             >  2:4355798/1‑93 (MQ=255)

                                                               aGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCa                            <  1:6730197/93‑1 (MQ=255)

                                                               aGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCa                            <  1:4099053/93‑1 (MQ=255)

                                                               aGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCa                            >  2:3067541/1‑93 (MQ=255)

                                                               aGCATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCa                            >  1:844417/1‑93 (MQ=255)

                                                                 cATTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATg                          <  1:3144223/93‑1 (MQ=255)

                                                                  aTTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATga                         <  1:5475322/93‑1 (MQ=255)

                                                                  aTTCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATga                         >  2:1435172/1‑93 (MQ=255)

                                                                  aTTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATga                         <  1:3722194/93‑1 (MQ=255)

                                                                    tCATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATgaga                       >  2:26211/1‑93 (MQ=255)

                                                                     cATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGAt                      <  1:505765/93‑1 (MQ=255)

                                                                     cATACCATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGAt                      <  2:6409700/93‑1 (MQ=255)

                                                                         nnaTTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTa                  <  2:88465/91‑1 (MQ=255)

                                                                         ccATTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTa                  <  2:5852328/93‑1 (MQ=255)

                                                                          cATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTAt                 >  1:1443826/1‑93 (MQ=255)

                                                                           aTTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATc                <  1:3612345/93‑1 (MQ=255)

                                                                           aTTGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATc                <  2:5748155/93‑1 (MQ=255)

                                                                            ttGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGCCGCATGAGATGGTATcc               >  2:6711683/1‑93 (MQ=255)

                                                                             tGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCa              <  1:2644236/93‑1 (MQ=255)

                                                                             tGATGCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCa              <  2:5762478/93‑1 (MQ=255)

                                                                                 gCCACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGt          <  1:1893300/93‑1 (MQ=255)

                                                                                  ccACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGtt         <  2:6123189/93‑1 (MQ=255)

                                                                                  ccACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGtt         <  1:2458246/93‑1 (MQ=255)

                                                                                   cACCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTa        <  1:5605426/93‑1 (MQ=255)

                                                                                    aCCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTAc       <  1:4391277/93‑1 (MQ=255)

                                                                                    aCCCAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTAc       <  1:3888220/93‑1 (MQ=255)

                                                                                      ccAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACGa     >  2:1787239/1‑93 (MQ=255)

                                                                                      ccAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACGa     >  1:930426/1‑93 (MQ=255)

                                                                                       cAGTCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACGAc    >  2:239333/1‑93 (MQ=255)

                                                                                        ggCCGCGGGAGGCCGTGATGAAGGCAACCCGCACTACGGTGACCCACTGGGTG‑AGGAGTTGGACGCATGAGATGGTATCCATGGTTACGACCt  <  2:267847/92‑1 (MQ=255)

GAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATGCAATCCGCACTACCGTGACCCACTGGGTGAAGGAGTTGGACGCATGAGATGGTATCCATGGTTACGACCT  >  NC_000913/1154830‑1155011

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