BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	804,347	T→C	V549A (GTG→GCG)	ybhJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	804,347	0	T	C	99.6%	719.1 / NA	222	V549A (GTG→GCG)	ybhJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (111/110); minor base G (1/0); total (112/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

AAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTC  >  NC_000913/804257‑804436

aaGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAgcgt                                                                                         <  1:2908872/93‑1 (MQ=255)

  ggggCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGc                                                                                       <  1:4125699/93‑1 (MQ=255)

    ggCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTc                                                                                     <  1:3952458/93‑1 (MQ=255)

     gCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCg                                                                                    >  2:6082698/1‑93 (MQ=255)

       aaCTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGaa                                                                                  >  2:1043715/1‑93 (MQ=255)

       aaCTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGaa                                                                                  <  1:3854838/93‑1 (MQ=255)

       aaCTCAGCAACCGCTGATTAACGGGCCGATGATGAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGaa                                                                                  <  2:5641034/93‑1 (MQ=255)

         cTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGa                                                                                <  2:3159094/93‑1 (MQ=255)

          tCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGAt                                                                               <  2:4629894/93‑1 (MQ=255)

            aGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATcc                                                                             <  1:3833063/93‑1 (MQ=255)

            aGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATcc                                                                             >  2:3901423/1‑93 (MQ=255)

            aGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATcc                                                                             >  1:2120455/1‑93 (MQ=255)

            aGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATcc                                                                             <  1:5400978/93‑1 (MQ=255)

             gCAACCGCTGATTTACGGGCCGAACATTTAAGACTGGCCGGGATTGGGTGCGCTGACTGACCATATCGTCGTGAACGCGTGCTCGAAGATACt                                                                            >  2:350314/1‑93 (MQ=255)

             gCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCt                                                                            >  2:2171742/1‑93 (MQ=255)

             gCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCt                                                                            >  2:1506735/1‑93 (MQ=255)

             gCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCt                                                                            <  1:4842587/93‑1 (MQ=255)

             gCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCt                                                                            <  1:496250/93‑1 (MQ=255)

             gCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCt                                                                            <  2:847525/93‑1 (MQ=255)

             gCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCt                                                                            >  1:189248/1‑93 (MQ=255)

              cAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTc                                                                           <  1:6064979/93‑1 (MQ=255)

                 ccGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTcgac                                                                        <  2:1949818/93‑1 (MQ=255)

                  cGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTcgacg                                                                       >  1:5116044/1‑93 (MQ=255)

                    cTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGaa                                                                     >  1:311133/1‑93 (MQ=255)

                    cTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGaa                                                                     >  1:1369564/1‑93 (MQ=255)

                     tGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAg                                                                    <  1:963082/93‑1 (MQ=255)

                     tGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAg                                                                    <  1:4194391/93‑1 (MQ=255)

                      gATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGt                                                                   >  2:658053/1‑93 (MQ=255)

                      gATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGt                                                                   >  1:3048573/1‑93 (MQ=255)

                       aTTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTg                                                                  <  2:1389743/93‑1 (MQ=255)

                       aTTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTg                                                                  <  2:2867465/93‑1 (MQ=255)

                        tttACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGa                                                                 >  2:4273176/1‑93 (MQ=255)

                         ttACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGac                                                                <  1:468825/93‑1 (MQ=255)

                           aCGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGGGTGCTCGAAGATACTCGACGAAGTGacca                                                              >  1:148787/1‑93 (MQ=255)

                            cGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGaccac                                                             <  2:417592/93‑1 (MQ=255)

                               gCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGaccaccac                                                          <  1:3533581/93‑1 (MQ=255)

                                ccGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGaccaccacc                                                         <  1:5437874/93‑1 (MQ=255)

                                 cGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGAACACCACcg                                                        >  1:1183364/1‑93 (MQ=255)

                                  gAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACcga                                                       <  1:84525/93‑1 (MQ=255)

                                   aaCATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACcgac                                                      <  2:1821037/93‑1 (MQ=255)

                                   aaCATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACcgac                                                      <  2:1873440/93‑1 (MQ=255)

                                     cATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACcgacga                                                    >  2:4798449/1‑93 (MQ=255)

                                     cATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACcgacga                                                    <  2:4548509/93‑1 (MQ=255)

                                     cATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACcgacga                                                    >  1:5548036/1‑93 (MQ=255)

                                      aTTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGaa                                                   <  1:3660688/93‑1 (MQ=255)

                                        tAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACt                                                 <  2:1800905/93‑1 (MQ=255)

                                        tAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACt                                                 >  1:5397657/1‑93 (MQ=255)

                                          aaGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGa                                               <  2:3101000/93‑1 (MQ=255)

                                          aaGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGa                                               >  2:4778030/1‑93 (MQ=255)

                                            gccTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGAtt                                             >  2:3295847/3‑93 (MQ=255)

                                            gACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGAtt                                             >  2:1781384/1‑93 (MQ=255)

                                             aCTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGAttc                                            <  2:2210899/93‑1 (MQ=255)

                                               tGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGAttgct                                          >  2:2961102/1‑93 (MQ=255)

                                                ggCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGAttcctt                                         >  1:626758/1‑93 (MQ=255)

                                                 gCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCAACGAAGTGACCACCACCGACGAACTGAttccttc                                        <  2:3783907/93‑1 (MQ=255)

                                                  ccGGAATTGGGTGAGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGAttccttcc                                       <  1:3604059/93‑1 (MQ=255)

                                                     gAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGt                                    >  1:6004405/1‑93 (MQ=255)

                                                     gAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCAACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGt                                    >  2:4445666/1‑93 (MQ=255)

                                                      aaTTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTg                                   >  1:2506507/1‑93 (MQ=255)

                                                      aaTTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGAACACCACCGACGAACTGATTCCTTCCGGTg                                   >  2:2099552/1‑93 (MQ=255)

                                                       aTTGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGa                                  >  1:5358816/1‑93 (MQ=255)

                                                        ttGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGaa                                 <  1:2023034/93‑1 (MQ=255)

                                                        ttGGGTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGaa                                 >  2:3235319/1‑93 (MQ=255)

                                                           ggTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAAcc                              <  2:5817155/93‑1 (MQ=255)

                                                           ggTGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGACCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAAcc                              <  2:6213194/93‑1 (MQ=255)

                                                             tGCGCTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctc                            <  2:1106681/93‑1 (MQ=255)

                                                              gcgcTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctct                           <  1:4170134/93‑1 (MQ=255)

                                                              gcgcTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctct                           <  2:3967218/93‑1 (MQ=255)

                                                              gcgcTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctct                           <  1:3691155/93‑1 (MQ=255)

                                                               cgcTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCtctt                          <  2:6050797/93‑1 (MQ=255)

                                                               cgcTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCtctt                          <  2:1612328/93‑1 (MQ=255)

                                                               cgcTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCtctt                          <  1:4814135/93‑1 (MQ=255)

                                                                gcTGACTGACAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGCCCACCACCGACGAACTGATTCCTTCCGGTGAAACCtcttc                         <  1:2600217/93‑1 (MQ=255)

                                                                 ctgactgaCAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCtcttct                        >  2:5559712/1‑93 (MQ=255)

                                                                  tgactgacAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTcttctt                       <  1:3337086/93‑1 (MQ=255)

                                                                  tgactgacAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTcttctt                       <  2:2368534/93‑1 (MQ=255)

                                                                  tgactgacAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTcttctt                       >  2:3672908/1‑93 (MQ=255)

                                                                   gactgacAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTa                      >  2:1431943/1‑93 (MQ=255)

                                                                       gacAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGt                  <  1:4223538/93‑1 (MQ=255)

                                                                       gacAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGt                  >  2:760434/1‑93 (MQ=255)

                                                                        acAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGtt                 <  2:3055779/93‑1 (MQ=255)

                                                                         cAATATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTc                >  1:5929397/1‑93 (MQ=255)

                                                                          aaTATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCa               <  1:1020240/93‑1 (MQ=255)

                                                                          aaTATCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCa               >  1:5707699/1‑93 (MQ=255)

                                                                           atatCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCaa              >  1:5627766/1‑93 (MQ=255)

                                                                             atCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAAt            <  2:4327147/93‑1 (MQ=255)

                                                                             atCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAAt            >  2:5435566/1‑93 (MQ=255)

                                                                             atCGTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAAt            <  2:2697578/93‑1 (MQ=255)

                                                                              tCGTCCTGAAAGCGTGCTCGAAGATCCGCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATc           <  2:2577493/93‑1 (MQ=255)

                                                                                gTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCg         >  1:2422271/1‑93 (MQ=255)

                                                                                gTCCTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCg         >  1:67076/1‑93 (MQ=255)

                                                                                   cTGAAAGGGTGCCCCAAGAAGCTCGACGAAGTGACCCCCACCCACGAACTGATTCCTTCCGGTGAAGCCACTTCTTATCGTTCAAATCCGGtt      <  2:5073669/93‑1 (MQ=255)

                                                                                   cTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGAtt      >  2:3461915/1‑93 (MQ=255)

                                                                                   cTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGAtt      >  1:4110319/1‑93 (MQ=255)

                                                                                   cTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGAtt      <  2:5727321/93‑1 (MQ=255)

                                                                                   cTGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGAtt      <  1:4591761/93‑1 (MQ=255)

                                                                                    tGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTg     >  2:5953258/1‑93 (MQ=255)

                                                                                    tGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTg     >  2:2531387/1‑93 (MQ=255)

                                                                                    tGAAAGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTg     >  1:700682/1‑93 (MQ=255)

                                                                                      aaaGCGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGt   >  2:5031026/1‑93 (MQ=255)

AAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTC  >  NC_000913/804257‑804436

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