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	99.3%	993.7 / ‑1.8	280	coding (397/501 nt)	ytfI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/0); new base A (140/138); total (142/138)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 305 total aligned reads displayed.

GAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑AA‑ACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAT  >  NC_000913/4435589‑4435768

 aGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑a                                                                                         <  1:2419753/93‑1 (MQ=255)

   gTAGTGGATTTGAATACTAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aa‑a                                                                                      <  1:2536172/93‑1 (MQ=255)

   gTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aa‑a                                                                                      <  2:4496812/93‑1 (MQ=255)

   gTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aa‑a                                                                                      <  1:4498864/93‑1 (MQ=255)

     aGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aa‑aac                                                                                    <  2:3275207/93‑3 (MQ=255)

      gTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aa‑aaca                                                                                   >  2:4731605/1‑90 (MQ=255)

      gTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aa‑aaca                                                                                   <  1:2769378/93‑4 (MQ=255)

       tGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGaac‑aa‑aacag                                                                                  <  2:899337/93‑5 (MQ=255)

         gATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGAc                                                                                 >  2:4102575/1‑93 (MQ=255)

              gAATACGAAGAAATACCCATAACAATAAATGATTATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTTGAACAAA‑ACAGATGATTg                                                                            <  2:264270/93‑1 (MQ=255)

               aaTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGt                                                                           >  2:1464567/1‑93 (MQ=255)

               aaTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGt                                                                           <  1:1539087/93‑1 (MQ=255)

               aaTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGt                                                                           <  2:4149082/93‑1 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTa                                                                          >  2:1477203/1‑93 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTa                                                                          >  2:1210361/1‑93 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTa                                                                          >  1:588522/1‑93 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTa                                                                          >  2:3058136/1‑93 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTa                                                                          <  1:3911904/93‑1 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTa                                                                          >  2:3492503/1‑93 (MQ=255)

                aTACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTa                                                                          <  1:2659982/93‑1 (MQ=255)

                 tACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTAt                                                                         >  1:1065157/1‑93 (MQ=255)

                  aCGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATg                                                                        >  1:1024379/1‑93 (MQ=255)

                   cGAAGAAATACCCATAACAATAAATGAATATAAGTTATTCATGGATTTCAAAAACCGTGAATGTCCGTTAGAACAAA‑ACAGATGATTGTATGc                                                                       >  2:2530866/1‑93 (MQ=255)

                    gaagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAGCCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGcc                                                                      <  1:2925631/93‑1 (MQ=255)

                    gaagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGcc                                                                      >  1:4135852/1‑93 (MQ=255)

                     aagaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCt                                                                     <  1:372589/93‑1 (MQ=255)

                      agaaATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTg                                                                    >  2:800162/1‑93 (MQ=255)

                      agaaATACCCATAAAAATAAATGAATATAACTCATTCATGGATTAAAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTg                                                                    <  2:2432990/93‑1 (MQ=255)

                         aaTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGa                                                                 <  2:3667927/93‑1 (MQ=255)

                          aTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGat                                                                <  1:2963288/93‑1 (MQ=255)

                          aTACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGat                                                                <  2:1644964/93‑1 (MQ=255)

                           tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGata                                                               >  2:2152193/1‑93 (MQ=255)

                           tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGata                                                               >  1:203523/1‑93 (MQ=255)

                           tACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGata                                                               <  2:1005204/93‑1 (MQ=255)

                           tACCCATAACAATAAATGAATATAACTCACTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGata                                                               <  1:3759301/93‑1 (MQ=255)

                              ccATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTc                                                            <  1:938577/93‑1 (MQ=255)

                               cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCt                                                           >  2:4336320/1‑93 (MQ=255)

                               cATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCt                                                           >  1:1553960/1‑93 (MQ=255)

                                 tAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTAc                                                         >  1:266978/1‑93 (MQ=255)

                                  aaCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACa                                                        <  1:4262028/93‑1 (MQ=255)

                                   aCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAg                                                       >  1:2912969/1‑93 (MQ=255)

                                   aCAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑AACACAGATGATTGCATGCCTGGGATATTCTACAg                                                       >  2:652956/1‑93 (MQ=255)

                                     aaTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGaa                                                     >  2:2860062/1‑93 (MQ=255)

                                      aTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGaaa                                                    >  1:4596215/1‑93 (MQ=255)

                                      aTAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGaaa                                                    <  2:2657953/93‑1 (MQ=255)

                                         aaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACATaaaaaa                                                 >  1:2320990/1‑93 (MQ=255)

                                         aaTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGaaaaaa                                                 <  1:763975/93‑1 (MQ=255)

                                          aTGAATATAACTCCTTCATGGATTTCAACAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTCCAGacaaacc                                                <  1:4209670/93‑3 (MQ=255)

                                          aTGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAAc                                                <  2:4007754/93‑1 (MQ=255)

                                           tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                               <  1:3116846/93‑1 (MQ=255)

                                           tGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACa                                               >  1:2231563/1‑93 (MQ=255)

                                            gAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAg                                              >  1:3381091/1‑93 (MQ=255)

                                             aaTATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGt                                             >  1:6528/1‑93 (MQ=255)

                                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                            <  2:1957425/93‑1 (MQ=255)

                                              atatAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTc                                            <  2:2952932/93‑1 (MQ=255)

                                               tataACTCATTCATGTATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                           <  2:3168704/93‑1 (MQ=255)

                                               tataACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCt                                           <  1:117293/93‑1 (MQ=255)

                                                   aCTCATTCATGGATTTCAAAAACCGTGAATGTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACa                                       <  1:3818843/93‑1 (MQ=255)

                                                    cTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                                      <  2:1390315/93‑1 (MQ=255)

                                                    cTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACAt                                      >  1:4013756/1‑93 (MQ=255)

                                                     tcattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATa                                     <  1:2279070/93‑1 (MQ=255)

                                                      cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAc                                    <  2:1769168/93‑1 (MQ=255)

                                                      cattcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAc                                    >  1:874295/1‑93 (MQ=255)

                                                       attcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATAcc                                   >  2:4110201/1‑93 (MQ=255)

                                                        ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATActg                                  >  2:3793587/1‑91 (MQ=255)

                                                        ttcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCg                                  <  1:3625840/93‑1 (MQ=255)

                                                         tcatGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGa                                 >  2:613093/1‑93 (MQ=255)

                                                          catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa                                >  1:3128262/1‑93 (MQ=255)

                                                          catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa                                >  2:1930235/1‑93 (MQ=255)

                                                          catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa                                <  1:1288737/93‑1 (MQ=255)

                                                          catGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa                                >  2:4642403/1‑93 (MQ=255)

                                                           atGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                               >  2:283349/1‑93 (MQ=255)

                                                           atGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                               >  1:1851286/1‑93 (MQ=255)

                                                           atGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa                               >  1:2566592/1‑93 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                              <  1:1726315/93‑1 (MQ=255)

                                                            tGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa                              <  2:2240655/93‑1 (MQ=255)

                                                               aTTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa                           <  2:1261515/93‑1 (MQ=255)

                                                               aTTTCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa                           >  2:665220/1‑93 (MQ=255)

                                                                 ttCAAAAACCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt                         >  1:4668931/1‑93 (MQ=255)

                                                                       aaCCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                   >  2:1139789/1‑93 (MQ=255)

                                                                       aaCCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                   <  2:2393929/93‑1 (MQ=255)

                                                                       aaCCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                   <  1:3303373/93‑1 (MQ=255)

                                                                       aaCCGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGACTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGt                   <  1:395665/93‑1 (MQ=255)

                                                                          cGTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtat                >  1:1629722/1‑93 (MQ=255)

                                                                           gTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtata               >  2:2027223/1‑93 (MQ=255)

                                                                           gTGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTtata               >  2:4318362/1‑93 (MQ=255)

                                                                            tGAATTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTatat              <  2:2470030/93‑1 (MQ=255)

                                                                              aaTTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATCCTCCAGAAAAACCAGCTTGCAATACCGAAAGGGTTTAAAGGTAATAtta            >  2:3097725/1‑92 (MQ=255)

                                                                               aTTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCa           <  2:3549529/93‑1 (MQ=255)

                                                                               aTTTCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCa           <  1:3019819/93‑1 (MQ=255)

                                                                                tttCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAc          <  1:2391765/93‑1 (MQ=255)

                                                                                tttCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAc          <  1:3434243/93‑1 (MQ=255)

                                                                                 ttCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAcc         >  1:1687944/1‑93 (MQ=255)

                                                                                 ttCCGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCAcc         <  2:1525772/93‑1 (MQ=255)

                                                                                   ccGTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAg       <  2:771510/93‑1 (MQ=255)

                                                                                     gTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGc     >  2:4714158/1‑93 (MQ=255)

                                                                                     gTTAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGc     <  1:2092031/93‑1 (MQ=255)

                                                                                      ttAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCg    >  1:1422722/1‑93 (MQ=255)

                                                                                      ttAGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCg    >  2:3378312/1‑93 (MQ=255)

                                                                                        aGAACAAA‑ACAGATGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAt  >  1:4739355/1‑93 (MQ=255)

GAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAAC‑AA‑ACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTAAAGGTTATATTCACCAGGCGAT  >  NC_000913/4435589‑4435768

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