BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,802	G→A	E57K (GAG→AAG)	yabB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,802	0	G	A	99.8%	1556.5 / ‑5.7	417	E57K (GAG→AAG)	yabB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (218/198); total (218/199)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 421 total aligned reads displayed.

CTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGG  >  NC_000913/89712‑89892

ctTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAg                                                                                          >  1:5288341/1‑93 (MQ=255)

ctTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAg                                                                                          >  2:5875731/1‑93 (MQ=255)

 ttGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGc                                                                                         <  1:125184/93‑1 (MQ=255)

    agaACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGCATGGGAAATTATCAAGCaaa                                                                                      >  2:5643224/1‑93 (MQ=255)

    agaACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCaaa                                                                                      <  2:1553477/93‑1 (MQ=255)

     gaACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGTAAATTATCAAGCaaaa                                                                                     >  2:6677276/1‑93 (MQ=255)

     gaACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCaaaa                                                                                     <  2:2343091/93‑1 (MQ=255)

     gaACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCaaaa                                                                                     <  1:4239096/93‑1 (MQ=255)

      aaCGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCaaaaa                                                                                    >  2:2017211/1‑93 (MQ=255)

      aaCGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCaaaaa                                                                                    <  2:6499575/93‑1 (MQ=255)

         gCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTa                                                                                 >  2:3276325/1‑93 (MQ=255)

         gCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTa                                                                                 >  2:1049227/1‑93 (MQ=255)

         ccTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTa                                                                                 <  1:8095150/92‑1 (MQ=255)

             ccGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                             <  2:5490635/93‑1 (MQ=255)

               ggTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGt                                                                           >  2:2889513/1‑93 (MQ=255)

               ggTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGt                                                                           >  1:7191503/1‑93 (MQ=255)

                gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                          >  2:4367985/1‑93 (MQ=255)

                gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                          <  2:112906/93‑1 (MQ=255)

                gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                          >  2:1867468/1‑93 (MQ=255)

                gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                          <  2:8548906/93‑1 (MQ=255)

                gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                          >  1:7102198/1‑93 (MQ=255)

                gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                          >  1:6235410/1‑93 (MQ=255)

                gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                          >  1:5696067/1‑93 (MQ=255)

                  cAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTg                                                                        >  1:5927240/1‑93 (MQ=255)

                  cAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTg                                                                        <  2:6046144/93‑1 (MQ=255)

                  cAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTg                                                                        <  2:7651186/93‑1 (MQ=255)

                   aaaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGt                                                                       >  2:3920373/1‑93 (MQ=255)

                    aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                      <  2:4600286/93‑1 (MQ=255)

                    aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                      >  2:8146364/1‑93 (MQ=255)

                      tGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGa                                                                    <  2:5743908/93‑1 (MQ=255)

                      tGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGa                                                                    <  2:7834197/93‑1 (MQ=255)

                        gTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGc                                                                  >  2:7289097/1‑93 (MQ=255)

                         tttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCa                                                                 <  1:4055297/93‑1 (MQ=255)

                         tttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCa                                                                 <  1:8176342/93‑1 (MQ=255)

                         tttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCa                                                                 <  2:6303471/93‑1 (MQ=255)

                         tttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCa                                                                 <  1:663386/93‑1 (MQ=255)

                          ttGCTCCAGTGACATTTATCGACCGTGCCGGCTGCGGTACCCCCTTCCTGAATGGGAAATTATCAAGCAAAAATTATGGCGTCTGTGGAGccg                                                                <  2:6805057/93‑3 (MQ=255)

                           tGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATg                                                               >  2:3151886/1‑93 (MQ=255)

                            gCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGa                                                              <  1:4827796/93‑1 (MQ=255)

                              aCCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAc                                                            >  1:5516389/1‑93 (MQ=255)

                               ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                           >  2:8364083/1‑93 (MQ=255)

                               ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                           <  2:4039787/93‑1 (MQ=255)

                                cATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAACATTATTGCGTCTGTCGAGCATGAAccc                                                          >  1:7504474/1‑93 (MQ=255)

                                  ttGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCgg                                                        >  1:2466361/1‑93 (MQ=255)

                                  ttGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCgg                                                        >  1:6161037/1‑93 (MQ=255)

                                   tGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGt                                                       <  1:9828637/93‑1 (MQ=255)

                                     aCATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTg                                                     <  2:8680215/93‑1 (MQ=255)

                                       aTTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAg                                                   >  1:3696143/1‑93 (MQ=255)

                                       aTTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAg                                                   >  1:4708610/1‑93 (MQ=255)

                                       aTTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGAAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAg                                                   <  2:623401/93‑1 (MQ=255)

                                        tttATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTAGAGCATGAACCCGGTTGAgc                                                  <  1:3108046/93‑1 (MQ=255)

                                         ttATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAgcg                                                 >  2:1659757/1‑93 (MQ=255)

                                            tCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCg                                              <  1:924387/93‑1 (MQ=255)

                                                ccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGtgtg                                          >  1:2863038/1‑93 (MQ=255)

                                                ccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGtgtg                                          >  2:4173495/1‑93 (MQ=255)

                                                 cGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGc                                         <  1:6859714/93‑1 (MQ=255)

                                                 cGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGc                                         >  2:7528037/1‑93 (MQ=255)

                                                 cGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGc                                         >  1:7401944/1‑93 (MQ=255)

                                                 cGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGc                                         <  1:8370987/93‑1 (MQ=255)

                                                      ctgctgCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgcgc                                    >  1:3410209/1‑93 (MQ=255)

                                                       tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                   >  1:6379218/1‑93 (MQ=255)

                                                           gctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACt                               >  1:6101001/1‑93 (MQ=255)

                                                            ctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTg                              <  1:4502894/93‑1 (MQ=255)

                                                            ctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTg                              <  1:9507560/93‑1 (MQ=255)

                                                             tttACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGt                             >  2:6948001/1‑93 (MQ=255)

                                                               tACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTa                           <  1:2201578/93‑1 (MQ=255)

                                                               tACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTa                           <  2:9830669/93‑1 (MQ=255)

                                                                aCCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAg                          <  1:1023017/93‑1 (MQ=255)

                                                                 cccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAgg                         <  1:7733472/93‑1 (MQ=255)

                                                                   cccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTc                       >  2:10042384/1‑93 (MQ=255)

                                                                   cccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTc                       >  1:5882827/1‑93 (MQ=255)

                                                                   cccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTc                       >  1:1799916/1‑93 (MQ=255)

                                                                    cctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCa                      >  1:1433074/1‑93 (MQ=255)

                                                                    cctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCa                      <  1:8900243/93‑1 (MQ=255)

                                                                        cctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGcc                  <  1:6518451/93‑1 (MQ=255)

                                                                         ctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCa                 >  2:1454436/1‑93 (MQ=255)

                                                                         ctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCa                 <  2:9026500/93‑1 (MQ=255)

                                                                           gAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGc               <  2:9600891/93‑1 (MQ=255)

                                                                            aaTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCg              >  2:8088352/1‑93 (MQ=255)

                                                                             aTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGa             <  2:377229/93‑1 (MQ=255)

                                                                             aTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGa             >  2:4982542/1‑93 (MQ=255)

                                                                              tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGaa            <  1:8437696/93‑1 (MQ=255)

                                                                               gggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAAt           <  2:9395121/93‑1 (MQ=255)

                                                                               gggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAAt           >  1:9142767/1‑93 (MQ=255)

                                                                               gggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAAt           >  2:1213729/1‑93 (MQ=255)

                                                                                ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg          >  2:3841456/1‑93 (MQ=255)

                                                                                ggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATg          <  1:9904426/93‑1 (MQ=255)

                                                                                 gAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGt         >  2:6621665/1‑93 (MQ=255)

                                                                                  aaaTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTc        <  1:6987445/93‑1 (MQ=255)

                                                                                   aaTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCa       <  1:3832069/93‑1 (MQ=255)

                                                                                    aTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAg      >  1:2022733/1‑93 (MQ=255)

                                                                                    aTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAg      >  2:4820386/1‑93 (MQ=255)

                                                                                     ttATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga     <  2:2558174/93‑1 (MQ=255)

                                                                                     ttATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga     >  2:3618547/1‑93 (MQ=255)

                                                                                     ttATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga     <  2:3526590/93‑1 (MQ=255)

                                                                                      tATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgat    >  2:7158730/1‑93 (MQ=255)

                                                                                      tATCAAGCAAAAATTAGTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgat    <  1:8687899/93‑1 (MQ=255)

                                                                                       aTCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatg   <  1:511535/93‑1 (MQ=255)

                                                                                       aTCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatg   >  2:5390480/1‑93 (MQ=255)

                                                                                       aTCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatg   <  1:1637849/93‑1 (MQ=255)

CTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGG  >  NC_000913/89712‑89892

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