BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,311,496	+C	coding (301/408 nt)	phnQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,311,495	1	.	C	99.5%	1281.8 / ‑1.9	370	coding (302/408 nt)	phnQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (193/175); total (193/177)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 405 total aligned reads displayed.

GAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑C‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311403‑4311582

   tAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑ccg                                                                                       >  1:636675/1‑91 (MQ=255)

     aTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑ccgtt                                                                                     >  2:6454515/1‑89 (MQ=255)

     aTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑ccgtt                                                                                     >  2:7427377/1‑89 (MQ=255)

     aTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑ccgtt                                                                                     >  2:8098702/1‑89 (MQ=255)

       gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTcc                                                                                   >  2:8444768/1‑93 (MQ=255)

        aTTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCg                                                                                  >  1:946573/1‑93 (MQ=255)

         ttAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTACTGAAACGCTTTCCGCTATTTCTTTTATTAACCCTGCTAACGCTGCCCACC‑GTTCCgc                                                                                 >  1:1656156/1‑92 (MQ=255)

         ttAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCgg                                                                                 <  2:2099050/93‑1 (MQ=255)

           aaTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGcc                                                                               >  1:5359097/1‑93 (MQ=255)

            aTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCCCCCTGCTCACGCTGTCCACC‑GTTCCGGccc                                                                              <  1:5212620/93‑1 (MQ=255)

             tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCa                                                                             <  2:7707022/93‑1 (MQ=255)

              gtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAg                                                                            >  1:4404354/1‑93 (MQ=255)

              gtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAg                                                                            <  2:1036308/93‑1 (MQ=255)

                 gATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGcc                                                                         >  2:82726/1‑93 (MQ=255)

                  aTAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCg                                                                        <  2:1599813/93‑1 (MQ=255)

                  aTAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCg                                                                        <  1:2042198/93‑1 (MQ=255)

                   tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGa                                                                       <  1:7684885/93‑1 (MQ=255)

                     gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGata                                                                     <  2:351615/93‑1 (MQ=255)

                      agTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGatat                                                                    <  2:4982549/93‑1 (MQ=255)

                      agTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGatat                                                                    >  1:7552908/1‑93 (MQ=255)

                      agTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGatat                                                                    <  2:4893058/93‑1 (MQ=255)

                       gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGAtatt                                                                   >  1:6775724/1‑93 (MQ=255)

                        tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGAtatta                                                                  <  2:6888659/93‑1 (MQ=255)

                          tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATc                                                                >  1:3820459/1‑93 (MQ=255)

                          tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATc                                                                >  2:5198500/1‑93 (MQ=255)

                          tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATc                                                                >  2:4274789/1‑93 (MQ=255)

                            tttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCTAGGCCGATATTATcga                                                              <  2:7253235/93‑1 (MQ=255)

                            tttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATcga                                                              >  1:7767129/1‑93 (MQ=255)

                               ttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATcgacga                                                           >  2:3013772/1‑93 (MQ=255)

                                tGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCgacgac                                                          >  1:3539823/1‑93 (MQ=255)

                                tGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTACCgacgac                                                          <  2:5323467/93‑1 (MQ=255)

                                  aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGacac                                                        <  2:3459797/93‑1 (MQ=255)

                                    gTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGACCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCa                                                      >  2:2877821/1‑93 (MQ=255)

                                    gTTGACTATAGATACGATTTCCGCTAATTCTTTTATACACCATGCTCACGCTGTCCACC‑GTTCCGGCACAGACCGGAATTATCGACGaaacat                                                      <  2:7844641/93‑3 (MQ=255)

                                     ttGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCAt                                                     <  2:5446266/93‑1 (MQ=255)

                                      tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                    >  2:1787798/1‑93 (MQ=255)

                                      tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                    <  1:802369/93‑1 (MQ=255)

                                       ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTACCGACGACACCATTg                                                   <  2:8665352/93‑1 (MQ=255)

                                        gCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                                  >  2:7554066/1‑93 (MQ=255)

                                         cTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGc                                                 >  2:3674452/1‑93 (MQ=255)

                                          tAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCa                                                >  1:3316290/1‑93 (MQ=255)

                                          tAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCa                                                >  2:1988869/1‑93 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               >  1:6151972/1‑93 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               <  2:3579298/93‑1 (MQ=255)

                                              gAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                                            <  2:8534644/93‑1 (MQ=255)

                                              gAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                                            <  1:5592304/93‑1 (MQ=255)

                                              gAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGAAatat                                            <  2:1250683/93‑1 (MQ=255)

                                                aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTc                                          <  1:7071322/93‑1 (MQ=255)

                                                 aCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTca                                         <  1:3928417/93‑1 (MQ=255)

                                                 aCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTca                                         <  1:4244963/93‑1 (MQ=255)

                                                  cGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcag                                        <  2:4448356/93‑1 (MQ=255)

                                                   gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                       <  1:624595/93‑1 (MQ=255)

                                                   gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                       >  1:8149158/1‑93 (MQ=255)

                                                    cTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagca                                      <  1:1609871/93‑1 (MQ=255)

                                                     tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                     >  1:4607516/1‑93 (MQ=255)

                                                     tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                     <  2:5863535/93‑1 (MQ=255)

                                                     tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                     >  2:7099457/1‑93 (MQ=255)

                                                       tCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCAATGGCAATATTcagcagcc                                   >  1:4855395/1‑91 (MQ=255)

                                                         cGCTATTTCTTTTATTCCCCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                                 <  1:226979/93‑1 (MQ=255)

                                                         cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                                 <  2:2638226/93‑1 (MQ=255)

                                                         cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTAAGCAGGCGa                                 >  1:2666065/1‑93 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                >  1:1905069/1‑93 (MQ=255)

                                                           ttATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAtt                               <  2:2522975/92‑1 (MQ=255)

                                                           cTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAtt                               >  1:1010414/1‑93 (MQ=255)

                                                           cTATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAtt                               <  2:6670411/93‑1 (MQ=255)

                                                            tATTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTa                              >  2:6779568/1‑93 (MQ=255)

                                                             aTTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                             >  2:6534178/1‑93 (MQ=255)

                                                             aTTTCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                             >  1:4465377/1‑93 (MQ=255)

                                                              tttCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAAt                            <  2:3899210/93‑1 (MQ=255)

                                                               ttCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGCTTAATg                           >  1:4085231/1‑93 (MQ=255)

                                                                tCTTTTATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa                          >  2:7977649/1‑93 (MQ=255)

                                                                 cTTTTATTCCCCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGAc                         <  1:6615244/93‑1 (MQ=255)

                                                                  ttttATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACg                        >  2:1360113/1‑93 (MQ=255)

                                                                   tttATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGc                       <  2:2347013/93‑1 (MQ=255)

                                                                    ttATTCCCCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                      <  1:7944795/93‑1 (MQ=255)

                                                                    ttATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                      <  1:8600762/93‑1 (MQ=255)

                                                                     tATTCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCtt                     >  2:1658977/1‑93 (MQ=255)

                                                                      tttCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTta                    <  1:7392275/92‑1 (MQ=255)

                                                                       ttCTCCCTGCTCCCGCTGTCCA‑CTGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACTCTtat                   <  2:3818942/93‑1 (MQ=255)

                                                                       ttCACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat                   >  1:5016447/1‑93 (MQ=255)

                                                                         cACCCTGCTCACGCTGTCGACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCCGCAGGCGATTAATGACGCTTATaa                 >  1:146774/1‑93 (MQ=255)

                                                                         cACCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                 <  1:2772832/93‑1 (MQ=255)

                                                                          aCCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                <  1:5985843/93‑1 (MQ=255)

                                                                          aCCCTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGAAGGCGATTAATGACGCTTATAAc                >  1:1787837/1‑93 (MQ=255)

                                                                           cccTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAcc               >  1:1198444/1‑93 (MQ=255)

                                                                           cccTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAcc               >  1:241942/1‑93 (MQ=255)

                                                                            ccTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc              <  2:6160539/93‑1 (MQ=255)

                                                                            ccTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc              >  1:4729499/1‑93 (MQ=255)

                                                                            ccTGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc              >  1:683290/1‑93 (MQ=255)

                                                                              tGCTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTg            <  1:5689788/93‑1 (MQ=255)

                                                                                cTCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          <  2:3136976/93‑1 (MQ=255)

                                                                                 tCACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAcc         >  1:2522936/1‑93 (MQ=255)

                                                                                  cACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCg        >  2:3285661/1‑93 (MQ=255)

                                                                                  cACGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCg        <  2:6217693/93‑1 (MQ=255)

                                                                                   aCGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGt       >  2:5648571/1‑93 (MQ=255)

                                                                                   aCGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGt       >  1:4945068/1‑93 (MQ=255)

                                                                                   aCGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGt       <  2:4520542/93‑1 (MQ=255)

                                                                                    cGCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAGTAATGACGCTTATAACCCTGACCGTg      <  2:1528469/93‑1 (MQ=255)

                                                                                     gCTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTgg     <  1:3070705/93‑1 (MQ=255)

                                                                                      cTGTCCACC‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTggg    <  2:2246481/93‑1 (MQ=255)

GAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCA‑C‑GTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311403‑4311582

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