BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,261,155	A→T	D186V (GAT→GTT)	yjbM →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,155	0	A	T	93.2%	387.4 / 15.9	120	D186V (GAT→GTT)	yjbM	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (3/5); new base T (53/56); total (58/62)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 120 total aligned reads displayed.

ACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATAC  >  NC_000913/4261059‑4261250

aCCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTcc                                                                                              >  2:2038504/1‑100 (MQ=255)

aCCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATcc                                                                                              >  1:1037580/1‑100 (MQ=255)

 ccTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCt                                                                                             <  1:815991/100‑1 (MQ=255)

 ccTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCt                                                                                             <  2:1165894/100‑1 (MQ=255)

  cTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTg                                                                                            <  1:59860/100‑1 (MQ=255)

    gggAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGat                                                                                          >  1:380647/1‑100 (MQ=255)

     ggAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGata                                                                                         >  2:1216504/1‑100 (MQ=255)

      gagagaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGatat                                                                                        <  1:239093/100‑1 (MQ=255)

       agagaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATAtt                                                                                       <  1:1368558/100‑1 (MQ=255)

          gaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTa                                                                                    >  1:1585226/1‑100 (MQ=255)

          gaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTa                                                                                    >  1:1630874/1‑100 (MQ=255)

          gaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTa                                                                                    >  1:1700801/1‑100 (MQ=255)

           aTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTaa                                                                                   >  2:76380/1‑100 (MQ=255)

           aTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTaa                                                                                   <  2:1321806/100‑1 (MQ=255)

           aTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTaa                                                                                   <  1:1547779/100‑1 (MQ=255)

            tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAg                                                                                  >  1:839321/1‑100 (MQ=255)

            tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAg                                                                                  >  2:1216586/1‑100 (MQ=255)

            tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAg                                                                                  <  2:417760/100‑1 (MQ=255)

            tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACCCTCGTTCCTGATATTTTAAg                                                                                  <  2:1429042/100‑1 (MQ=255)

             gAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGc                                                                                 >  1:1720140/1‑100 (MQ=255)

             gAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGc                                                                                 >  2:1212851/1‑100 (MQ=255)

             gAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGc                                                                                 >  2:17867/1‑100 (MQ=255)

             gAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGc                                                                                 >  2:156466/1‑100 (MQ=255)

               gcggccTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCaa                                                                               <  1:945174/95‑1 (MQ=255)

                gCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCaaa                                                                              <  2:1037578/100‑1 (MQ=255)

                 cGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCaaaa                                                                             >  1:2024436/1‑100 (MQ=255)

                 cGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCaaaa                                                                             >  1:1868344/1‑100 (MQ=255)

                  gTCTGGCTGATGGGGGAGAAGTGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCaaaa                                                                             >  2:1425737/1‑99 (MQ=255)

                  gTCTGGCTGATGGGGGAGAAATGGATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAAt                                                                            <  1:1421361/100‑1 (MQ=255)

                  gTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAAt                                                                            <  1:1371228/100‑1 (MQ=255)

                     tGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCg                                                                         <  1:238614/100‑1 (MQ=255)

                     tGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCg                                                                         <  1:1402581/100‑1 (MQ=255)

                     tGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCg                                                                         <  1:1402873/100‑1 (MQ=255)

                      ggCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCga                                                                        <  1:2115787/100‑1 (MQ=255)

                      ggCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCga                                                                        >  2:204201/1‑100 (MQ=255)

                       gCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCgag                                                                       <  2:160628/100‑1 (MQ=255)

                         tGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGAt                                                                     >  1:1015724/1‑100 (MQ=255)

                         ccgatctggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGAt                                                                     <  1:37721/93‑1 (MQ=255)

                          gATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGAtg                                                                    >  1:670455/1‑100 (MQ=255)

                              ggggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGACGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGAtgctgc                                                                >  1:1550025/1‑100 (MQ=255)

                               gggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGca                                                               <  1:648585/100‑1 (MQ=255)

                               gggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGca                                                               >  1:2010135/1‑100 (MQ=255)

                                ggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGAtagatcgg                                                              >  2:37719/1‑93 (MQ=255)

                                ggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                                                              >  1:743110/1‑100 (MQ=255)

                                 gagaAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGcacc                                                             >  2:873543/1‑100 (MQ=255)

                                   gaAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGcaccac                                                           <  2:1720135/100‑1 (MQ=255)

                                    aaaTGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACcaca                                                          <  1:265829/100‑1 (MQ=255)

                                    aaaTGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACcaca                                                          <  1:603582/100‑1 (MQ=255)

                                     aaTGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaa                                                         >  1:1286677/1‑100 (MQ=255)

                                      aTGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaat                                                        <  2:1753205/100‑1 (MQ=255)

                                        gAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaatta                                                      >  2:601917/1‑100 (MQ=255)

                                         aaTTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaa                                                     <  2:1015722/100‑1 (MQ=255)

                                          aTTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                    >  2:684374/1‑100 (MQ=255)

                                          aTTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                    >  2:54778/1‑100 (MQ=255)

                                          aTTTTGATGGCTATCTGGATAAGATGACACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                    >  1:1425858/1‑100 (MQ=255)

                                           ttttGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaatt                                                   <  2:1585221/100‑1 (MQ=255)

                                             ttGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAg                                                 <  2:335861/100‑1 (MQ=255)

                                               gATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCt                                               >  1:1859880/1‑100 (MQ=255)

                                                   gCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcgc                                           <  2:1559726/100‑1 (MQ=255)

                                                    cTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcgcg                                          >  2:1714109/1‑100 (MQ=255)

                                                    cTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcgcg                                          <  2:220963/100‑1 (MQ=255)

                                                      aacTGGATAAGATATCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGaa                                        >  1:1574790/3‑100 (MQ=255)

                                                      aTCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGaa                                        >  2:2118288/1‑100 (MQ=255)

                                                      aTCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGaa                                        >  1:1880568/1‑100 (MQ=255)

                                                         tGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTc                                     >  1:734268/1‑100 (MQ=255)

                                                           gATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGt                                   <  1:952802/100‑1 (MQ=255)

                                                            aTAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGtt                                  <  1:1936409/100‑1 (MQ=255)

                                                            aTAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTCAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGtt                                  <  1:1749453/100‑1 (MQ=255)

                                                             tAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGttt                                 >  1:275075/1‑100 (MQ=255)

                                                              aaGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTg                                <  2:925073/100‑1 (MQ=255)

                                                              aaGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTg                                <  2:807347/100‑1 (MQ=255)

                                                              aaGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTg                                <  2:2010128/100‑1 (MQ=255)

                                                               aGATGTCACATCTGGTGAATGAAGGAACGCTCGGTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCTCAATTAATTAGCTTCGCGAAGTCGTTTGa                               >  2:1411051/1‑100 (MQ=255)

                                                                gATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAc                              >  2:320741/1‑100 (MQ=255)

                                                                gATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAc                              <  1:1212854/100‑1 (MQ=255)

                                                                    tCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCaa                          >  1:277518/1‑100 (MQ=255)

                                                                     cacaTCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                         >  2:345184/1‑100 (MQ=255)

                                                                           tGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCaaa                   >  2:2199456/1‑100 (MQ=255)

                                                                           tGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCaaa                   >  2:1466602/1‑100 (MQ=255)

                                                                            ggTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAg                  <  1:320743/100‑1 (MQ=255)

                                                                            ggTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAg                  >  2:1976670/1‑100 (MQ=255)

                                                                             gTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAgg                 <  2:275073/100‑1 (MQ=255)

                                                                              tgaatgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGa                <  2:839319/100‑1 (MQ=255)

                                                                               gaatgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGaa               >  2:1790780/1‑100 (MQ=255)

                                                                                 atgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAggagga             >  2:633260/1‑100 (MQ=255)

                                                                                  tgacggAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAagag            <  2:1550021/96‑1 (MQ=255)

                                                                                   gaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAgaga           <  2:2068258/100‑1 (MQ=255)

                                                                                   gaaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAgaga           >  2:1270100/1‑100 (MQ=255)

                                                                                    aaGGAACGCTCGTTCCTGATATTTTACGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAt          >  1:1425786/1‑100 (MQ=255)

                                                                                    aaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAt          <  1:1976677/100‑1 (MQ=255)

                                                                                    aaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAt          <  1:967085/100‑1 (MQ=255)

                                                                                     aGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCGTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAtt         <  2:380898/100‑1 (MQ=255)

                                                                                      ggAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTa        >  1:473830/1‑100 (MQ=255)

                                                                                          cGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAAt    <  1:2042482/100‑1 (MQ=255)

                                                                                          cGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAAt    <  1:1094828/100‑1 (MQ=255)

                                                                                          cGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAAt    <  1:76382/100‑1 (MQ=255)

                                                                                           gCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATa   <  1:345186/100‑1 (MQ=255)

                                                                                            cTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATAc  <  1:17868/100‑1 (MQ=255)

                                                                                            cTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATAc  >  2:1092472/1‑100 (MQ=255)

                                                                                            cTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATAc  >  2:1883173/1‑100 (MQ=255)

ACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATAC  >  NC_000913/4261059‑4261250

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