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	99.1%	407.3 / NA	116	K286* (AAA→TAA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand): ref base A (0/0); major base T (61/54); minor base C (0/1); total (61/55)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAG  >  NC_000913/1756457‑1756649

gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaa                                                                                               >  1:915187/1‑100 (MQ=255)

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

 aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaaa                                                                                              >  2:1096478/1‑100 (MQ=255)

 aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaaa                                                                                              >  2:558845/1‑100 (MQ=255)

  agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAc                                                                                             <  1:2980592/100‑1 (MQ=255)

   gaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAcc                                                                                            <  1:1071547/100‑1 (MQ=255)

   gaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAcc                                                                                            <  2:3336153/100‑1 (MQ=255)

    aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAccc                                                                                           <  2:3263799/100‑1 (MQ=255)

     agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCa                                                                                          >  2:2037942/1‑100 (MQ=255)

     agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCa                                                                                          >  1:1483939/1‑100 (MQ=255)

      gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCAc                                                                                         <  1:1185059/100‑1 (MQ=255)

       ttATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACg                                                                                        >  2:2533826/1‑100 (MQ=255)

        tATCTTCGCCCAGCAGATGATGATCGAAAAATGTATCTGTGCACGAAAAGTGGTTATTACTGCGAGCCAGATGCTGGATTCCATGATCTAacacac                                                                                           >  1:324091/1‑94 (MQ=255)

        tATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGc                                                                                       >  1:581754/1‑100 (MQ=255)

         aTCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGcc                                                                                      <  1:3406981/100‑1 (MQ=255)

          tCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGccc                                                                                     >  1:1086809/1‑100 (MQ=255)

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

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

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

                  cAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgcg                                                                             <  2:581756/100‑1 (MQ=255)

                  cAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgcg                                                                             <  1:2207865/100‑1 (MQ=255)

                     aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAg                                                                          <  1:1312589/100‑1 (MQ=255)

                        atgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAg                                                                       >  1:1620752/1‑100 (MQ=255)

                        atgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAg                                                                       <  1:1843791/100‑1 (MQ=255)

                           atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgcgagatcggaa                                                                    >  2:60067/1‑91 (MQ=255)

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

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

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

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

                             gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGATGGATTCCACGATCTAAAACCCACACCCGACTCGCGCAGAAGCCGGt                                                                  >  1:2376787/1‑100 (MQ=255)

                               tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                                 <  1:1087420/99‑1 (MQ=255)

                               tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                                                <  1:2015568/100‑1 (MQ=255)

                                cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                                               >  1:2985779/1‑100 (MQ=255)

                                cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                                               <  2:788889/100‑1 (MQ=255)

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

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

                                  aaaaaTGTATTCGTGCACGTAAAGTCGTTAGCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                                 <  1:3056761/96‑1 (MQ=255)

                                  aaaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt                                                             <  2:2494547/100‑1 (MQ=255)

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

                                     aaTGTATCCGTGAACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                                          <  1:2867596/100‑1 (MQ=255)

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

                                       tGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCaa                                                        >  2:1808684/1‑100 (MQ=255)

                                         tATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCagat                                                      >  2:1432500/1‑97 (MQ=255)

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

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

                                            ccGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGcc                                                   >  1:1167720/1‑100 (MQ=255)

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

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

                                                 cACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCt                                              <  1:514716/100‑1 (MQ=255)

                                                  aCGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTc                                             >  1:1209868/1‑100 (MQ=255)

                                                   cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                            <  1:433897/100‑1 (MQ=255)

                                                   cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                            >  2:659679/1‑100 (MQ=255)

                                                   cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                            <  1:1592964/100‑1 (MQ=255)

                                                   cGTAAAGTCATTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                            >  1:453963/1‑100 (MQ=255)

                                                     tAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                          >  1:3483452/1‑100 (MQ=255)

                                                     tAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                          >  2:2669698/1‑100 (MQ=255)

                                                      aaaGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACg                                         <  2:187415/100‑1 (MQ=255)

                                                          tCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTan                                     <  2:883066/100‑2 (MQ=255)

                                                           cGTTATCACTGCGACCCAGATGCTGGATTCCATGACCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt                                    >  1:592707/1‑100 (MQ=255)

                                                            gTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTg                                   >  2:271052/1‑100 (MQ=255)

                                                              tATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc                                 <  1:1826508/100‑1 (MQ=255)

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

                                                                  aCTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                             >  2:427591/1‑100 (MQ=255)

                                                                  aCTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                             <  1:2989856/100‑1 (MQ=255)

                                                                  aCTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                             <  1:2037927/100‑1 (MQ=255)

                                                                  aCTCCGCCCCAGATCATGGCTTCCATGCTCCACAACAAACGCCCACCTCGCGAAGACGCCGGTAACGATGCAAACGCCATCCTCGACGGTACTTACGCa                              <  1:288240/99‑1 (MQ=255)

                                                                   cTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCGCGACTCGCGCAGAAGCCGGGGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt                            >  1:1921790/1‑100 (MQ=255)

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

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

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

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

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

                                                                          ccAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTg                     <  2:2874240/100‑1 (MQ=255)

                                                                           cAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGt                    >  1:955356/1‑100 (MQ=255)

                                                                           cAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGt                    >  2:2216978/1‑100 (MQ=255)

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

                                                                             gATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCt                  <  2:3426832/100‑1 (MQ=255)

                                                                              aTGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg                 >  2:3395666/1‑100 (MQ=255)

                                                                              aTGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg                 >  2:215589/1‑100 (MQ=255)

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

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

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

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

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

                                                                                 cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg              >  2:545816/1‑100 (MQ=255)

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

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

                                                                                   ggATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa            >  2:1415219/1‑100 (MQ=255)

                                                                                   ggATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa            >  1:2435015/1‑100 (MQ=255)

                                                                                   ggATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa            >  1:2405289/1‑100 (MQ=255)

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

                                                                                        ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGc       >  1:588327/1‑100 (MQ=255)

                                                                                        ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGc       >  1:641434/1‑100 (MQ=255)

                                                                                        ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGc       >  2:318497/1‑100 (MQ=255)

                                                                                        ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGc       >  2:1721939/1‑100 (MQ=255)

                                                                                        ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGc       >  2:85867/1‑100 (MQ=255)

                                                                                        ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCCACGGTACTGACGCAGTGATGCTGTCTGGTGAAccc         >  1:145198/1‑98 (MQ=255)

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

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

                                                                                             aTCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAg  >  1:2314408/1‑100 (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: `‑`