BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,456,536	A→G	S83S (TCA→TCG)	b1394 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,456,536	0	A	G	99.5%	754.7 / ‑4.6	219	S83S (TCA→TCG)	b1394	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (107/111); total (108/111)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 223 total aligned reads displayed.

TACAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGC  >  NC_000913/1456446‑1456625

tacAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGt                                                                                         <  2:1200158/93‑1 (MQ=255)

   aGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGa                                                                                      >  2:5899271/1‑93 (MQ=255)

    ggTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGaa                                                                                     >  1:1020850/1‑93 (MQ=255)

      tGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACg                                                                                   <  1:3281753/93‑1 (MQ=255)

       gCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGt                                                                                  <  1:3838342/93‑1 (MQ=255)

        ccGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGtt                                                                                 <  1:689012/93‑1 (MQ=255)

          gggCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTc                                                                               <  2:2820329/93‑1 (MQ=255)

          gggCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTc                                                                               <  1:617634/93‑1 (MQ=255)

           ggCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCt                                                                              >  1:6223174/1‑93 (MQ=255)

           ggCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCt                                                                              <  2:2283755/93‑1 (MQ=255)

            gcgcGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCta                                                                             <  1:5007777/93‑1 (MQ=255)

              gcgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGCTTTAGGGATGTCGGTTGAACGTTTCtata                                                                           >  1:1572534/1‑93 (MQ=255)

               cgGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATaa                                                                          <  2:185574/93‑1 (MQ=255)

                 gATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAAcc                                                                        <  2:2110652/93‑1 (MQ=255)

                  aTTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAAccc                                                                       >  2:4314038/1‑93 (MQ=255)

                  aTTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAAccc                                                                       >  1:2880774/1‑93 (MQ=255)

                  aTTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAAccc                                                                       <  2:516937/93‑1 (MQ=255)

                   ttttGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCa                                                                      >  1:1720746/1‑93 (MQ=255)

                      tgtgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTg                                                                   <  1:1092980/93‑1 (MQ=255)

                       gtgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTgg                                                                  >  1:3135810/1‑93 (MQ=255)

                       gtgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTgg                                                                  <  1:5396260/93‑1 (MQ=255)

                        tgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGt                                                                 <  2:1518200/93‑1 (MQ=255)

                        tgCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGt                                                                 <  1:5282317/93‑1 (MQ=255)

                          cTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTAc                                                               <  1:5547565/93‑1 (MQ=255)

                          cTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTAc                                                               <  2:3060927/93‑1 (MQ=255)

                          cTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTAc                                                               >  2:3831255/1‑93 (MQ=255)

                           tGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACg                                                              >  1:4269351/1‑93 (MQ=255)

                            ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTTGTACGt                                                             <  1:5362265/93‑1 (MQ=255)

                            ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGt                                                             <  2:2648257/93‑1 (MQ=255)

                            ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGt                                                             <  1:6287047/93‑1 (MQ=255)

                            ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGt                                                             <  2:538690/93‑1 (MQ=255)

                            ggTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGt                                                             >  2:344850/1‑93 (MQ=255)

                             gTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTc                                                            <  2:1275769/93‑1 (MQ=255)

                             gTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTc                                                            >  1:704377/1‑93 (MQ=255)

                             gTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGCCGGTTGAACGTTTCTATAACCCACTGGTACGTc                                                            <  2:1754553/93‑1 (MQ=255)

                                aGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGcc                                                         >  2:4805429/1‑93 (MQ=255)

                                aGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGcc                                                         >  1:4817388/1‑93 (MQ=255)

                                 ggATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCt                                                        <  2:2272208/93‑1 (MQ=255)

                                 ggATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCt                                                        >  2:1064626/1‑93 (MQ=255)

                                    tCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGGTGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGc                                                     >  2:3274170/1‑93 (MQ=255)

                                        aaCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCaaaa                                                 >  2:3033427/1‑93 (MQ=255)

                                         aCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGTATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAAc                                                <  1:2649831/93‑1 (MQ=255)

                                         aCGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAAc                                                >  1:2715366/1‑93 (MQ=255)

                                           gATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTg                                              >  2:304634/1‑93 (MQ=255)

                                            aTCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGc                                             <  2:4632272/93‑1 (MQ=255)

                                            aTCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGc                                             >  1:4145740/1‑93 (MQ=255)

                                             tCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGcc                                            <  1:2336670/93‑1 (MQ=255)

                                              cGTAACGTTGATCCCACATGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCa                                           <  2:832701/93‑1 (MQ=255)

                                                 aaCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCaaaa                                        >  2:170673/1‑93 (MQ=255)

                                                  aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAc                                       >  2:3294987/1‑93 (MQ=255)

                                                  aCGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAc                                       >  1:2191854/1‑93 (MQ=255)

                                                     ttGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCgg                                    <  1:5986347/93‑1 (MQ=255)

                                                      tGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGt                                   <  2:5052015/93‑1 (MQ=255)

                                                       gATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTg                                  <  1:236467/93‑1 (MQ=255)

                                                       gATCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTg                                  <  1:2140338/93‑1 (MQ=255)

                                                        aTCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGa                                 >  2:6119863/1‑93 (MQ=255)

                                                        aTCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGa                                 <  1:1999877/93‑1 (MQ=255)

                                                        aTCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGa                                 >  1:246418/1‑93 (MQ=255)

                                                        aTCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGa                                 >  1:4500342/1‑93 (MQ=255)

                                                        aTCCCACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGa                                 >  1:3361806/1‑93 (MQ=255)

                                                          cccACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATc                               >  1:6377684/1‑93 (MQ=255)

                                                           ccACAGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCt                              <  2:3123306/93‑1 (MQ=255)

                                                            cacaGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtg                             >  1:6410691/1‑93 (MQ=255)

                                                            cacaGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtg                             <  2:4969651/93‑1 (MQ=255)

                                                              caGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCtgtg                           <  2:5754874/93‑1 (MQ=255)

                                                               aGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGc                          >  2:952291/1‑93 (MQ=255)

                                                               aGGCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGc                          <  1:344850/93‑1 (MQ=255)

                                                               aGGCCCCGCACCGGATTTAGGGATGTCGGGTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGc                          >  1:4019615/1‑93 (MQ=255)

                                                                ggCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCa                         <  2:2465479/93‑1 (MQ=255)

                                                                ggCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCa                         >  1:2534853/1‑93 (MQ=255)

                                                                ggCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCa                         <  1:6094016/93‑1 (MQ=255)

                                                                ggCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCa                         >  2:3899039/1‑93 (MQ=255)

                                                                 gCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAg                        <  1:2428989/93‑1 (MQ=255)

                                                                 gCCCCGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAg                        <  1:3323136/93‑1 (MQ=255)

                                                                   cccGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTc                      >  1:1918320/1‑93 (MQ=255)

                                                                    ccGCACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCa                     <  2:1974513/93‑1 (MQ=255)

                                                                       cACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtg                  <  2:5799307/93‑1 (MQ=255)

                                                                       cACCGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtg                  <  1:6291822/93‑1 (MQ=255)

                                                                         ccGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggc                <  1:2907725/93‑1 (MQ=255)

                                                                         ccGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggc                >  1:3180201/1‑93 (MQ=255)

                                                                         ccGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggc                >  2:4914201/1‑93 (MQ=255)

                                                                          cGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcg               >  2:2008154/1‑93 (MQ=255)

                                                                          cGGATTTAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcg               >  2:4185343/1‑93 (MQ=255)

                                                                              tttAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgtggc           >  1:1204741/1‑93 (MQ=255)

                                                                              tttAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgtggc           <  2:2926366/93‑1 (MQ=255)

                                                                              tttAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAAtggcgtggc           <  2:5543879/93‑1 (MQ=255)

                                                                              tttAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAAACGGTGATCTGTGCAGTCAAtggcgtggc           >  2:286795/1‑93 (MQ=255)

                                                                                tAGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTggcgg         >  1:864665/1‑93 (MQ=255)

                                                                                 aGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTggcggc        <  2:6278575/93‑1 (MQ=255)

                                                                                 aGGGATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTggcggc        >  1:4222598/1‑93 (MQ=255)

                                                                                  gggATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggca       >  2:5691823/1‑93 (MQ=255)

                                                                                  gggATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggca       <  2:1926072/93‑1 (MQ=255)

                                                                                  gggATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggca       >  1:3800974/1‑93 (MQ=255)

                                                                                    gATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagg     <  2:5504501/93‑1 (MQ=255)

                                                                                    gATGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagg     >  1:4554764/1‑93 (MQ=255)

                                                                                     aTGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggg    >  2:1580687/1‑93 (MQ=255)

                                                                                      tGTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcagggg   <  2:3591796/93‑1 (MQ=255)

                                                                                       gTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggggc  <  1:3580123/93‑1 (MQ=255)

                                                                                       gTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggggc  >  1:5824539/1‑93 (MQ=255)

                                                                                       gTCGGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCggcaggggc  >  2:1026872/1‑93 (MQ=255)

TACAGGTGCCGGGCGCGGATTTTGTGCTGGTCAGGATCTTAACGATCGTAACGTTGATCCCACAGGCCCCGCACCGGATTTAGGGATGTCAGTTGAACGTTTCTATAACCCACTGGTACGTCGCCTGGCAAAACTGCCAAAACCGGTGATCTGTGCAGTCAATGGCGTGGCGGCAGGGGC  >  NC_000913/1456446‑1456625

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