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.6%	741.9 / ‑4.0	227	coding (497/1458 nt)	ascF	PTS system enzyme IIBC
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (109/117); total (109/118)
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.93e-01

Only 100 of 252 total aligned reads displayed.

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑GCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACT  >  NC_000913/2837951‑2838130

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

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

gaTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                                          <  2:653644/93‑1 (MQ=255)

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

   gAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggc                                                                                      >  2:2501836/1‑91 (MQ=255)

    aGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑ggca                                                                                     >  1:2081934/1‑90 (MQ=255)

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

       ggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGATGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCt                                                                                   <  1:5117083/93‑1 (MQ=255)

       ggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCt                                                                                   <  2:2226433/93‑1 (MQ=255)

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

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

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

          gTGCTGACAAAAGGCTCGCCGCCCTTAACCATTCTGAATATGATTGGTGACGGTGCTTTCTTCTTACTGCCGCTGATGGTCGCGGCATCtctc                                                                                >  1:5779781/1‑90 (MQ=255)

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

                aCAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgccatc                                                                          >  2:4959228/1‑89 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCa                                                                        <  2:720479/93‑1 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCa                                                                        <  1:4327826/93‑1 (MQ=255)

                   aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGGGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCaaa                                                                      >  1:3488323/1‑93 (MQ=255)

                    aaGGATCGCGGCCCTTAACCATTCTGAATGTGATTGGTGACGGTCCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAt                                                                     <  1:2317057/93‑1 (MQ=255)

                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAt                                                                     >  1:400918/1‑93 (MQ=255)

                     aGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAtt                                                                    >  1:5330371/1‑93 (MQ=255)

                     aGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAAtt                                                                    >  1:4755790/1‑93 (MQ=255)

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

                       gATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTa                                                                  >  2:3346147/1‑93 (MQ=255)

                       gATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTa                                                                  <  1:1909654/93‑1 (MQ=255)

                        aTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaa                                                                 >  2:4569118/1‑93 (MQ=255)

                          cGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaaa                                                               >  1:5910245/1‑93 (MQ=255)

                          cGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTaaaa                                                               <  2:3717501/93‑1 (MQ=255)

                           gcggACTTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTCCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAAc                                                              <  2:433352/90‑1 (MQ=255)

                             cGACCTTAACCATTCTGAATGTGATTGTTGACGTTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCa                                                            <  1:937938/93‑1 (MQ=255)

                             cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCa                                                            >  1:3146043/1‑93 (MQ=255)

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

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

                               aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAat                                                          <  1:1244473/93‑1 (MQ=255)

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

                                  ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                       <  2:2809991/93‑1 (MQ=255)

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

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

                                    aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                     >  2:5227840/1‑93 (MQ=255)

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

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

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

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

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

                                      ccATTCTGAATGGGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                   <  1:5367450/93‑1 (MQ=255)

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

                                       cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                  >  1:5994260/1‑93 (MQ=255)

                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                 >  1:1790544/1‑93 (MQ=255)

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                <  2:5483223/93‑1 (MQ=255)

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                <  2:295026/93‑1 (MQ=255)

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                >  1:6247381/1‑93 (MQ=255)

                                          tCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTAGc                                               >  1:5754113/1‑93 (MQ=255)

                                           cTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCg                                              >  1:5496435/1‑93 (MQ=255)

                                            tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTCCCGC‑CATCAAATTTAAAACCAATATGTCGCTggggt                                             >  1:6444607/1‑90 (MQ=255)

                                             gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAt                                            <  1:3311171/93‑1 (MQ=255)

                                             gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAt                                            <  1:3205784/93‑1 (MQ=255)

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

                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                           >  2:5125059/1‑93 (MQ=255)

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

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

                                                   gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                      <  2:5832343/93‑1 (MQ=255)

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

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

                                                    aTTGGTGACGGTGCTGTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                     <  2:1869310/93‑1 (MQ=255)

                                                     ttGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                                    <  2:2706205/93‑1 (MQ=255)

                                                       ggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgcg                                  <  1:2466249/93‑1 (MQ=255)

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

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

                                                           tcGGTCCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                              <  1:2870814/92‑1 (MQ=255)

                                                            cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                             >  1:4423248/1‑93 (MQ=255)

                                                            cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGAGATTGCCATCGCGGgtgg                             >  1:1100862/1‑92 (MQ=255)

                                                             ggTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGtgtg                            >  2:5399694/1‑93 (MQ=255)

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

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

                                                               tccTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt                          <  2:648829/91‑1 (MQ=255)

                                                               ggCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt                          <  1:5439912/92‑1 (MQ=255)

                                                                gCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                         >  2:5668198/1‑93 (MQ=255)

                                                                gCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTg                         >  1:770819/1‑93 (MQ=255)

                                                                  tttCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                       <  2:1595616/93‑1 (MQ=255)

                                                                  tttCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGt                       <  1:3405394/93‑1 (MQ=255)

                                                                   ttcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                      >  1:1402430/1‑93 (MQ=255)

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

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

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

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

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

                                                                        cttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATcc                 <  1:2730340/93‑1 (MQ=255)

                                                                         ttcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCg                <  2:1389802/93‑1 (MQ=255)

                                                                          tcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa               >  1:4661792/1‑93 (MQ=255)

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

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

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

                                                                          tcCTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa               >  2:4506264/1‑93 (MQ=255)

                                                                            cTGCCGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt             >  1:2631737/1‑93 (MQ=255)

                                                                               ccGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtttt          >  2:76932/1‑93 (MQ=255)

                                                                                cGCTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa         <  1:5257839/93‑1 (MQ=255)

                                                                                  cTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt       >  2:5713574/1‑93 (MQ=255)

                                                                                  cTGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt       <  2:1468162/93‑1 (MQ=255)

                                                                                   tGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTg      <  2:1133322/93‑1 (MQ=255)

                                                                                   tGATGGTCGCGGCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTg      <  1:3691564/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: `‑`