BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,838,044	+G	coding (498/1458 nt)	ascF →	PTS system enzyme IIBC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,838,043	1	.	G	99.2%	411.9 / ‑4.9	128	A166G (GCG→GGG)	ascF	PTS system enzyme IIBC
Reads supporting (aligned to +/- strand): ref base . (0/0); major base G (65/62); minor base A (1/0); total (66/62)
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.52e-01

Only 100 of 138 total aligned reads displayed.

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACT  >  NC_000913/2837951‑2838130

gatgatCGGTGTGCTGACAAAAGGATCGCCGACCTTACCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                                          <  2:1076088/93‑1 (MQ=255)

gaTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                                          >  2:2678490/1‑93 (MQ=255)

gaTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                                          >  1:3916775/1‑93 (MQ=255)

  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑gg                                                                                       <  1:2906708/93‑2 (MQ=255)

  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑gg                                                                                       <  1:2251608/93‑2 (MQ=255)

    aGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggca                                                                                     <  1:260413/93‑4 (MQ=255)

     gcggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggcat                                                                                    >  2:1431812/1‑89 (MQ=255)

     gcggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggcat                                                                                    >  2:2104926/1‑89 (MQ=255)

       ggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCt                                                                                   >  2:4042345/1‑93 (MQ=255)

        gcgTGGTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTg                                                                                  <  1:1894195/93‑1 (MQ=255)

        gcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTg                                                                                  <  1:1894196/93‑1 (MQ=255)

        gcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCAGCCGCTGATGGTCGCGGCATCTg                                                                                  >  2:3390693/1‑93 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGc                                                                                 <  1:274472/93‑1 (MQ=255)

           tGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcg                                                                               >  2:809413/1‑93 (MQ=255)

           tGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcg                                                                               >  2:3439156/1‑93 (MQ=255)

             cTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgcc                                                                             <  1:2158757/93‑2 (MQ=255)

               gACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgccat                                                                           >  1:4499081/1‑90 (MQ=255)

               gACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgccat                                                                           <  1:688041/93‑4 (MQ=255)

                aCAAAAGGATCGCCTTCCTTAACCATTCTGAATGTGATTGGTAACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgccatc                                                                          <  2:1622260/93‑5 (MQ=255)

                aCAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgccatc                                                                          <  2:3330769/93‑5 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCa                                                                        >  2:4339260/1‑93 (MQ=255)

                   aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaaa                                                                      <  2:1764994/93‑1 (MQ=255)

                   aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaaa                                                                      >  1:2998620/1‑93 (MQ=255)

                     aGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAtt                                                                    <  1:4342880/93‑1 (MQ=255)

                      ggATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAttt                                                                   >  1:1057500/1‑93 (MQ=255)

                         tCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaa                                                                >  2:2890595/1‑93 (MQ=255)

                         tCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaa                                                                >  2:2890596/1‑93 (MQ=255)

                         tCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaa                                                                <  2:2822705/93‑1 (MQ=255)

                           gCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAAc                                                              <  1:525738/93‑1 (MQ=255)

                              gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCaa                                                           <  2:3326029/93‑1 (MQ=255)

                              gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCaa                                                           >  2:2224932/1‑93 (MQ=255)

                              gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCaa                                                           >  2:4137724/1‑93 (MQ=255)

                              gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCaa                                                           >  2:2508389/1‑93 (MQ=255)

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAata                                                         >  1:1655408/1‑93 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGCCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAatat                                                        <  1:4357623/93‑1 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAatat                                                        <  2:3647135/93‑1 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAatat                                                        >  2:646978/1‑93 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAatat                                                        <  2:797924/93‑1 (MQ=255)

                                  ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                       >  1:2518689/1‑93 (MQ=255)

                                   tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGt                                                      >  1:3643110/1‑93 (MQ=255)

                                   tAACCATTCGGAATGTGATTGGTGACGGTGCTTTCTTCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGt                                                      <  1:856050/93‑1 (MQ=255)

                                    aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                     <  2:505634/93‑1 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                    >  2:2080888/1‑93 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                    <  2:620011/93‑1 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                    >  1:1076158/1‑93 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGAATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                    >  2:3503096/1‑93 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                   <  1:2406693/93‑1 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                   >  1:229542/1‑93 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCC‑CCCATCAAATTTAAAACCCATATGTccc                                                   >  1:229543/1‑91 (MQ=255)

                                       cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                  <  1:2605719/93‑1 (MQ=255)

                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                 <  1:681381/93‑1 (MQ=255)

                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                 <  1:1578679/93‑1 (MQ=255)

                                            tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGa                                             >  1:2866219/1‑93 (MQ=255)

                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                           <  2:4517299/93‑1 (MQ=255)

                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                           <  2:1210978/93‑1 (MQ=255)

                                               aTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTg                                          <  1:3651635/93‑1 (MQ=255)

                                                 gtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGcc                                        <  1:2018617/93‑1 (MQ=255)

                                                  tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                       <  2:1455234/93‑1 (MQ=255)

                                                   gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                      >  1:1526979/1‑93 (MQ=255)

                                                   gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                      >  1:825209/1‑93 (MQ=255)

                                                   gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                      >  2:2556621/1‑93 (MQ=255)

                                                    aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                     >  1:4039464/1‑93 (MQ=255)

                                                    aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                     <  2:1884240/93‑1 (MQ=255)

                                                    aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                     <  2:2242506/93‑1 (MQ=255)

                                                       ggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgcg                                  >  1:580898/1‑93 (MQ=255)

                                                         tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                                <  1:2352257/93‑1 (MQ=255)

                                                         tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                                >  2:2615327/1‑93 (MQ=255)

                                                         tGACGGTGCTTTCTTCGTCCTGCCGCTGATGGGCGCGGCATCTGCCG‑ACATCAGATTTAAAACCAATAGGTCGGTGGCGATTGGCATCGCggg                                >  2:2579225/1‑93 (MQ=255)

                                                          gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgt                               >  2:2846195/1‑93 (MQ=255)

                                                           aCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                              <  2:3011103/93‑1 (MQ=255)

                                                           aCGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCG‑ACATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                              <  1:1112078/93‑1 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                           >  2:3192737/1‑93 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                           <  1:3068977/93‑1 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                           >  2:1909451/1‑93 (MQ=255)

                                                                gCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                         <  2:633524/93‑1 (MQ=255)

                                                                gCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGGGTGCTg                         >  2:1066143/1‑93 (MQ=255)

                                                                 cTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTgg                        >  1:4021603/1‑93 (MQ=255)

                                                                 cTTTCTTCTTCCTGCCGCTGATGGTCGCAGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTgg                        >  2:1894099/1‑93 (MQ=255)

                                                                   ttcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                      <  2:333927/93‑1 (MQ=255)

                                                                    tcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                     >  1:3846950/1‑93 (MQ=255)

                                                                     cttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa                    <  2:1710831/93‑1 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                   >  1:2950079/1‑93 (MQ=255)

                                                                       tcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                  <  1:4736435/93‑1 (MQ=255)

                                                                       tcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                  <  1:4075890/93‑1 (MQ=255)

                                                                       tcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                  <  2:4579037/93‑1 (MQ=255)

                                                                       tcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                  >  1:2758552/1‑93 (MQ=255)

                                                                         ttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCg                >  1:1338801/1‑93 (MQ=255)

                                                                          tcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa               <  2:2998620/93‑1 (MQ=255)

                                                                           ccTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAg              >  2:481300/1‑93 (MQ=255)

                                                                                cGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa         >  1:2953477/1‑93 (MQ=255)

                                                                                cGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa         >  1:2880786/1‑93 (MQ=255)

                                                                                 gCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTATATCCGAGTTTTAt        >  1:3907259/1‑93 (MQ=255)

                                                                                 gCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt        <  1:128028/93‑1 (MQ=255)

                                                                                   tGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTg      >  1:3068918/1‑93 (MQ=255)

                                                                                   tGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTg      >  2:4310763/1‑93 (MQ=255)

                                                                                    gATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGa     >  1:2590076/1‑93 (MQ=255)

                                                                                    gATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGa     >  2:1274065/1‑93 (MQ=255)

                                                                                    gATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGa     <  2:4535479/93‑1 (MQ=255)

                                                                                     aTGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGaa    <  2:2668541/93‑1 (MQ=255)

                                                                                      tGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAAc   <  1:1432658/93‑1 (MQ=255)

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACT  >  NC_000913/2837951‑2838130

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