BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,684	+A	coding (400/501 nt)	ytfI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,681	1	.	A	98.8%	929.0 / 1.2	248	coding (397/501 nt)	ytfI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/1); new base A (117/128); total (119/129)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 267 total aligned reads displayed.

GAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AAC‑AAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAT  >  NC_000913/4435589‑4435768

gAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑aac                                                                                          <  1:2385044/93‑1 (MQ=255)

gAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑aac                                                                                          <  2:2204839/93‑1 (MQ=255)

  ggTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑aac‑aa                                                                                       >  1:2058613/1‑93 (MQ=255)

    tAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑aac‑aaaa                                                                                     <  2:6014908/93‑2 (MQ=255)

     aGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑aac‑aaaac                                                                                    <  1:1104294/93‑3 (MQ=255)

             tGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGAtt                                                                             <  2:5787925/93‑1 (MQ=255)

              gAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTg                                                                            >  2:2575357/1‑93 (MQ=255)

              gAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTg                                                                            >  2:3093375/1‑93 (MQ=255)

              gAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTg                                                                            <  2:6726705/93‑1 (MQ=255)

               aaTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGt                                                                           <  2:2093906/93‑1 (MQ=255)

               aaTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGt                                                                           >  1:885375/1‑93 (MQ=255)

               aaTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGt                                                                           >  1:4430439/1‑93 (MQ=255)

                    gaagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATgac                                                                      <  1:4203238/93‑3 (MQ=255)

                    gaagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGcc                                                                      <  2:2404273/93‑1 (MQ=255)

                    gaagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGcc                                                                      >  1:5899788/1‑93 (MQ=255)

                    gaagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGcc                                                                      <  1:3307729/93‑1 (MQ=255)

                     aagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCt                                                                     <  2:3388343/93‑1 (MQ=255)

                     aagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCt                                                                     >  2:5119270/1‑93 (MQ=255)

                     aagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCt                                                                     <  2:1275065/93‑1 (MQ=255)

                       gaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTgg                                                                   <  1:3056028/93‑1 (MQ=255)

                       gaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTgg                                                                   <  1:2263598/93‑1 (MQ=255)

                         aaTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGa                                                                 >  1:5020870/1‑93 (MQ=255)

                          aTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGat                                                                <  2:4135251/93‑1 (MQ=255)

                              ccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTc                                                            <  1:4662248/93‑1 (MQ=255)

                               cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCt                                                           >  1:3947334/1‑93 (MQ=255)

                                aTAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTa                                                          >  2:3141562/1‑93 (MQ=255)

                                aTAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTa                                                          >  2:5683508/1‑93 (MQ=255)

                                 tAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTAc                                                         <  1:1912456/93‑1 (MQ=255)

                                  aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACa                                                        <  1:6294759/93‑1 (MQ=255)

                                  aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACa                                                        >  1:1364098/1‑93 (MQ=255)

                                  aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACa                                                        <  2:4875420/93‑1 (MQ=255)

                                    cAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGa                                                      >  1:6321441/1‑93 (MQ=255)

                                     aaTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaa                                                     <  2:1561526/93‑1 (MQ=255)

                                      aTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaa                                                    <  2:767659/93‑1 (MQ=255)

                                      aTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaa                                                    >  2:5733937/1‑93 (MQ=255)

                                      aTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaa                                                    >  1:5808658/1‑93 (MQ=255)

                                       tAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaa                                                   >  1:758161/1‑93 (MQ=255)

                                        aaaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaaa                                                  >  1:2922077/1‑93 (MQ=255)

                                         aaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGaaaaaa                                                 <  2:2477478/93‑1 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                <  1:4892918/93‑1 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                <  2:6287025/93‑1 (MQ=255)

                                           tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                               <  2:1600960/93‑1 (MQ=255)

                                           tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                               <  1:5307793/93‑1 (MQ=255)

                                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                            <  1:6039051/93‑1 (MQ=255)

                                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                            <  1:1613438/93‑1 (MQ=255)

                                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                            <  1:4267204/93‑1 (MQ=255)

                                                ataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTg                                          >  2:6566277/1‑93 (MQ=255)

                                                ataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTg                                          >  2:4532445/1‑93 (MQ=255)

                                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGAc                                        >  1:2603886/1‑93 (MQ=255)

                                                  aaCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGAc                                        <  2:181230/93‑1 (MQ=255)

                                                   aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                       <  1:493138/93‑1 (MQ=255)

                                                   aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                       <  2:4488008/93‑1 (MQ=255)

                                                   aCTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                       >  1:3738172/1‑93 (MQ=255)

                                                    cTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                                      >  1:5660105/1‑93 (MQ=255)

                                                    cTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                                      >  2:121100/1‑93 (MQ=255)

                                                         tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAAAA‑AAACAGATGATTGTATGCCTGGGATATTCTACCGAAAAAACAGTCTGACATAGCGa                                 >  2:392220/1‑93 (MQ=255)

                                                         tcatGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGa                                 >  2:3007322/1‑93 (MQ=255)

                                                          catGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa                                >  1:6446128/1‑93 (MQ=255)

                                                           atGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                               <  2:473252/93‑1 (MQ=255)

                                                           atGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                               >  1:412796/1‑93 (MQ=255)

                                                           atGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                               >  2:3512532/1‑93 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                              >  2:5193166/1‑93 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                              >  2:3777299/1‑93 (MQ=255)

                                                             ggATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg                             <  2:6556308/93‑1 (MQ=255)

                                                             ggATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg                             >  1:1006087/1‑93 (MQ=255)

                                                             ggATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg                             <  1:5235922/93‑1 (MQ=255)

                                                             ggATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg                             <  2:6681244/93‑1 (MQ=255)

                                                              gATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTCCAGAAAAAACAGTCTGACATACCGAAAAgg                            <  2:3655326/93‑1 (MQ=255)

                                                              gATTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg                            <  1:3141562/93‑1 (MQ=255)

                                                               aTTTCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa                           <  2:2149552/93‑1 (MQ=255)

                                                                 ttCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGACAAAACAGTCTGACATACCGAAAAGGAtt                         >  2:2745069/1‑93 (MQ=255)

                                                                 ttCAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt                         <  2:1409393/93‑1 (MQ=255)

                                                                   cAAAAACCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaa                       <  2:2838101/93‑1 (MQ=255)

                                                                   cAAAAACCGTGAATTTCCGTGAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaa                       <  2:3272797/93‑1 (MQ=255)

                                                                    aaaaaCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTaaa                      <  1:5148870/93‑1 (MQ=255)

                                                                     aaaaCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAg                     <  2:5887411/93‑1 (MQ=255)

                                                                     aaaaCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAg                     <  1:3412955/93‑1 (MQ=255)

                                                                     aaaaCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAg                     <  1:361483/93‑1 (MQ=255)

                                                                      aaaCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGATATACCGAAAAGGATTAAAgg                    >  2:5034268/1‑93 (MQ=255)

                                                                      aaaCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAgg                    <  1:320135/93‑1 (MQ=255)

                                                                       aaCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                   >  1:4518987/1‑93 (MQ=255)

                                                                       aaCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                   <  2:6502179/93‑1 (MQ=255)

                                                                       aaCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGACTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                   >  1:1802973/1‑93 (MQ=255)

                                                                        aCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGtt                  <  1:454676/93‑1 (MQ=255)

                                                                        aCCGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGtt                  >  1:2796481/1‑93 (MQ=255)

                                                                          cGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtat                >  1:3876825/1‑93 (MQ=255)

                                                                          cGTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtat                <  1:2145377/93‑1 (MQ=255)

                                                                           gTGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtata               <  1:2223806/93‑1 (MQ=255)

                                                                            tGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTatat              <  1:5082389/93‑1 (MQ=255)

                                                                            tGAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTatat              >  1:5367147/1‑93 (MQ=255)

                                                                             gAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATAtt             <  2:4429546/93‑1 (MQ=255)

                                                                             gAATTTCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATAtt             >  2:54619/1‑93 (MQ=255)

                                                                                tttCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAc          >  1:2787620/1‑93 (MQ=255)

                                                                                  tCCGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCa        <  1:2523510/93‑1 (MQ=255)

                                                                                   ccGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAg       <  2:1070141/93‑1 (MQ=255)

                                                                                   ccGTTAG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAg       <  1:343128/93‑1 (MQ=255)

                                                                                        aG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAt  <  1:2753557/93‑1 (MQ=255)

                                                                                        aG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAt  >  1:379479/1‑93 (MQ=255)

                                                                                        aG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAt  <  1:1525709/93‑1 (MQ=255)

                                                                                        aG‑AACAAAACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAt  <  2:873816/93‑1 (MQ=255)

GAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAG‑AAC‑AAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAT  >  NC_000913/4435589‑4435768

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