BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,553	A→T	K286* (AAA→TAA)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,553	0	A	T	98.3%	401.7 / ‑2.6	117	K286* (AAA→TAA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand): ref base A (1/1); new base T (50/65); total (51/66)
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.76e-01

Only 100 of 120 total aligned reads displayed.

GAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAG  >  NC_000913/1756457‑1756649

gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaa                                                                                               <  1:2017950/100‑1 (MQ=255)

  agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAc                                                                                             >  2:556470/1‑100 (MQ=255)

    aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAccc                                                                                           >  2:2393370/1‑100 (MQ=255)

      gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCAc                                                                                         >  1:2416740/1‑100 (MQ=255)

          tCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaagatcggaaga                                                                                     >  2:1512508/1‑90 (MQ=255)

          tCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGccc                                                                                     <  2:2381027/100‑1 (MQ=255)

          tCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGccc                                                                                     <  1:1443535/100‑1 (MQ=255)

           cTTCGCCCAGAAGATGATGATCGAAAAATGTCTCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCg                                                                                    >  2:365281/1‑100 (MQ=255)

           cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCg                                                                                    <  1:505812/100‑1 (MQ=255)

           cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCg                                                                                    >  1:780080/1‑100 (MQ=255)

             tCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGAc                                                                                  <  1:819869/100‑1 (MQ=255)

              cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCATTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACt                                                                                 <  2:1499292/100‑1 (MQ=255)

               gCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTc                                                                                >  2:1428049/1‑100 (MQ=255)

                cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcg                                                                               <  1:2077439/100‑1 (MQ=255)

                cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcg                                                                               <  2:2564527/100‑1 (MQ=255)

                cccAGAAGATGATGATCGAAAAATGTATCCATGGACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcg                                                                               >  2:1768345/1‑100 (MQ=255)

                 ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgc                                                                              >  2:463131/1‑100 (MQ=255)

                  cAGAAGATGATTATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgcn                                                                             <  2:2548193/100‑2 (MQ=255)

                  cAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgcg                                                                             >  2:1797287/1‑100 (MQ=255)

                  cAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcacg                                                                             <  1:1083254/100‑1 (MQ=255)

                   agaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCgcgc                                                                            <  2:1960269/100‑1 (MQ=255)

                   agaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCgcgc                                                                            <  1:2585028/100‑1 (MQ=255)

                         tgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGc                                                                      <  2:2503546/100‑1 (MQ=255)

                         tgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGc                                                                      <  2:962071/100‑1 (MQ=255)

                          gatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGcc                                                                     >  1:1386128/1‑100 (MQ=255)

                           atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                                    >  2:334595/1‑100 (MQ=255)

                           atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                                    >  2:561896/1‑100 (MQ=255)

                            tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                                   >  2:1695222/1‑100 (MQ=255)

                            tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                                   <  1:1229643/100‑1 (MQ=255)

                             gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                                                  >  1:2234526/1‑100 (MQ=255)

                             gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                                                  <  1:1355901/100‑1 (MQ=255)

                             gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                                                  >  1:28521/1‑100 (MQ=255)

                              atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                                 >  2:1263307/1‑100 (MQ=255)

                              atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                                 >  2:132203/1‑100 (MQ=255)

                              atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                                 >  1:125364/1‑100 (MQ=255)

                              atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                                 >  2:46458/1‑100 (MQ=255)

                                 gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg                                                              >  2:1680987/1‑100 (MQ=255)

                                   aaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGtt                                                            <  2:1314346/100‑1 (MQ=255)

                                    aaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                                           <  1:1327209/100‑1 (MQ=255)

                                    aaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                                           <  2:2562853/100‑1 (MQ=255)

                                     aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                                          >  2:2441029/1‑100 (MQ=255)

                                      aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCa                                                         <  2:118946/100‑1 (MQ=255)

                                      aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCa                                                         <  1:393720/100‑1 (MQ=255)

                                         tATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAAc                                                      >  1:1363717/1‑100 (MQ=255)

                                          aTCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACg                                                     <  1:2393394/100‑1 (MQ=255)

                                           tCCGTGCACGTAAAGTAGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTGAAACACACGCCCGACTCGCGCAGACGCCGGTGACGTTGCAAACGc                                                    <  2:2416716/100‑1 (MQ=255)

                                            ccGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGcc                                                   <  1:1439162/100‑1 (MQ=255)

                                             cGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCa                                                  >  2:2569896/1‑100 (MQ=255)

                                              gTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCAt                                                 <  1:1428063/100‑1 (MQ=255)

                                              gTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCAt                                                 >  2:852102/1‑100 (MQ=255)

                                               tGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATc                                                >  1:2087718/1‑100 (MQ=255)

                                                gCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                                               <  2:304022/100‑1 (MQ=255)

                                                gCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                                               >  1:2648489/1‑100 (MQ=255)

                                                   ngTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                            >  2:1182239/2‑100 (MQ=255)

                                                    gTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGa                                           <  1:1695237/100‑1 (MQ=255)

                                                     tAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                          <  1:836070/100‑1 (MQ=255)

                                                      aaaGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACg                                         >  2:329560/1‑100 (MQ=255)

                                                      aaaGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACg                                         >  2:1408606/1‑100 (MQ=255)

                                                         gTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTa                                      >  1:761402/1‑100 (MQ=255)

                                                           cGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt                                    <  2:426065/100‑1 (MQ=255)

                                                           cGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt                                    <  1:2530893/100‑1 (MQ=255)

                                                            gTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTg                                   <  2:659976/100‑1 (MQ=255)

                                                                tCACTGCGACCCAGATGCTGGATTCCATGATATAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                               <  1:365286/100‑1 (MQ=255)

                                                                 cACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCa                              <  1:733248/100‑1 (MQ=255)

                                                                 cACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCa                              >  1:1613991/1‑100 (MQ=255)

                                                                   cTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                             >  1:207411/1‑99 (MQ=255)

                                                                   cTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt                            <  1:2298165/100‑1 (MQ=255)

                                                                    tGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTg                           <  1:1845511/100‑1 (MQ=255)

                                                                     gCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGa                          <  2:109219/100‑1 (MQ=255)

                                                                       gACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                        <  2:2299732/100‑1 (MQ=255)

                                                                        aCCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc                       <  1:2539259/100‑1 (MQ=255)

                                                                         cccAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCt                      >  2:2492362/1‑100 (MQ=255)

                                                                            aGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTc                   <  2:2178097/100‑1 (MQ=255)

                                                                               tGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg                <  1:1050799/100‑1 (MQ=255)

                                                                                gCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGt               >  1:2285370/1‑100 (MQ=255)

                                                                                gCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGt               <  1:2556807/100‑1 (MQ=255)

                                                                                gCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGt               >  1:1803648/1‑100 (MQ=255)

                                                                                 cTGGATTCCATGATCTAAAGCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg              <  2:2069478/100‑1 (MQ=255)

                                                                                 cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg              <  2:2016299/100‑1 (MQ=255)

                                                                                 cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg              <  2:853266/100‑1 (MQ=255)

                                                                                 cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg              <  1:1797307/100‑1 (MQ=255)

                                                                                 cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg              <  1:986592/100‑1 (MQ=255)

                                                                                   ggATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa            <  1:579442/100‑1 (MQ=255)

                                                                                    gATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAAt           <  1:561903/100‑1 (MQ=255)

                                                                                     aTTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATc          >  1:300511/1‑100 (MQ=255)

                                                                                      ttCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATcc         <  1:1249449/100‑1 (MQ=255)

                                                                                      ttCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATcc         <  2:1363706/100‑1 (MQ=255)

                                                                                         cATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCGATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGc       >  1:1605037/1‑99 (MQ=255)

                                                                                         cATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa      >  2:1785067/1‑100 (MQ=255)

                                                                                         cATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa      <  1:42573/100‑1 (MQ=255)

                                                                                         cATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa      >  2:38613/1‑100 (MQ=255)

                                                                                         cATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa      >  1:136298/1‑100 (MQ=255)

                                                                                         cATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa      <  2:1166491/100‑1 (MQ=255)

                                                                                          aTGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAACCGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaa     >  2:1266610/1‑100 (MQ=255)

                                                                                            gATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaaa   <  1:2047167/100‑1 (MQ=255)

                                                                                            gATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaaa   <  1:1785087/100‑1 (MQ=255)

                                                                                            gATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaaa   <  1:2492387/100‑1 (MQ=255)

                                                                                            gATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGACATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaaa   <  2:1298551/100‑1 (MQ=255)

                                                                                             aTCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAg  >  1:319418/1‑100 (MQ=255)

                                                                                             aTCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAg  <  2:888855/100‑1 (MQ=255)

GAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAG  >  NC_000913/1756457‑1756649

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