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%	972.4 / ‑4.3	280	coding (497/1458 nt)	ascF	PTS system enzyme IIBC
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (141/138); total (142/138)
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.99e-01

Only 100 of 313 total aligned reads displayed.

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑G‑CATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTG  >  NC_000913/2837951‑2838131

 aTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑g                                                                                          <  1:8321974/93‑1 (MQ=255)

  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc‑gg                                                                                         >  2:5337043/1‑92 (MQ=255)

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

       ggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCt                                                                                    >  2:1472681/1‑93 (MQ=255)

         cgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGc                                                                                  >  1:1464173/1‑93 (MQ=255)

          gTGCTGACAAAAGGCTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTCCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGcc                                                                                 <  2:2570469/93‑1 (MQ=255)

          gTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGcc                                                                                 >  1:6842275/1‑93 (MQ=255)

          gTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGcc                                                                                 <  2:7185120/93‑1 (MQ=255)

           tGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCcg                                                                                <  1:7852891/93‑1 (MQ=255)

           tGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCcg                                                                                <  1:3522571/93‑1 (MQ=255)

            gCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCcgc                                                                               <  2:8483878/93‑1 (MQ=255)

            gCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCcgc                                                                               <  2:6860858/93‑1 (MQ=255)

              tGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCcgccc                                                                             >  2:6744835/1‑91 (MQ=255)

                aCAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTCCCTGCCGCTGATGGTCGCGG‑CATCTGCcgccatc                                                                           <  2:3465987/93‑5 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCa                                                                         >  2:611949/1‑93 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCa                                                                         <  1:3783010/93‑1 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCa                                                                         >  1:8697441/1‑93 (MQ=255)

                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCa                                                                         >  1:5859123/1‑93 (MQ=255)

                  aaaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCaa                                                                        >  1:8182756/1‑93 (MQ=255)

                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAAt                                                                      <  2:134985/93‑1 (MQ=255)

                      ggATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAAttt                                                                    >  2:5393643/1‑93 (MQ=255)

                       gATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTa                                                                   <  1:4486569/93‑1 (MQ=255)

                        aTCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTaa                                                                  >  1:1213849/1‑93 (MQ=255)

                         tCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTaaa                                                                 >  1:747394/1‑93 (MQ=255)

                         tCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTaaa                                                                 <  1:6945542/93‑1 (MQ=255)

                         tCGCCGACCTTAACCACTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTTCTGCCGCTGATGGTCGCGG‑CATCCGCCGC‑CATCAAATTTaaa                                                                 >  2:5794749/1‑93 (MQ=255)

                            ccGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAAcc                                                              <  1:8382468/93‑1 (MQ=255)

                             cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCa                                                             >  2:8733691/1‑93 (MQ=255)

                             cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCa                                                             >  2:778885/1‑93 (MQ=255)

                             cGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCa                                                             >  2:7428849/1‑93 (MQ=255)

                              gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCaa                                                            >  1:8042579/1‑93 (MQ=255)

                              gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCaa                                                            >  1:5648975/1‑93 (MQ=255)

                              gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CACCTGCCGC‑CATCAAATTTAAAACCaa                                                            >  1:1869876/1‑93 (MQ=255)

                               aCCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAat                                                           >  1:3027507/1‑93 (MQ=255)

                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAata                                                          >  2:5864182/1‑93 (MQ=255)

                                ccTTAACCATTCTGAATGTGATTGGCGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAata                                                          <  1:1066316/93‑1 (MQ=255)

                                ccTTAACCATTCTGAAAGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAata                                                          >  2:1264348/1‑93 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAatat                                                         >  2:6704110/1‑93 (MQ=255)

                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAatat                                                         <  2:901514/93‑1 (MQ=255)

                                  ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATg                                                        <  2:3756550/93‑1 (MQ=255)

                                    aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                      <  2:439317/93‑1 (MQ=255)

                                    aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                      <  2:4824836/93‑1 (MQ=255)

                                    aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTc                                                      >  2:8144250/1‑93 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                     >  1:408921/1‑93 (MQ=255)

                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCg                                                     >  1:4906138/1‑93 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                    >  1:4356507/1‑93 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                    >  2:4040382/1‑93 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                    >  1:2583489/1‑93 (MQ=255)

                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGc                                                    <  1:3116826/93‑1 (MQ=255)

                                       cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                   >  1:7759249/1‑93 (MQ=255)

                                       cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCt                                                   >  2:152056/1‑93 (MQ=255)

                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                  >  2:6253810/1‑93 (MQ=255)

                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                  >  2:5030632/1‑93 (MQ=255)

                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTg                                                  <  1:4244973/93‑1 (MQ=255)

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                 <  1:6545920/93‑1 (MQ=255)

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                 >  2:8383586/1‑93 (MQ=255)

                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTgg                                                 >  1:859163/1‑93 (MQ=255)

                                            tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGa                                              <  2:594475/93‑1 (MQ=255)

                                            tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGa                                              <  2:6356040/93‑1 (MQ=255)

                                                tgtgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAACCCAATATGTCGCTGGCGATTGc                                          >  1:7348129/1‑93 (MQ=255)

                                                  tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                        <  1:220619/93‑1 (MQ=255)

                                                  tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                        >  1:5304216/1‑93 (MQ=255)

                                                  tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                        <  2:1139053/93‑1 (MQ=255)

                                                  tgATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCAGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCa                                        <  2:2997418/93‑1 (MQ=255)

                                                     ttGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcg                                     <  2:1793013/93‑1 (MQ=255)

                                                       ggTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgcg                                   <  1:2159668/93‑1 (MQ=255)

                                                         tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                                 <  2:3294497/93‑1 (MQ=255)

                                                          gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCgggg                                >  2:5536977/1‑92 (MQ=255)

                                                           aCGTTGCTGTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtg                               <  2:6462532/93‑1 (MQ=255)

                                                            cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                              <  1:5462223/93‑1 (MQ=255)

                                                            cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                              <  1:6096464/93‑1 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                            >  1:1657192/1‑93 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                            >  2:8345115/1‑93 (MQ=255)

                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                            >  2:2638659/1‑93 (MQ=255)

                                                               tGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt                           >  1:7270520/1‑93 (MQ=255)

                                                               tGCTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt                           <  1:944154/93‑1 (MQ=255)

                                                                 cTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTgg                         >  1:3741252/1‑93 (MQ=255)

                                                                 cTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTgg                         <  2:1618604/93‑1 (MQ=255)

                                                                 cTTTCTTCTTCCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGCGTGCTgg                         >  2:1923035/1‑93 (MQ=255)

                                                                   ttcttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                       >  1:2043198/1‑93 (MQ=255)

                                                                    tcttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                      <  1:6365881/93‑1 (MQ=255)

                                                                    tcttcttcCTGCCGCTGATGGTCGCGG‑CAGCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                      <  1:2932618/93‑1 (MQ=255)

                                                                     cttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa                     >  2:7623984/1‑93 (MQ=255)

                                                                     cttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa                     >  2:4577583/1‑93 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                    >  2:3570107/1‑93 (MQ=255)

                                                                      ttcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACAt                    >  2:7991515/1‑93 (MQ=255)

                                                                       tcttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                   >  1:7048621/1‑93 (MQ=255)

                                                                        cttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATcc                  >  2:2240884/1‑93 (MQ=255)

                                                                         ttcCTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCg                 <  1:3862016/93‑1 (MQ=255)

                                                                          tcCTGCCGATGATGGCCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                <  2:4310244/93‑1 (MQ=255)

                                                                            cTGCCGCTGATGTTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGt              <  2:4986881/93‑1 (MQ=255)

                                                                            cTGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCACGGGTGTGCTGGTACATCCCAGt              >  1:4738968/1‑93 (MQ=255)

                                                                             tGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt             >  2:5735293/1‑93 (MQ=255)

                                                                             tGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt             >  2:2937197/1‑93 (MQ=255)

                                                                             tGCCGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt             <  1:2826440/93‑1 (MQ=255)

                                                                                cGCTGATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa          <  1:5357542/93‑1 (MQ=255)

                                                                                  cTGATGGTCGC‑GTCATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt        <  2:3983076/93‑1 (MQ=255)

                                                                                    gATGGTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGa      >  1:6726702/1‑93 (MQ=255)

                                                                                       ggTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGAGTGCTGGTACATCCGAGGTTTATTGAACt   >  2:6805324/1‑93 (MQ=255)

                                                                                        gTCGCGG‑CATCTGCCGC‑CATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTg  >  1:5156667/1‑93 (MQ=255)

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGC‑G‑CATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTG  >  NC_000913/2837951‑2838131

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: `‑`