BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,229,642	A→G	L391L (TTG→CTG)	lysC ←	aspartokinase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,229,642	0	A	G	99.6%	808.7 / ‑4.5	226	L391L (TTG→CTG)	lysC	aspartokinase III
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (130/95); total (131/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

Only 100 of 239 total aligned reads displayed.

CATAACAAATCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATCGCCAG  >  NC_000913/4229552‑4229731

cATAACAAATCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGcg                                                                                         >  2:29490/1‑93 (MQ=255)

cATAACAAATCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGcg                                                                                         >  1:3100333/1‑93 (MQ=255)

cATAACAAATCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGcg                                                                                         >  1:1021263/1‑93 (MQ=255)

cATAACAAATCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGcg                                                                                         >  2:684962/1‑93 (MQ=255)

 aTAACAAATCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAgcgc                                                                                        <  2:4684507/93‑1 (MQ=255)

 aTAACAAATCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAgcgc                                                                                        <  1:2750286/93‑1 (MQ=255)

  tAACAAATCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGTCAATCAGcgcg                                                                                       >  2:6076722/1‑93 (MQ=255)

      aaaTCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCa                                                                                   >  2:6009382/1‑93 (MQ=255)

      aaaTCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCa                                                                                   <  2:5510745/93‑1 (MQ=255)

      aaaTCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCa                                                                                   <  1:1361082/93‑1 (MQ=255)

       aaTCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAg                                                                                  >  2:509296/1‑93 (MQ=255)

       aaTCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAg                                                                                  >  2:4502913/1‑93 (MQ=255)

       aaTCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAg                                                                                  >  1:2693045/1‑93 (MQ=255)

       aaTCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAg                                                                                  >  1:5676147/1‑93 (MQ=255)

         tCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAgcg                                                                                >  2:2620393/1‑93 (MQ=255)

            tGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCa                                                                             >  1:1462994/1‑93 (MQ=255)

            tGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCa                                                                             >  1:5020338/1‑93 (MQ=255)

            tGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCa                                                                             <  1:4118438/93‑1 (MQ=255)

            tGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCa                                                                             <  2:1240461/93‑1 (MQ=255)

               gAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGAc                                                                          >  2:4997451/1‑93 (MQ=255)

                  tGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACctt                                                                       <  2:38552/93‑1 (MQ=255)

                   gTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACcttc                                                                      <  1:2053399/93‑1 (MQ=255)

                      gAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCttcttc                                                                   <  1:4328146/93‑1 (MQ=255)

                      gAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCttcttc                                                                   <  2:4651592/93‑1 (MQ=255)

                       aaCGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTtcc                                                                  <  1:3146458/93‑1 (MQ=255)

                        aCGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTtcca                                                                 >  1:4396169/1‑93 (MQ=255)

                        aCGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTtcca                                                                 <  2:203118/93‑1 (MQ=255)

                        aCGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTtcca                                                                 >  2:3700269/1‑93 (MQ=255)

                          ggTTCCAGTACGCCGAGTACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTtccacc                                                               >  1:5856616/1‑93 (MQ=255)

                           gTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTtccacct                                                              >  2:3142223/1‑93 (MQ=255)

                           gTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTtccacct                                                              >  1:6076841/1‑93 (MQ=255)

                             tCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTtccacctcc                                                            >  1:4943657/1‑93 (MQ=255)

                              ccAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTtccacctcca                                                           >  1:1165652/1‑93 (MQ=255)

                                  tACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCg                                                       >  2:1093759/1‑93 (MQ=255)

                                      ccGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGacac                                                   <  2:5958799/93‑1 (MQ=255)

                                      ccGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGacac                                                   <  1:5101258/93‑1 (MQ=255)

                                       cGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCCCCTCCACCCGAcaca                                                  <  2:845186/93‑1 (MQ=255)

                                        gagtACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAg                                                 <  1:4608722/90‑1 (MQ=255)

                                        gAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAg                                                 >  2:1722315/1‑93 (MQ=255)

                                        gAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAg                                                 <  2:5654057/93‑1 (MQ=255)

                                         aaTACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCGCCTCCACCCGACACAGt                                                <  1:3375747/93‑1 (MQ=255)

                                         aaTACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGt                                                <  1:231986/93‑1 (MQ=255)

                                         aaTACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGt                                                <  1:3068679/93‑1 (MQ=255)

                                          aTACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTg                                               >  1:1109471/1‑93 (MQ=255)

                                           tACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGc                                              >  2:2308377/1‑93 (MQ=255)

                                           tACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGc                                              >  1:1859352/1‑93 (MQ=255)

                                            tccTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCg                                             <  1:4912941/92‑1 (MQ=255)

                                            aCCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCg                                             >  1:618085/1‑93 (MQ=255)

                                            aCCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCg                                             >  1:3899927/1‑93 (MQ=255)

                                             ccTCTTTTCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCgg                                            <  2:2256138/93‑1 (MQ=255)

                                               tctTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGaa                                          >  2:1778891/1‑93 (MQ=255)

                                               tctTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGaa                                          >  2:2555749/1‑93 (MQ=255)

                                                ctTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGaaa                                         <  2:4056373/93‑1 (MQ=255)

                                                 tttGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAg                                        <  2:3228726/93‑1 (MQ=255)

                                                  ttGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGc                                       <  1:4613034/93‑1 (MQ=255)

                                                    gCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTc                                     <  2:260631/93‑1 (MQ=255)

                                                    gCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTc                                     <  2:1497970/93‑1 (MQ=255)

                                                     ccAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                    >  1:5160973/1‑93 (MQ=255)

                                                     ccAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                    >  1:2754620/1‑93 (MQ=255)

                                                     ccAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTcc                                    <  1:3014144/93‑1 (MQ=255)

                                                      cAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCa                                   >  1:600400/1‑93 (MQ=255)

                                                        aCGCCGCAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATc                                 <  2:1276327/93‑1 (MQ=255)

                                                         cgccgcAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATca                                <  2:3404360/93‑1 (MQ=255)

                                                           ccgcAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcagc                              <  2:1359737/93‑1 (MQ=255)

                                                           ccgcAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcagc                              >  2:1497755/1‑93 (MQ=255)

                                                            cgcAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcagca                             >  2:4420643/1‑93 (MQ=255)

                                                            cgcAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcagca                             >  2:5291364/1‑93 (MQ=255)

                                                            cgcAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcagca                             <  1:932248/93‑1 (MQ=255)

                                                             gcAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcagcag                            >  1:829053/1‑93 (MQ=255)

                                                             gcAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcagcag                            <  1:3950829/93‑1 (MQ=255)

                                                             gcAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATcagcag                            <  1:4533385/93‑1 (MQ=255)

                                                              cAGGCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCaga                           >  1:758648/1‑93 (MQ=255)

                                                                ggCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAgaga                         >  1:3115460/1‑93 (MQ=255)

                                                                ggCTTTTGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCAACAGCAgaga                         <  1:3600261/93‑1 (MQ=255)

                                                                ggCTTTTGACAGGGCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAgaga                         >  1:5551957/1‑93 (MQ=255)

                                                                   ttttGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTg                      >  2:1654411/1‑93 (MQ=255)

                                                                   ttttGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTg                      >  1:5591919/1‑93 (MQ=255)

                                                                    tttGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGc                     >  2:3909397/1‑93 (MQ=255)

                                                                    tgtGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGc                     <  1:172398/91‑1 (MQ=255)

                                                                      tGACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGt                   <  2:285968/93‑1 (MQ=255)

                                                                       gACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTc                  <  2:4850875/93‑1 (MQ=255)

                                                                       gACAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTc                  <  2:3117284/93‑1 (MQ=255)

                                                                        aCAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCa                 <  2:1901087/93‑1 (MQ=255)

                                                                        aCAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCa                 >  2:3430692/1‑93 (MQ=255)

                                                                         cAGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAg                >  1:2151645/1‑93 (MQ=255)

                                                                          agctCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGCCACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGc               <  1:4717248/90‑1 (MQ=255)

                                                                          aGGTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGc               >  2:5168974/1‑93 (MQ=255)

                                                                            gTCATTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCaa             <  1:1950467/93‑1 (MQ=255)

                                                                               aTTGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGt          >  1:5378352/1‑93 (MQ=255)

                                                                                ttGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTa         <  1:3720738/93‑1 (MQ=255)

                                                                                ttGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTa         >  1:5808631/1‑93 (MQ=255)

                                                                                ttGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTa         <  1:2462088/93‑1 (MQ=255)

                                                                                ttGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTa         >  1:1069022/1‑93 (MQ=255)

                                                                                 tGCCAATCAGCGCGACCAGCGCCAGCCCTTCTTCCACCTCCCCCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTAt        <  2:1696210/93‑1 (MQ=255)

                                                                                 tGCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTAt        <  2:3782956/93‑1 (MQ=255)

                                                                                  gCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATc       <  2:4083877/93‑1 (MQ=255)

                                                                                  gCCAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATc       <  2:5722574/93‑1 (MQ=255)

                                                                                    cAATCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATCGc     >  1:4046453/1‑93 (MQ=255)

                                                                                      aTCAGCGCGACCAGCGCCAGACCGTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATCGCCa   <  1:5036160/93‑1 (MQ=255)

                                                                                       tCAGCGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATCGCCAg  >  1:3576881/1‑93 (MQ=255)

CATAACAAATCATGCGAATGTTGAACGGTTCCAGTACGCCGAATACCTCTTTGCCAACGCCGCAGGCTTTTGACAGGTCATTGCCAATCAACGCGACCAGCGCCAGACCTTCTTCCACCTCCACCCGACACAGTGCGGAAAGCTCCATCAGCAGAGATTGCGTCAGCAACGTATCGCCAG  >  NC_000913/4229552‑4229731

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