BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	677,597	A→G	K320K (AAA→AAG)	ybeS →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	677,597	0	A	G	95.4%	472.9 / 9.0	173	K320K (AAA→AAG)	ybeS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/5); new base G (76/89); total (79/94)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

Only 100 of 175 total aligned reads displayed.

TGCAATGGTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAACCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCTGGGGC  >  NC_000913/677509‑677685

   aaTGGTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAACCGCTat                                                                                   >  1:3476213/1‑92 (MQ=255)

    aTGGTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACt                                                                                  <  2:3824778/93‑1 (MQ=255)

    aTGGTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACt                                                                                  >  2:1044544/1‑93 (MQ=255)

     tGGTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACtt                                                                                 <  2:4257559/93‑1 (MQ=255)

      ggTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAACCGCTACttt                                                                                <  1:4118554/93‑1 (MQ=255)

       gTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTa                                                                               <  2:3570547/93‑1 (MQ=255)

       gTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTa                                                                               >  2:3114271/1‑93 (MQ=255)

       gTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTa                                                                               <  1:4498347/93‑1 (MQ=255)

        tGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAACCGCTACTTTAt                                                                              <  1:83481/93‑1 (MQ=255)

          ccGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGt                                                                            <  2:103400/93‑1 (MQ=255)

            gTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTcc                                                                          <  2:2401577/93‑1 (MQ=255)

            gTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTcc                                                                          <  1:1497168/93‑1 (MQ=255)

            gTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTcc                                                                          >  1:4063655/1‑93 (MQ=255)

              cAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTcccc                                                                        <  1:562887/93‑1 (MQ=255)

              cAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTcccc                                                                        <  2:972009/93‑1 (MQ=255)

              cAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTcccc                                                                        <  2:1802721/93‑1 (MQ=255)

               aGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCg                                                                       <  2:3314178/93‑1 (MQ=255)

                gCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCgg                                                                      >  2:1069174/1‑93 (MQ=255)

                gCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCgg                                                                      >  1:3059483/1‑93 (MQ=255)

                 cATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGc                                                                     >  2:1019812/1‑93 (MQ=255)

                 cATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGc                                                                     <  2:1411921/93‑1 (MQ=255)

                  aTCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCg                                                                    >  1:2610379/1‑93 (MQ=255)

                   tCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGa                                                                   >  1:2169472/1‑93 (MQ=255)

                     gggTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACa                                                                 <  1:1509226/93‑1 (MQ=255)

                     gggTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACa                                                                 <  2:436788/93‑1 (MQ=255)

                      ggTCGAATATTTCCCATTACTGGGGATTGTGGATGAAGCCGCGTGTTC‑GGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAg                                                                <  2:524665/93‑1 (MQ=255)

                        tCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAgcg                                                              >  1:1138823/1‑93 (MQ=255)

                          gATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGcc                                                            >  1:1740616/1‑93 (MQ=255)

                          gATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGcc                                                            <  1:1150471/93‑1 (MQ=255)

                            ttattaCCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACCTTATGTCCCCGGCGACAGCGCCCg                                                          >  2:2458366/1‑93 (MQ=255)

                             tattaCCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCAGCTACTTTATGTCCCCGGCGACAGCGCCCGt                                                         >  1:4633203/1‑93 (MQ=255)

                              attaCCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTa                                                        <  1:681481/93‑1 (MQ=255)

                              attaCCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTa                                                        <  1:4523564/93‑1 (MQ=255)

                               ttaCCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTAc                                                       >  1:2869289/1‑93 (MQ=255)

                               ttaCCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTAc                                                       >  1:2650695/1‑93 (MQ=255)

                                 aCCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTAcgc                                                     >  1:3829111/1‑93 (MQ=255)

                                   ccATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGt                                                   >  2:2419747/1‑93 (MQ=255)

                                   ccATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGt                                                   >  2:1262904/1‑93 (MQ=255)

                                   ccATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGCACGCGt                                                   >  2:38643/1‑93 (MQ=255)

                                    cATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGtt                                                  <  2:2921034/93‑1 (MQ=255)

                                    cATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGtt                                                  >  2:904220/1‑93 (MQ=255)

                                     aTTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGttt                                                 >  2:3276287/1‑93 (MQ=255)

                                     aTTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGttt                                                 <  2:736301/93‑1 (MQ=255)

                                     aTTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGttt                                                 >  1:2242281/1‑93 (MQ=255)

                                        aCTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCa                                              <  2:3917364/93‑1 (MQ=255)

                                          tGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCAtt                                            <  1:4376544/93‑1 (MQ=255)

                                          tGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCAtt                                            >  1:1178791/1‑93 (MQ=255)

                                          tGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCAtt                                            >  1:795486/1‑93 (MQ=255)

                                          tGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCAtt                                            >  1:1620109/1‑93 (MQ=255)

                                           ggTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCAttt                                           <  1:3272157/93‑1 (MQ=255)

                                           ggTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCAttt                                           <  2:2228869/93‑1 (MQ=255)

                                           ggTGATGATAATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCAttt                                           <  2:3670557/93‑1 (MQ=255)

                                            gTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTta                                          <  1:3568367/93‑1 (MQ=255)

                                             tgatgatgaTTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTtat                                         <  1:2048996/93‑1 (MQ=255)

                                               atgatgatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATAc                                       <  1:3837991/93‑1 (MQ=255)

                                               atgatgatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATAc                                       <  2:1247702/93‑1 (MQ=255)

                                               atgatgatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATAc                                       <  1:3860441/93‑1 (MQ=255)

                                                tgatgatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACg                                      >  1:513290/1‑93 (MQ=255)

                                                tgatgatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAACCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACg                                      >  1:108414/1‑93 (MQ=255)

                                                 gatgatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAACCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGa                                     <  1:3127747/93‑1 (MQ=255)

                                                    gatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAAt                                  >  1:3316975/1‑93 (MQ=255)

                                                    gatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAAt                                  <  2:4635292/93‑1 (MQ=255)

                                                    gatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAAt                                  >  2:585507/1‑93 (MQ=255)

                                                    gatTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAAt                                  >  1:3591676/1‑93 (MQ=255)

                                                        gAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATtctc                              <  1:3429441/93‑1 (MQ=255)

                                                        gAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTc                              >  2:1546482/1‑93 (MQ=255)

                                                        gAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTc                              <  2:2520339/93‑1 (MQ=255)

                                                        gAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGACCCCGTACGCGTTTTCATTTATACGAAATACTc                              >  2:25731/1‑93 (MQ=255)

                                                           g‑cgcGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGc                           <  1:1049085/93‑1 (MQ=255)

                                                           g‑cgcGTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGc                           <  1:2431432/93‑1 (MQ=255)

                                                              gcgTGATCTGGTTAACGATCAGGGTAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGCTTTCATTTATACGAAATACACAGCGa                         >  2:1120156/1‑93 (MQ=255)

                                                               cgTGATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCCCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGAt                        <  2:4035063/93‑1 (MQ=255)

                                                                  gATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGaa                     <  1:2009001/93‑1 (MQ=255)

                                                                  gATCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGaa                     <  2:1311952/93‑1 (MQ=255)

                                                                   aTCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGaaa                    >  1:4188738/1‑93 (MQ=255)

                                                                    tCTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGaaaa                   <  1:1165089/93‑1 (MQ=255)

                                                                     cTGGTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGaaaaa                  <  1:3224816/93‑1 (MQ=255)

                                                                       ggTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAAct                >  1:1555576/1‑93 (MQ=255)

                                                                        gTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAActc               >  2:4201249/1‑93 (MQ=255)

                                                                        gTTAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAActc               <  1:3858287/93‑1 (MQ=255)

                                                                         ttAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAActct              <  1:1753315/93‑1 (MQ=255)

                                                                          tAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAActctc             <  1:1406708/93‑1 (MQ=255)

                                                                          tAACGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAActctc             <  2:1753751/93‑1 (MQ=255)

                                                                           aaCGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAActctct            <  2:247634/93‑1 (MQ=255)

                                                                           aaCGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAActctct            <  2:1443777/93‑1 (MQ=255)

                                                                           aaCGATCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAActctct            <  1:4685717/93‑1 (MQ=255)

                                                                              gATCAGGGAAAGCCGCTACTTTATGTCCCCGGCTACAGCGTCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTgcg         >  2:4680648/1‑93 (MQ=255)

                                                                               aTCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTgcgc        >  2:2336160/1‑93 (MQ=255)

                                                                                tCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCTTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCt       <  1:3221059/93‑1 (MQ=255)

                                                                                tCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGTGCt       >  1:3921576/1‑93 (MQ=255)

                                                                                tCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCt       <  1:291716/93‑1 (MQ=255)

                                                                                tCAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCt       >  2:1484920/1‑93 (MQ=255)

                                                                                 cAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTTATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCTg      <  2:2060489/93‑1 (MQ=255)

                                                                                 cAGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAACCTCTCTGCGCTg      <  1:4077541/93‑1 (MQ=255)

                                                                                  aGGGAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCTgg     <  2:725942/93‑1 (MQ=255)

                                                                                   gggAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCTggg    >  1:702231/1‑93 (MQ=255)

                                                                                   gggAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCTggg    >  1:1344482/1‑93 (MQ=255)

                                                                                   gggAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCTggg    <  1:2768658/93‑1 (MQ=255)

                                                                                   gggAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCTggg    >  1:4398828/1‑93 (MQ=255)

                                                                                    ggAAAGCCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCTgggg   <  2:827111/93‑1 (MQ=255)

TGCAATGGTGCCGTCAGCATCGGGTCGATTATTACCCATTACTGGTGATGATGATTGAAG‑CGCGTGATCTGGTTAACGATCAGGGAAAACCGCTACTTTATGTCCCCGGCGACAGCGCCCGTACGCGTTTTCATTTATACGAAATACTCAGCGATGAAAAACTCTCTGCGCTGGGGC  >  NC_000913/677509‑677685

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